Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,439,845 T→A 7.3% intergenic (+52/+27) ytfK → / ← ytfL DUF1107 family protein/UPF0053 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,439,8450TA 7.3% 310.0 / 11.7 105intergenic (+52/+27)ytfK/ytfLDUF1107 family protein/UPF0053 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (42/55);  new base A (4/4);  total (46/59)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

Only 100 of 108 total aligned reads displayed.

TAAGATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATCCACCACTT  >  NC_000913/4439697‑4439993
                                                                                                                                                    |                                                                                                                                                    
tAAGATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAt                                                                                                                                                      >  6:295499/1‑149 (MQ=255)
  aGATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTc                                                                                                                                                    <  8:273862/149‑1 (MQ=255)
   gATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAtt                                                                                                                                                     >  7:150664/1‑147 (MQ=255)
     ttCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAgg                                                                                                                                                 <  4:216427/149‑1 (MQ=255)
     ttCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAgg                                                                                                                                                 <  7:391982/149‑1 (MQ=255)
       cTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGa                                                                                                                                                <  3:65385/148‑1 (MQ=255)
          aTCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCg                                                                                                                                            >  7:23905/1‑149 (MQ=255)
           tCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGt                                                                                                                                           <  5:338592/149‑1 (MQ=255)
             ccAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGGCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGtc                                                                                                                                          <  7:84294/148‑2 (MQ=255)
             ccAAAAGGGAAGGACAAGCTGCATTTGTCAGTGATGGCCGACGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGGGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTg                                                                                                                                         <  4:29953/149‑1 (MQ=255)
              cAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGa                                                                                                                                        >  3:152626/1‑149 (MQ=255)
               aaaaGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGAt                                                                                                                                       <  6:83892/149‑1 (MQ=255)
               aaaaGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGAt                                                                                                                                       <  2:59985/149‑1 (MQ=255)
                  aGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGtcttatta                                                                                                                                    <  7:224435/149‑8 (MQ=255)
                    tGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTc                                                                                                                                  <  6:121756/149‑1 (MQ=255)
                        ggATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTggg                                                                                                                              <  1:77087/149‑1 (MQ=255)
                              gCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGTTTATGAGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGCTTCAGGCGCCGTTGATGGTTCTGGGGGTTAc                                                                                                                        >  4:32276/1‑149 (MQ=255)
                              gCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGcc                                                                                                                                             >  7:353410/1‑128 (MQ=255)
                              gCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGcc                                                                                                                                             <  8:353410/128‑1 (MQ=255)
                                tGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAcgc                                                                                                                      >  8:109470/1‑149 (MQ=255)
                                tGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAcgc                                                                                                                      >  1:217394/1‑149 (MQ=255)
                                  cATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGttt                                                                                                                                   >  5:106160/1‑134 (MQ=255)
                                  cATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGttt                                                                                                                                   <  6:106160/134‑1 (MQ=255)
                                  cATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAcgcg                                                                                                                     <  1:3449/148‑1 (MQ=255)
                                  cATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGa                                                                                                                    >  3:52155/1‑149 (MQ=255)
                                  cATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGa                                                                                                                    >  3:57429/1‑149 (MQ=255)
                                               tGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTcttctt                                                                                                        >  7:425192/1‑148 (MQ=255)
                                                 tCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttat                                                                                                     >  2:136439/1‑149 (MQ=255)
                                                 tCCGAAGTTAACCGGCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCGTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtcttta                                                                                                      <  5:36498/148‑1 (MQ=255)
                                                   cGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGCCGCTTTCTtctttatct                                                                                                   >  5:24843/1‑149 (MQ=255)
                                                   cGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTATCTtctttatct                                                                                                   >  8:237365/1‑149 (MQ=255)
                                                    gAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatctt                                                                                                  <  4:57429/149‑1 (MQ=255)
                                                    gAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatctt                                                                                                  >  3:9663/1‑149 (MQ=255)
                                                    gAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatctt                                                                                                  <  4:52155/149‑1 (MQ=255)
                                                      aGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttta                                                                                                <  2:217394/149‑1 (MQ=255)
                                                           aCCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGt                                                                                            <  7:214112/148‑1 (MQ=255)
                                                               tCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgg                                                                                                                               <  4:227647/109‑1 (MQ=255)
                                                               tCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgg                                                                                                                               >  3:227647/1‑109 (MQ=255)
                                                               tCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCGTTATCTTTAGCGTCAgg                                                                                        >  3:269715/1‑148 (MQ=255)
                                                                       tGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGGTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTgg                                                                               <  2:7976/149‑1 (MQ=255)
                                                                           tCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttta                                                                                                >  8:34239/1‑128 (MQ=255)
                                                                           tCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttta                                                                                                <  7:34239/128‑1 (MQ=255)
                                                                            cTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGc                                                                                       <  8:275996/136‑1 (MQ=255)
                                                                            cTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGc                                                                                       >  7:275996/1‑136 (MQ=255)
                                                                                  aCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCgg                                                                    >  6:165105/1‑149 (MQ=255)
                                                                                  aCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCgg                                                                    >  7:66723/1‑149 (MQ=255)
                                                                                     gagaTGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTg                                                                  <  1:413866/148‑1 (MQ=255)
                                                                                     gagaTGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTg                                                                  >  2:413866/1‑148 (MQ=255)
                                                                                       gaTGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGcc                                                               >  8:11228/1‑149 (MQ=255)
                                                                                         tGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCtt                                                             <  2:54833/149‑1 (MQ=255)
                                                                                           gCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTg                                                            <  8:425192/148‑1 (MQ=255)
                                                                                           gCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGc                                                           <  4:9663/149‑1 (MQ=255)
                                                                                             ttAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTg                                                         >  5:231819/1‑149 (MQ=255)
                                                                                             ttAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGGTTTGGCGAAAGGGCGGTGGCCTTGCTg                                                         >  4:46900/1‑149 (MQ=255)
                                                                                             ttAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGGTTTGGCGAAAGGGCGGTGGCCTTGCTg                                                         >  5:87547/1‑149 (MQ=255)
                                                                                               aaCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTGAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGt                                                        <  7:237365/148‑1 (MQ=255)
                                                                                                aCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAggg                                                                       <  1:234089/132‑1 (MQ=255)
                                                                                                aCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAggg                                                                       >  2:234089/1‑132 (MQ=255)
                                                                                                 ccAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGa                                                     >  5:12403/1‑149 (MQ=255)
                                                                                                      gtgctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCGGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCgg                                                <  7:109470/149‑1 (MQ=255)
                                                                                                      gtgctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGGTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCgg                                                <  6:112281/149‑1 (MQ=255)
                                                                                                      gggctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTTATTCACGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCgg                                                <  1:283848/147‑1 (MQ=255)
                                                                                                            tgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCAc                                           <  1:136439/148‑1 (MQ=255)
                                                                                                              cagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACca                                         <  8:23905/148‑1 (MQ=255)
                                                                                                               agtaaaaagtgctatgcagtaaTAAGACGGCTTCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGc                                                                              >  4:92588/1‑110 (MQ=255)
                                                                                                               agtaaaaagtgctatgcagtaaTAAGACGGCTTCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGc                                                                              <  3:92588/110‑1 (MQ=255)
                                                                                                               agtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACcagc                                       <  7:91243/149‑1 (MQ=255)
                                                                                                                 taaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACcagcag                                     >  7:377631/1‑149 (MQ=255)
                                                                                                                    aaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTg                                  >  7:136038/1‑149 (MQ=255)
                                                                                                                    aaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTGTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTg                                  <  5:147340/149‑1 (MQ=255)
                                                                                                                        tgctatgcagtaaTAAGACGGCTCCTGAATCAGGAGCGGTTGAGGTTTATGGGGGTTAAGCGACGCTTTCTTCTTTATCTTTAGAGTCAGGCAGGTTTGGCGAAAGAGCGGTGGCCTTGCGGTCGATCCGGGACACCAGCAGCTGGTct                              >  4:52471/1‑148 (MQ=255)
                                                                                                                          ctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGAt                            <  7:11228/149‑1 (MQ=255)
                                                                                                                           tatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttta                                                                                                >  3:58429/1‑80 (MQ=255)
                                                                                                                           tatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttta                                                                                                <  4:58429/80‑1 (MQ=255)
                                                                                                                           tatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGa                             >  5:373489/1‑147 (MQ=255)
                                                                                                                           tatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGa                             <  6:373489/147‑1 (MQ=255)
                                                                                                                            atgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGtt                                                                                   <  3:52132/92‑1 (MQ=255)
                                                                                                                              gcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGaaa                                                                          >  7:88080/1‑99 (MQ=255)
                                                                                                                              gcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGaaa                                                                          <  8:88080/99‑1 (MQ=255)
                                                                                                                               cagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCgg                        <  5:295499/148‑1 (MQ=255)
                                                                                                                                   aataaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATg                           >  7:329022/1‑141 (MQ=255)
                                                                                                                                   aataaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATg                           <  8:329022/141‑1 (MQ=255)
                                                                                                                                     taaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTAt                 >  4:243542/1‑149 (MQ=255)
                                                                                                                                      tcaaCGGCTCCTGAATCAGGAGCCGTCGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCgg                                                <  1:320002/114‑1 (MQ=255)
                                                                                                                                      tcaaCGGCTCCTGAATCAGGAGCCGTCGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCgg                                                >  2:320002/4‑117 (MQ=255)
                                                                                                                                       aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGtttt                                                                                 <  6:305551/83‑1 (MQ=255)
                                                                                                                                       aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATc                <  8:377631/148‑1 (MQ=255)
                                                                                                                                        gACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATc                <  8:232475/147‑1 (MQ=255)
                                                                                                                                        gACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATc                >  7:232475/1‑147 (MQ=255)
                                                                                                                                        gACGGCTCCTCATTCCGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGGCGATCCGGGGCACCAGCAGCTGGTCGATGCGGTAGTTATCGa              <  6:266327/149‑1 (MQ=255)
                                                                                                                                           ggCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGatat           <  3:243542/149‑1 (MQ=255)
                                                                                                                                               ccTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGa                             >  1:337257/1‑127 (MQ=255)
                                                                                                                                               ccTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGa                             <  2:337257/127‑1 (MQ=255)
                                                                                                                                                cTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATccac       >  1:296000/1‑148 (MQ=255)
                                                                                                                                                 tGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGGTGGTCGATGCGGTAGTTATCGATATgcgcca     >  6:318065/1‑149 (MQ=255)
                                                                                                                                                  gATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACcagcag                                     >  7:53691/1‑116 (MQ=255)
                                                                                                                                                  gATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACcagcag                                     <  8:53691/116‑1 (MQ=255)
                                                                                                                                                   aTTCGGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATCCACCACt   >  3:39561/1‑149 (MQ=255)
                                                                                                                                                   aTTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATccacca     <  6:12403/147‑1 (MQ=255)
                                                                                                                                                    ttCAGGAGGCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGTCCTTGCGGTCGATCCGGGTCACCAGCAGCTTGTCGATGCGGTAGTTATCGATATCCACCACtt  <  3:32276/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAAGATTCTTATCCCAAAAGTGAAGGATAAGCTGCATTTGTCAGTGATGTCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATCCGGGTCACCAGCAGCTGGTCGATGCGGTAGTTATCGATATCCACCACTT  >  NC_000913/4439697‑4439993

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: