Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,450:1 +A 5.6% intergenic (+22/+57) 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,4501.A 5.6% 762.8 / 39.6 180intergenic (+22/+57)ftnB/yecJferritin B, putative ferrous iron reservoir/DUF2766 family protein
Reads supporting (aligned to +/- strand):  ref base . (80/90);  new base A (4/6);  total (84/96)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

Only 100 of 187 total aligned reads displayed.

CAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAAC  >  NC_000913/1987186‑1987713
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                       
                  tACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCAACCCAGCTGCATTGATCATCATCGGAGCTAATGCATT‑GCGCCGGTGAAGGTTTTGAGAACCCGCt                                                                                                                                                                                                                                             >  2:35822/1‑275 (MQ=255)
                                        tGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCACTTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGACGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGCTCATAACCGGCTCTAATGCATT‑GCGCCGATCAAGGTTTTGAGAACCCGCTGCCTCATCTGTTTGAagccgc                                                                                                                                                                                                                        >  4:146405/1‑274 (MQ=255)
                                                        tGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCCCACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTTAGAAACCGCTGCCTCCTCTGTTTGAAGCAGCGGTTTTTTTATTGGGAt                                                                                                                                                                                                       >  1:10303/1‑275 (MQ=255)
                                                                   aaCTGAATGATTATTCAACCGTCAATTTCCCGCGCGATCTGGCAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtgtg                                                                                                                                                                                             <  1:48421/274‑1 (MQ=255)
                                                                            aggatTGAACCTTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCGTCGGCGCTAATGCATT‑GCGCCGATGAAGTTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACt                                                                                                                                                                                    <  4:32161/272‑1 (MQ=255)
                                                                                              tCCTGCGCGATCTGGAAAAAGAACAGCAGCATGAGGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGa                                                                                                                                                                  <  2:99387/274‑1 (MQ=255)
                                                                                                         cTGGAAAAAGAACAGCAGCATGAGTGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGGGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGa                                                                                                                                                      <  2:273249/275‑1 (MQ=255)
                                                                                                                    aCAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCAtt                                                                                                                                                                                                                                                                          >  5:339238/1‑149 (MQ=255)
                                                                                                                      agcagcATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑gc                                                                                                                                                                                                                                                                       >  5:1479064/1‑149 (MQ=255)
                                                                                                                       gcagcaTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑gcg                                                                                                                                                                                                                                                                      >  5:1157340/1‑149 (MQ=255)
                                                                                                                          gcaTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGC‑AATGCATTAGCGCCg                                                                                                                                                                                                                                                                   <  8:971645/149‑1 (MQ=255)
                                                                                                                           caTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGa                                                                                                                                                                                                                                                                  <  6:339238/149‑1 (MQ=255)
                                                                                                                           caTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGC‑AATGCATTAGCGCCGa                                                                                                                                                                                                                                                                  <  7:719143/149‑1 (MQ=255)
                                                                                                                            atgatgGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGAGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGC‑AATGCATTAGCGCCGAt                                                                                                                                                                                                                                                                 <  5:1461904/149‑1 (MQ=255)
                                                                                                                                  gTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGt                                                                                                                                                                                                                                                           >  7:322845/1‑149 (MQ=255)
                                                                                                                                    cTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGttt                                                                                                                                                                                                                                                         <  8:322845/149‑1 (MQ=255)
                                                                                                                                     tGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGtttt                                                                                                                                                                                                                                                        <  5:951761/149‑1 (MQ=255)
                                                                                                                                     tGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGtttt                                                                                                                                                                                                                                                        <  8:165842/149‑1 (MQ=255)
                                                                                                                                           tgcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGaaa                                                                                                                                                                                                                                                  >  5:363214/1‑149 (MQ=255)
                                                                                                                                               aaCCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGc                                                                                                                                                                                                                                              <  6:58980/149‑1 (MQ=255)
                                                                                                                                                aCCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCt                                                                                                                                                                                                                                             <  6:1401835/149‑1 (MQ=255)
                                                                                                                                                  cATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGc                                                                                                                                                                                                                                           <  8:948230/149‑1 (MQ=255)
                                                                                                                                                   aTTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGcc                                                                                                                                                                                                                                          >  6:346043/1‑149 (MQ=255)
                                                                                                                                                   aTTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGcc                                                                                                                                                                                                                                          >  8:677494/1‑149 (MQ=255)
                                                                                                                                                       ttGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGttt                                                                                                                                                                                                                                                         <  7:1326379/130‑1 (MQ=255)
                                                                                                                                                       ttGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGttt                                                                                                                                                                                                                                                         >  8:1326379/1‑130 (MQ=255)
                                                                                                                                                         gatgaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCt                                                                                                                                                                                                                                    <  7:1240989/149‑1 (MQ=255)
                                                                                                                                                          atgaAGTGCGGAGTGCGAAACTTGCGGGTATGTGACCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGACGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGt                                                                                                      <  2:149586/274‑1 (MQ=255)
                                                                                                                                                            gaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGtt                                                                                                                                                                                                                                 <  5:192413/149‑1 (MQ=255)
                                                                                                                                                               gTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGatat                                                                                                >  2:45790/1‑275 (MQ=255)
                                                                                                                                                                 gcgcAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAag                                                                                                                                                                                                                            <  8:389353/149‑1 (MQ=255)
                                                                                                                                                                 gcgcAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAag                                                                                                                                                                                                                            <  7:1184336/149‑1 (MQ=255)
                                                                                                                                                                 gcgcAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAag                                                                                                                                                                                                                            >  7:389353/1‑149 (MQ=255)
                                                                                                                                                                         cGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtt                                                                                                                                                                                                                    <  6:363214/149‑1 (MQ=255)
                                                                                                                                                                         cGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGcc                                                                                                                                                 <  2:207368/216‑1 (MQ=255)
                                                                                                                                                                          gAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttt                                                                                                                                                                                                                   <  6:27518/149‑1 (MQ=255)
                                                                                                                                                                           aaaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttt                                                                                                                                                                                                                  >  7:1084610/1‑149 (MQ=255)
                                                                                                                                                                               ttGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTa                                                                                                                                                                                                              >  5:1358095/1‑149 (MQ=255)
                                                                                                                                                                                tGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTaa                                                                                                                                                                                                             <  7:785250/149‑1 (MQ=255)
                                                                                                                                                                                 gCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATggg                                                                                                                                                                                                         >  2:364192/1‑152 (MQ=255)
                                                                                                                                                                                     gTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGa                                                                                                                                                                                                        <  6:47472/149‑1 (MQ=255)
                                                                                                                                                                                         gtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCa                                                                                                                                                                                                    >  6:749226/1‑149 (MQ=255)
                                                                                                                                                                                              cTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtg                                                                                                                                                                                               >  8:414364/1‑149 (MQ=255)
                                                                                                                                                                                                gtgCAGACCGACCAACATGTTCTGAATGTCGTGTCACACGAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtgtg                                                                                                                                                                                             <  7:960150/149‑1 (MQ=255)
                                                                                                                                                                                                gtgCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtgtg                                                                                                                                                                                             >  5:22688/1‑149 (MQ=255)
                                                                                                                                                                                                 tgCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTgg                                                                                                                                                                                            <  5:247436/149‑1 (MQ=255)
                                                                                                                                                                                                  gCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTggg                                                                                                                                                                                           <  8:1360107/149‑1 (MQ=255)
                                                                                                                                                                                                  gCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTggg                                                                                                                                                                                           >  6:787595/1‑149 (MQ=255)
                                                                                                                                                                                                      accgaccAACATGTTCTGAATGTCGTGTCACACCAGCTTCATTGATCATCATCGGCGC‑AATGCATTAGCGCCGATGatg                                                                                                                                                                                                                                                             >  6:1020458/1‑77 (MQ=25)
                                                                                                                                                                                                      accgaccAACATGTTCTGAATGTCGTGTCACACCAGCTTCATTGATCATCATCGGCGC‑AATGCATTAGCGCCGATGatg                                                                                                                                                                                                                                                             <  5:1020458/79‑3 (MQ=25)
                                                                                                                                                                                                           ccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAcc                                                                                                                                                                                                  <  8:884577/133‑1 (MQ=255)
                                                                                                                                                                                                           ccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAAAGGgcggg                                                    >  2:123747/1‑273 (MQ=255)
                                                                                                                                                                                                           ccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGAGGGCTTTTTGGATTGCCCGCATGGTGACCGGGTAGGCGGCGGCGGGGTGGGTAGTTAAGTCGAGACGCTTCAGCTGACTGGACCTTTGTTCACGAAGTTCAACgggggac                                                    >  4:245512/1‑271 (MQ=255)
                                                                                                                                                                                                            cAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAg                                                                                                                                                                                                                                                             >  7:840393/1‑73 (MQ=255)
                                                                                                                                                                                                             aaCATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAg                                                                                                                                                                                >  7:604789/1‑149 (MQ=255)
                                                                                                                                                                                                             aaCATGTTCTGAATGTCGTGGCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAg                                                                                                                                                                                <  8:855890/149‑1 (MQ=255)
                                                                                                                                                                                                               cATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGtt                                                                                                                                                                              >  5:886477/1‑149 (MQ=255)
                                                                                                                                                                                                                aTGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa                                                                                                                                                                                 >  7:1012209/1‑145 (MQ=255)
                                                                                                                                                                                                                aTGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa                                                                                                                                                                                 <  8:1012209/145‑1 (MQ=255)
                                                                                                                                                                                                                   ttCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                                                                                                                                                            <  6:608871/115‑1 (MQ=255)
                                                                                                                                                                                                                   ttCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                                                                                                                                                            >  5:608871/1‑115 (MQ=255)
                                                                                                                                                                                                                   ttCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAAt                                                                                                                                                                          >  7:1328223/1‑149 (MQ=255)
                                                                                                                                                                                                                   ttCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAAt                                                                                                                                                                          >  5:856868/1‑149 (MQ=255)
                                                                                                                                                                                                                   ttCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAGGGGATTCACCCTGTGGGGTAAACTTGa                                                                                                                                                                                 >  7:873461/1‑142 (MQ=255)
                                                                                                                                                                                                                       gAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttttt                                                                                                                                                                                                                >  7:1347821/1‑107 (MQ=255)
                                                                                                                                                                                                                       gAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttttt                                                                                                                                                                                                                <  8:1347821/107‑1 (MQ=255)
                                                                                                                                                                                                                            tCGTGTCAAACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATc                                   <  4:45790/275‑1 (MQ=255)
                                                                                                                                                                                                                             cGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTgg                                                                                                                                                                                            >  4:359701/1‑121 (MQ=255)
                                                                                                                                                                                                                                   acacCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCtttt                                                                                                                                                          >  7:485642/1‑149 (MQ=255)
                                                                                                                                                                                                                                   acacCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCtttt                                                                                                                                                          >  7:211254/1‑149 (MQ=255)
                                                                                                                                                                                                                                    cacCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCttttt                                                                                                                                                         >  5:1255279/1‑149 (MQ=255)
                                                                                                                                                                                                                                    cacCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCttttt                                                                                                                                                         <  8:621630/149‑1 (MQ=255)
                                                                                                                                                                                                                                    cacCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCttttt                                                                                                                                                         <  5:305613/149‑1 (MQ=255)
                                                                                                                                                                                                                                    cacCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCttttt                                                                                                                                                         <  7:1435822/149‑1 (MQ=255)
                                                                                                                                                                                                                                     acCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGt                                                                                                                                                                                                                                  <  8:1223317/75‑1 (MQ=255)
                                                                                                                                                                                                                                     acCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTg                                                                                                                                                        <  7:299721/149‑1 (MQ=255)
                                                                                                                                                                                                                                     acCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTg                                                                                                                                                        <  8:448496/149‑1 (MQ=255)
                                                                                                                                                                                                                                     acCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTg                                                                                                                                                        <  6:123778/149‑1 (MQ=255)
                                                                                                                                                                                                                                       cAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCACCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGa                                                                                                                                                      >  8:732207/1‑149 (MQ=255)
                                                                                                                                                                                                                                           tGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGc                                                                                                                                                  <  7:677494/149‑1 (MQ=255)
                                                                                                                                                                                                                                             cATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAgcgc                                                                                                                                                                    >  7:11208/1‑129 (MQ=255)
                                                                                                                                                                                                                                             cATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAgcgc                                                                                                                                                                    <  8:11208/129‑1 (MQ=255)
                                                                                                                                                                                                                                                  atcatcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATg                                                                                                                                           >  5:645057/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    catcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGt                                                                                                                                         <  7:179425/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    catcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGt                                                                                                                                         >  8:258882/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    catcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTGTAATGGGATTCTCCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTg                                                                                                                                                        >  6:1053632/1‑134 (MQ=255)
                                                                                                                                                                                                                                                     atcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTg                                                                                                                                        >  6:238943/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     atcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTg                                                                                                                                        >  6:682212/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       catcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGAc                                                                                                                                      >  7:1196937/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         tcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGa                                                                                                                                                      >  5:223511/1‑131 (MQ=255)
                                                                                                                                                                                                                                                         tcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGa                                                                                                                                                      <  6:223511/131‑1 (MQ=255)
                                                                                                                                                                                                                                                         tcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCg                                                                                                                                    <  8:604789/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         tcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCg                                                                                                                                    >  7:925408/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         tcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                                                                                                                                                            >  6:1475484/1‑77 (MQ=38)
                                                                                                                                                                                                                                                          cGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCgg                                                                                                                                   <  7:414364/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          cGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCgg                                                                                                                                   <  5:542007/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             cgcTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCACCGGGCGCGATTACATCGAGACAATCCATAATTTGTTTGATAACCAACTGGCAAGCAAc  >  4:114820/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                tAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgc                                                                                                                             <  5:682212/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aaTGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgcg                                                                                                                            <  5:442943/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tt‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgcggcggcg                                                                                                                      >  5:54434/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                       
CAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAAC  >  NC_000913/1987186‑1987713

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