Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,600,390 G→A Q132* (CAG→TAG)  rpoH ← RNA polymerase sigma factor RpoH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,600,3900GA93.3% 141.8 / ‑0.7 46Q132* (CAG→TAG) rpoHRNA polymerase sigma factor RpoH
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (17/25);  total (19/27)
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.50e-01

CGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTGCGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGTGAGACAGGATCAGCG  >  NC_000913/3600153‑3600629
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
cGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTc                                                                                                                                                                                                                                    >  1:59128/1‑251 (MQ=255)
 gTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTGCGCTTTGGTGGTCg                                                                                                                                                                                                                                   >  2:53991/1‑251 (MQ=255)
                    ttATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGcc                                                                                                                                                                                                                 <  1:95881/250‑1 (MQ=255)
                        ccTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCTTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCCTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTa                                                                                                                                                                                                            <  1:118556/251‑1 (MQ=255)
                        ccTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGtt                                                                                                                                                                                                             >  1:183597/1‑250 (MQ=255)
                           gcacaTAGAGCACCGGAGCCATCGGGTGGCTGTCGGACTCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCGGTTCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCTCCATTTCGACTTCATCCTGGTTAAACCAGCCCATACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTAc                                                                                                                                                                                                           <  2:12791/245‑1 (MQ=255)
                                 agagCACCGGAGCCATCGGCTGGCTGTCGGAATCGTCGTCGGACGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACg                                                                                                                                                                                                   >  1:139461/1‑251 (MQ=255)
                                  gagCACCGGAGCCATCGGCTGGCTGTCGGAATCGACGTCGGAAGACAGGTAAAAGGTCAAGTCCTGTGCCGCCATACGTGATTCCATCTCACGCCCGACTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTAAGCTTTGGTGGTCGCAACTTTGACGATACGCCAATTACGAAGAACg                                                                                                                                                                                                   >  1:118586/1‑250 (MQ=255)
                                    gCACCGGAGCCATCGGCTGGCTGTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCATAACGTAt                                                                                                                                                                                                <  1:193378/251‑1 (MQ=255)
                                      aCCGGAGCCATCGGCTGGCTGTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTc                                                                                                                                                                                              >  1:36455/1‑251 (MQ=255)
                                       ccgtaGCCATCGGCTGGCTGTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCg                                                                                                                                                                                             <  1:149062/247‑1 (MQ=255)
                                        cGGAGCCATCGGCTGGCTGTCGGACTCGCCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCAGCCATACGAGATTCCATATCACGTCCGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGCAGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTCCGCTTTGGTGGTCGCACCTTTGCCGCTACGCCAGTTACGCAGAACGTATTCGt                                                                                                                                                                                            >  2:233120/1‑251 (MQ=255)
                                         ggAGCCATCGGCTGGCTGTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTg                                                                                                                                                                                           <  2:59128/251‑1 (MQ=255)
                                                     gggctgGCGGAATCGTCGGCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGGTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTGGCGCTGCGCTTTGGGGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATGCGTGGTTCTCTGCttt                                                                                                                                                                               <  2:259193/250‑1 (MQ=255)
                                                      ggctgTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTGTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCttt                                                                                                                                                                               <  2:277589/250‑1 (MQ=255)
                                                            cGGAATCGTCGGCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCGCTGCTTTGATCCAg                                                                                                                                                                        <  2:43817/251‑1 (MQ=255)
                                                                   gtcgtcgGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACgg                                                                                                                                                                 >  2:110058/1‑251 (MQ=255)
                                                                     cgtcgGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATGCGTGGATCTCTGCTTTGATCCAGTGAACGGc                                                                                                                                                                <  1:23054/250‑1 (MQ=255)
                                                                             gACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGAc                                                                                                                                                       <  1:178148/251‑1 (MQ=255)
                                                                               aaGGTCAAAGTGCATGTCCTGTGCCGCCATACGTGAGTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTGTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATGCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCa                                                                                                                                                     <  2:64671/250‑1 (MQ=255)
                                                                                  gTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGgc                                                                                                                                                  >  2:142036/1‑251 (MQ=255)
                                                                                      aaGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGcac                                                                                                                                              >  1:217610/1‑251 (MQ=255)
                                                                                       aGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGcaca                                                                                                                                             <  1:19678/251‑1 (MQ=255)
                                                                                       aGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGcac                                                                                                                                              >  1:254311/1‑250 (MQ=255)
                                                                                          tCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACAccc                                                                                                                                          >  2:18326/1‑251 (MQ=255)
                                                                                               tCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCCGTTCACGGGCCATCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGATCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTc                                                                                                                                     >  2:28846/1‑251 (MQ=255)
                                                                                                ccTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCCTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTc                                                                                                                                     <  1:39654/250‑1 (MQ=255)
                                                                                                ccTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTCTGCTGGTTCCGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGCGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTcc                                                                                                                                    >  1:177004/1‑251 (MQ=255)
                                                                                                ccTGGGCCGCCATTCGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCATGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGACCGTATGCGTGGATCTCTGCTTTGATCCAGTGAACGTCGAAGTAGACCAGGCGCACACCCTCTTcc                                                                                                                                    <  2:273756/251‑1 (MQ=255)
                                                                                                                 gATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCAc                                                                                                                   >  2:245646/1‑251 (MQ=255)
                                                                                                                     ccATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGc                                                                                                                <  2:42312/250‑1 (MQ=255)
                                                                                                                           cACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGGCCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTtc                                                                                                            <  1:228471/248‑1 (MQ=255)
                                                                                                                             cgtacgtCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTGCGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAg                                                                                                       <  1:53991/251‑1 (MQ=255)
                                                                                                                                acgtCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGTAGACCAGGCGCACACCCGCTTCCGGTTTGAAGCGGCGCACTGCTTTCATCAGGcc                                                                                                    <  2:139461/251‑1 (MQ=255)
                                                                                                                                 cgtCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGACCg                                                                                                   >  1:176764/1‑251 (MQ=255)
                                                                                                                                      ttGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATg                                                                                                <  2:224083/249‑1 (MQ=255)
                                                                                                                                       tGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGtt                                                                                              <  1:175374/250‑1 (MQ=255)
                                                                                                                                                  cccAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAAt                                                                                  >  2:123890/1‑251 (MQ=255)
                                                                                                                                                   ccAGTTCACGGGCCACCATAGCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACGGTTTGCGCTACGCTTTGTTGGTCGCAACTTTGACGATACTTCAGTTACGCAGAACGTATTCGTGGATCTCTGCTGTGATCCAGTGAACGGCGAAGGAGACCATGCGCACACCCTCTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGGCGATGTTACCTTCCTGAATc                                                                                 <  2:246173/251‑1 (MQ=255)
                                                                                                                                                                   aaTTTCGACTTCATCCGGGTTAAACCAGCCCAGACGCTGTTGGGTTTTACGCAGTTTAAAGAACAGTTTGCGCTACGCTTTTGTGGTCGCAACTTGGGCGATCCGCCTGTTAGGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACAAGGCGCACACCCACTCCCGGGTTGAATCGGCGCACTGCTTTCATCGGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCa                                                                 <  2:255953/250‑1 (MQ=255)
                                                                                                                                                                        cGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGcc                                                             >  1:123918/1‑250 (MQ=255)
                                                                                                                                                                                        aaaCCAGCCCAGACGCTGTTTGGTTTTACGCAGGTGGAAGAACAGTTTGGGCGACGTTTTGTTGGTCGTAACTTTGACGATAAGCCAGTTACGCAGAACGTATTCGTGGATCGCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCTGGTTGAAACGTCGCTCTGCTTGCATCAGGCCGATGTTACCTTCCTGATTCAAATCCGCCTGTGGCAGGCCAGAGCCCGCATAAttt                                             <  2:151762/250‑2 (MQ=255)
                                                                                                                                                                                          aCCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGa                                          <  1:183776/251‑1 (MQ=255)
                                                                                                                                                                                                                    cGCAGGTTGAAGAACAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGTg                <  1:91022/251‑1 (MQ=255)
                                                                                                                                                                                                                                aaCAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGTGAGACAGGATCAg    <  1:134330/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cAGTTTGCGCTACGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACACACCGCAGGTGAGACCGGATCAGCg  >  1:191957/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cgctacgctTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATcc                                                                           >  2:109512/1‑171 (MQ=255)
                                                                                                                                                                                                                                         cgctacgctTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATcc                                                                           <  1:109512/171‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
CGTCGGCAAAGTTAGATGATTTATCCTGCAGATAGAGCACCGGAGCCATCGGCTGGCTGTCGGAATCGTCGTCGGAAGACAGGTCAAAGGTCATGTCCTGTGCCGCCATACGTGATTCCATCTCACGTACGTCTTTGCTGGTTACGCCCAGTTCACGGGCCACCATTTCGACTTCATCCTGGTTAAACCAGCCCAGACGCTGCTTGGTTTTACGCAGGTTGAAGAACAGTTTGCGCTGCGCTTTGGTGGTCGCAACTTTGACGATACGCCAGTTACGCAGAACGTATTCGTGGATCTCTGCTTTGATCCAGTGAACGGCGAAGGAGACCAGGCGCACACCCACTTCCGGGTTGAAACGGCGCACTGCTTTCATCAGGCCGATGTTACCTTCCTGAATCAAATCCGCCTGTGGCAGGCCATAGCCCGCATAATTACGAGCAATATGAACAACAAACCGCAGGTGAGACAGGATCAGCG  >  NC_000913/3600153‑3600629

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: