Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,442 Δ1 bp 6.1% intergenic (+14/+65) ftnB → / ← yecJ ferritin B, putative ferrous iron reservoir/DUF2766 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,4420T. 6.1% 292.5 / 11.3 82intergenic (+14/+65)ftnB/yecJferritin B, putative ferrous iron reservoir/DUF2766 family protein
Reads supporting (aligned to +/- strand):  ref base T (33/44);  new base . (2/3);  total (35/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

AAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGC  >  NC_000913/1987296‑1987590
                                                                                                                                                  |                                                                                                                                                     
aaaaGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa                                                                                                                                                     <  4:891347/149‑1 (MQ=255)
aaaaGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa                                                                                                                                                     >  4:867913/1‑149 (MQ=255)
 aaaGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAAt                                                                                                                                                    >  2:41070/1‑149 (MQ=255)
     aaCAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCAt                                                                                                                                                >  4:37961/1‑149 (MQ=255)
      tttgcagcaTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCAtt                                                                                                                                               >  2:639351/4‑149 (MQ=255)
       cagcagCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑g                                                                                                                                             <  2:468849/149‑1 (MQ=255)
       cagcagCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑g                                                                                                                                             <  3:233499/149‑1 (MQ=255)
          cagcaTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑gcgc                                                                                                                                          <  4:433600/149‑1 (MQ=255)
          cagcaTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGC‑AATGCAttagcgc                                                                                                                                          <  1:949327/149‑6 (MQ=255)
              atgatgGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGAt                                                                                                                                      <  3:65692/149‑1 (MQ=255)
              atgatgGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGAt                                                                                                                                      <  3:37961/149‑1 (MQ=255)
                  tgGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGC‑AATGCATTAGCGCCGATGatg                                                                                                                                  <  1:506739/149‑3 (MQ=255)
                     tCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGtt                                                                                                                               <  2:1068791/149‑1 (MQ=255)
                       tGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGtttt                                                                                                                             <  3:1121577/149‑1 (MQ=255)
                       tGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGtttt                                                                                                                             <  2:989355/149‑1 (MQ=255)
                             tgcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGaaa                                                                                                                       >  1:1131476/1‑149 (MQ=255)
                             tacAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGaaa                                                                                                                       <  3:935902/147‑1 (MQ=255)
                              gcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAAc                                                                                                                      <  2:972888/149‑1 (MQ=255)
                              gcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAAc                                                                                                                      >  4:430845/1‑149 (MQ=255)
                                 aaCCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGGCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGc                                                                                                                   >  2:66594/1‑149 (MQ=255)
                                 aaCCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGc                                                                                                                   >  2:1069392/1‑149 (MQ=255)
                                  aCCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCt                                                                                                                  <  1:1069392/149‑1 (MQ=255)
                                  aCCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCt                                                                                                                  <  4:238159/149‑1 (MQ=255)
                                    cATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGc                                                                                                                >  3:313225/1‑149 (MQ=255)
                                         ttGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAAt                                                                                                                                                    <  2:166101/109‑1 (MQ=255)
                                         ttGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAAt                                                                                                                                                    >  1:166101/1‑109 (MQ=255)
                                         ttGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCAt                                                                                                           >  3:502917/1‑149 (MQ=255)
                                           gatgaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCTTCaaac                                                                                                         >  3:159391/1‑146 (MQ=255)
                                           gatgaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCt                                                                                                         <  2:316113/149‑1 (MQ=255)
                                              gaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGtt                                                                                                      <  4:966872/149‑1 (MQ=255)
                                               aaGTGCGCAGTGCGAAACTTGCGGGTATGTGGCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGttt                                                                                                     <  1:66594/149‑1 (MQ=255)
                                                  tGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGaa                                                                                                  <  4:505084/149‑1 (MQ=255)
                                                   gcgcAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAag                                                                                                 <  3:867913/149‑1 (MQ=255)
                                                     gcAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAagca                                                                                               >  3:818718/1‑149 (MQ=255)
                                                       aGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAagcagc                                                                                             <  2:1131476/149‑1 (MQ=255)
                                                            gAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGtttt                                                                                                                             >  1:689872/1‑112 (MQ=255)
                                                            gAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGtttt                                                                                                                             <  2:689872/112‑1 (MQ=255)
                                                              aaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttttt                                                                                      >  2:453606/1‑149 (MQ=255)
                                                                 ttGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTa                                                                                   >  3:490131/1‑149 (MQ=255)
                                                                  tGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTaa                                                                                  <  4:313225/149‑1 (MQ=255)
                                                                  tGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTaa                                                                                  <  2:930470/149‑1 (MQ=255)
                                                                     gggTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGTGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATgg                                                                               >  4:751982/1‑149 (MQ=255)
                                                                      ggTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTaa                                                                                  <  1:526064/145‑1 (MQ=255)
                                                                      ggTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTaa                                                                                  >  2:526064/1‑145 (MQ=255)
                                                                           gtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCa                                                                         >  1:581697/1‑149 (MQ=255)
                                                                             gCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAcc                                                                       <  3:289826/149‑1 (MQ=255)
                                                                                cTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtg                                                                    >  3:103473/1‑149 (MQ=255)
                                                                                   tgCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtgtg                                                                  <  2:162697/148‑1 (MQ=255)
                                                                                   tgCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtgtg                                                                  >  1:162697/1‑148 (MQ=255)
                                                                                      aGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGttt                                                                                                     <  2:203094/110‑1 (MQ=255)
                                                                                      aGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGttt                                                                                                     >  1:203094/1‑110 (MQ=255)
                                                                                      aGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGt                                                              <  1:639351/149‑1 (MQ=255)
                                                                                       gaccgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttttt                                                                                     >  2:491704/1‑125 (MQ=255)
                                                                                       gaccgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCa                                                                         >  3:287449/1‑137 (MQ=255)
                                                                                       gaccgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCa                                                                         <  4:287449/137‑1 (MQ=255)
                                                                                       gaccgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGAGCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttttt                                                                                     <  1:491704/125‑1 (MQ=255)
                                                                                              cAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa                                                      <  3:751982/149‑1 (MQ=255)
                                                                                               aaCATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttttt                                                                                     <  3:1137228/117‑1 (MQ=255)
                                                                                               aaCATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttttt                                                                                     >  4:1137228/1‑117 (MQ=255)
                                                                                                aCATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa                                                      <  3:176271/147‑1 (MQ=255)
                                                                                                aCATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa                                                      >  4:176271/1‑147 (MQ=255)
                                                                                                 cATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                                 <  3:264385/119‑1 (MQ=255)
                                                                                                 cATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                                 >  4:264385/1‑119 (MQ=255)
                                                                                                  aTGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATCCTCCCTGTGGGGTAACCTTGAGtac                                                  <  4:490131/149‑3 (MQ=255)
                                                                                                     ttCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa                                                      >  3:422410/1‑142 (MQ=255)
                                                                                                     ttCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa                                                      <  4:422410/142‑1 (MQ=255)
                                                                                                     ttCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAAt                                               >  2:42761/1‑149 (MQ=255)
                                                                                                       cTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATaa                                             <  4:502917/149‑1 (MQ=255)
                                                                                                              tCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGAt                                      <  4:103473/149‑1 (MQ=255)
                                                                                                                     acacCAGCTGCATTGAACATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAc                                                                        >  4:64196/1‑108 (MQ=255)
                                                                                                                     acacCAGCTGCATTGAACATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAc                                                                        <  3:64196/108‑1 (MQ=255)
                                                                                                                      cacCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCttttt                              <  4:818718/149‑1 (MQ=255)
                                                                                                                      cacCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCttttt                              <  1:1086775/149‑1 (MQ=255)
                                                                                                                       acCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTg                             <  1:42761/149‑1 (MQ=255)
                                                                                                                         cAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGa                           <  3:430845/149‑1 (MQ=255)
                                                                                                                                    atcatcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCa                  >  1:510093/1‑147 (MQ=255)
                                                                                                                                    atcatcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCa                  <  2:510093/147‑1 (MQ=255)
                                                                                                                                       ttcatcGGCGC‑AATGCATTAGCGCCGATGATGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCttttt                              >  1:764170/2‑132 (MQ=39)
                                                                                                                                       ttcatcGGCGC‑AATGCATTAGCGCCGATGATGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCttttt                              <  2:764170/131‑1 (MQ=39)
                                                                                                                                           tcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCg         >  2:271663/1‑149 (MQ=255)
                                                                                                                                                gcTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCt    >  1:136156/1‑149 (MQ=255)
                                                                                                                                                gcTAATGCATT‑GCGACGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGCTTCACCCTGTGGGGTAAAATTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGACCGCATGGTGACCGGGTCt    >  4:263033/1‑149 (MQ=255)
                                                                                                                                                  tAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgc  >  1:1067345/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGC  >  NC_000913/1987296‑1987590

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