Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,228,227 G→C 16.2% pseudogene (381/1287 nt) mdtQ ← pseudogene, putative channel/filament proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2270GC16.2% 124.4 / 2.8 59pseudogene (381/1287 nt)mdtQpseudogene, putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base G (19/28);  new base C (8/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 3.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

TTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAG  >  NC_000913/2227994‑2228459
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
ttcattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTc                                                                                                                                                                                                                                          <  2:197293/239‑1 (MQ=255)
        ccccGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGccc                                                                                                                                                                                                                                  <  1:31047/239‑1 (MQ=255)
         cccGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGACAGCCAGCAATTGGCGTGTTTGCTCGCCGTTCCCCGccc                                                                                                                                                                                                                                  >  1:1389657/1‑238 (MQ=255)
            gcgACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTACCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGc                                                                                                                                                                                                                              <  1:443207/239‑1 (MQ=255)
              gACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCtt                                                                                                                                                                                                                            <  2:105326/239‑1 (MQ=255)
                 aTCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTaa                                                                                                                                                                                                                         <  2:789846/239‑1 (MQ=255)
                 aTCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTaa                                                                                                                                                                                                                         <  2:548242/239‑1 (MQ=255)
                 aTCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGACGCCCGTGCTTTaa                                                                                                                                                                                                                         <  1:151826/239‑1 (MQ=255)
                         gctgctGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTAcc                                                                                                                                                                                                                  <  2:1334211/238‑1 (MQ=255)
                          ctgctgCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTTCCCCGTTTCCCCCCCCGGCGTTTAACCGCAcccc                                                                                                                                                                                                                >  1:418316/1‑238 (MQ=255)
                               gcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCgg                                                                                                                                                                                                           <  2:744170/239‑1 (MQ=255)
                               gcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCCCCCCCCCGGCTTTTACCGGAACCcggcgg                                                                                                                                                                                                           >  1:322603/1‑239 (MQ=255)
                                  gggTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCtt                                                                                                                                                                                                                            <  2:1212908/219‑1 (MQ=255)
                                  gggTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCtt                                                                                                                                                                                                                            >  1:1212908/1‑219 (MQ=255)
                                   ggTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGc                                                                                                                                                                                                             <  1:894355/233‑1 (MQ=255)
                                     ttttGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCGAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGt                                                                                                                                                                                                      <  2:126143/238‑1 (MQ=255)
                                         gCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCt                                                                                                                                                                                                 <  1:316942/239‑1 (MQ=255)
                                             gCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCCGTTCGCCGCCCCGGGTTTTACCGACCCCCGGCGGGGAGtcaacccct                                                                                                                                                                                              >  2:894355/1‑230 (MQ=255)
                                               aTTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGccc                                                                                                                                                                                           <  1:670773/239‑1 (MQ=255)
                                                     atCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCGGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCt                                                                                                                                                                                     <  2:828557/239‑1 (MQ=255)
                                                         gTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCTCGTTCCCCCGCCCGTGCTTTAACCGTCCCCCGGCGGGCAGGAAACTCCGCCCCCTTCtttt                                                                                                                                                                                  >  1:643030/1‑237 (MQ=255)
                                                           ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTcccc                                                                                                                                                                               <  1:425302/239‑1 (MQ=255)
                                                           ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTccc                                                                                                                                                                                <  1:12253/238‑1 (MQ=255)
                                                           ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTccc                                                                                                                                                                                <  1:827047/238‑1 (MQ=255)
                                                                        aCTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTCCGCCCCCCCTGCCTTTACCCCACCCAGGCGGGCAGGAACCTCCCCCCCCTTCTTTcccccagtatcgcgaa                                                                                                                                                                   >  2:250724/1‑226 (MQ=255)
                                                                            aagaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCCTGCTTTTACCCGACCCAGGGGGGCAGTAACCTCCCCCCCGTTCTTTCCCCCGATATCGAGATGccc                                                                                                                                                              >  1:503147/1‑237 (MQ=255)
                                                                             agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGGTCGCGGTCCGCCCCCCGGGCTTTTACCGGACCCCGGCGGGCAGGAAACCCCCCCCGGTTTTTTTCCCCGATATCGAGATACCAg                                                                                                                                                             >  1:292781/1‑238 (MQ=255)
                                                                             agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCCCCGCCCGTGCTTTTACCGTACCCCGGCCGGCAGTAACCTCCCCCCCCTTCTTTCCCCCAATATCGCGATGcccg                                                                                                                                                             >  2:1336744/1‑236 (MQ=255)
                                                                             agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCCCGTTCCCCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTTCCCAGATATCGAGATGCCAg                                                                                                                                                             >  1:1066356/1‑238 (MQ=255)
                                                                              gaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCCCGTTCCGCCCCCCCTGCTTTTACCGTACCCAGGCGGGCAGTAACCCCCCCCCGGTTCTTTTCCCCGAAAACGGGCGGCCAGc                                                                                                                                                           >  1:870078/1‑238 (MQ=255)
                                                                                 gTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGGACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCATATATCGAGATGCCAGCccg                                                                                                                                                        >  1:866180/1‑238 (MQ=255)
                                                                                         gTTCTGATAGAGCGGGCGACGGGCCGCGATAATGGTTCTCTTCTTTTTTTCTATTTGCGGCAAGACCGTGTTTAACGCCCCCTGGGTTTCCCACCCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGc                                                                                                                                               <  2:292781/239‑1 (MQ=255)
                                                                                          ttCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTTCCCCCCCCGTGCTTTAACCGTAACCAGGCGGGCAGTAACCCCCGCCCCATTCTTTTCCCCGATATCGGGATGCCAGCCCGCCGTTTAGc                                                                                                                                               >  1:261383/1‑238 (MQ=255)
                                                                                               atatggcggggtgtgggtggggttAAGGGTGTTGTGTTTTTTTTCTATGGGCTTAAAGACCGTGTTAAGCGCCCCGTGGGTTGGCCTCCCCCAGTACAGGCGGGTTACGCGGCCAGCCAGCATTGGGCGGTTTTGCTCGCGTCCCCCCGCCCGTTCTTTAACCGTCCCCAGGGGGGCAGTACCCCCCGCCCTATTTTTCCCCCAGATGTAGAGATGCCAGCCCGCCGTTAAGCCAAAAg                                                                                                                                         <  2:503147/216‑1 (MQ=255)
                                                                                               atatAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCGCCCCTGCCTTTACCAGACCCAGGCGGGCCGTAACCTCCCCCCGATTCTTTCCCCCGGTATCGAGAAGCCAGCCCGCCGTTAAACCcaaaa                                                                                                                                         >  1:116749/1‑238 (MQ=255)
                                                                                                      tGGCGATGGGTCGCTTTATTGGTTTTTCTCTTTTTTTCATTTTGGTGCAGGCCCTTGTTTAACCCCGCCGGGGTTTGCCCCCCCCAGTACGGGCGGGCTACCCTGCCAGCCAGAAATTGGGGGGTTTGCTCGCGTTCCGCCCCCCTTGCTTAACCCGTACCCAGGCGGGCAGAACCCCCCGCCGGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGt                                                                                                                                   <  2:870078/238‑1 (MQ=255)
                                                                                                       ggCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCCGGGTCCGCCGCCCGTGCTTTTAACGTTCCCAGGCGGGCAGTAACCTCCCCCCGATTTTTTCCCCCGATATTGAGATGCCAGCCCGCCCGTAAACCCAAAGTACCGGTg                                                                                                                                 >  1:1119971/1‑239 (MQ=255)
                                                                                                       ggCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTTCCCTTTCCCCGCCCGCTGTTTAACCGGACCCCGGAGGGGCGTAACCTCCCCCCGCTTCTTTCCCCAGCTATCGGTGTGTCCGCCCCCCGCTCAGCCAAAAGTAACGTTg                                                                                                                                 >  2:979058/1‑239 (MQ=255)
                                                                                                       ggCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCCTGCTTTAACCCTACCCCGGCGGGCAGTAACCTCCCCCCCGTTTTTTTCCCCGCTATCCAAATGCCAACCCGCCGGTTAGCCAAAAAGACCCTTg                                                                                                                                 >  1:787541/1‑239 (MQ=255)
                                                                                                         cGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTTCCCCCCCCGCGCTTTTACCGTACCCAGCCGGGCAGTAACCCCCCCCCGATTCTTTTCCCCGATATTCAGATCACACCCCCCCCTTAAGACCAAAGAAcccctttg                                                                                                                               >  2:869279/1‑233 (MQ=255)
                                                                                                              ggTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCGCCGCCCCTGCTTTAACCGTACCCCGGCGGGCAGTAACCTCCCCCCGATTCTTTCCCCCGATATCCAGATGCCAGCCCGCCGTTTAGCCAAAAGTACCCTTGGGGTAcc                                                                                                                          >  1:531076/1‑239 (MQ=255)
                                                                                                                    taTAATGGTTTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTTTTTAACCCCGCCTGTGTTTGCCCCTCCCAGTACAGGCGGGCTCCGCTGCCAGCCAGCAATGGGCGGGTTTGCTCGCGTCCCGCCGCCCGTGCTTTACCCGTACCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATAGCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACggg                                                                                                                     <  1:314184/237‑1 (MQ=255)
                                                                                                                     aTACTGGTGTTCTGCTTTTTTTCTATTTGCGGCAAGACGGTTTTTACCGCCGCCTGGGTTTGCCCCCCCCAGTCCGGGCGGGCTACGCTGCCAGCCAGAAATTGGGGGGTTTGCCCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATACCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGc                                                                                                                    <  2:322603/238‑1 (MQ=255)
                                                                                                                              ttCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCCGGCGGGGAGGAACCTCCCCCCGCTTCTTTCCCCCGCAATCGAGATGCCCACCCCCCCTTAAACCAAAAGTACCGTTGGGGTACCCCGGGCCGGTCGTa                                                                                                            >  1:1212588/1‑237 (MQ=255)
                                                                                                                                         ttttttttgccggcaggccGTGTTTAACGCCCGCCGGGGTTCGCACCTCCCAGTACGGGTGGGCGCCGCTGCACGCAAGAATTTGGCGGTTTTCCTCGCTTCCCGCCGCCCGTGTTTTAACCGTCCCCGGGGGGGCAGTACCCCCCTCCCGATTCTTTCCCCAGATATCGAGCTCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGa                                                                                               <  2:835471/222‑1 (MQ=255)
                                                                                                                                                          gccggtgttaccgccggtgtggTTTCCCACCCCCCGGTACGGGCGGTCTGCGCGGCAAGCCAGCATGGGGCGGTTTTGCGCGCGTCCCGCCCCCCGTGTTTAAACCGTCCCCAGGCGGGCATAACCCCCCGCCCGATTCTTCCCCCAGATACCGAGAGGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCaaa                                                                                <  2:787541/217‑1 (MQ=255)
                                                                                                                                                             tgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTTGCGGGGTTTCTCGCGTTTCGCCCCCCCTGCTTTAACCCTACCCACGCGGGGCGCAACCCCCCCCCCATTCTTTTCCCCGATATCGCGCTGGCACACCCCCGTTTAGCCCAAAATACCGTTGGGGTACCACCGGGCTGTCGTACCCGCTGCCCGGACCGTTCACGGCAACGGccc                                                                         >  2:1172020/1‑238 (MQ=255)
                                                                                                                                                                                          ccagtaacgggCGGGCTACGCTCCCGCCAAGCAATGGGCGGTTTTCCTCCCTTCCCGCCCCCCGTGCTTTACCCGTCCCCAGGCGGGAAGTACCCTCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGc                                             <  2:1066356/230‑1 (MQ=255)
                                                                                                                                                                                           cccgTACAGGCCGGCTTCCCCTCCCGCCAAGCAATGGCGGGTTTTGCCCCCGTCCCCCCCCCCTTCCTTTACCCGTCCCCGGGCGGGAAGTACCCCCCCCCCGTTTTTTCCCCCAGATACCGGAGGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGc                                             <  2:418316/236‑1 (MQ=255)
                                                                                                                                                                                                   ggCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTTCGCCGCCCCGGCTTTTACCCCACCCCGGCGGGCAGTAACCCCCCCCCCATTCTTTTCCCCGATATCCAGGTGCCCGCCCCCCCTTTAACCCAAAAGACCGCTTGGGTACCCCCGGGCGGGCGTACCTGCCGGCGGGTCCTTCCATGGCAAAGGCCCCCCTAAGCCCTCTGCCGGCCttttttttcggcccc                                      >  2:530239/1‑228 (MQ=255)
                                                                                                                                                                                                    gCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGGGCTTTAAACCTACCCAGGCGGGCAGTAACCCCCCCCCGCTTTTTTCCCCCGATATCGAGATGCCAGCCCGCCGCTAAGCCAAAAAGACCGTTGGGGTACCCCCGGGCGGGCGTACCCCCGGCCGGGTCGTTTCGAGCAAAAGGCCCCCTTTCGCCCTCTTCCGCCATTTTTTGCCGcccccaat                                   >  2:149400/1‑235 (MQ=255)
                                                                                                                                                                                                    gCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGTTTCGCCCCCCGGGCTTTAACCGTACCCCGGCGGGCAGTAACCTCCGCCCGGTTCTTTTCCCCGATATCCGGATGCCCGCCCCCCGTTTAGCCCAAAAGAACGTTGGGGTACCCCCGGCCGGGCGTAACTGCGGCCGGGTCGGTCAGGGGAAACGGCCCCCATTAACCTTTCGCCGACCTTTTTTGCCGccccctt                                    >  2:461439/1‑235 (MQ=255)
                                                                                                                                                                                                      ggcctccGCTGCCAGCCAGCAATTGGCGGGTTGGCTCGCGTTCCGCCGCCCGTGTTTTACCCGTACCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATAGCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATc                                  <  1:1034106/232‑1 (MQ=255)
                                                                                                                                                                                                       ggCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTTCGCCGCCCCTGCTTTTACCGCACCCCGGCGGGGCGTAAACCCCCCCCCGTTCTTTTCCCCGAAATCCAGATGCCCGGCCCCCCGTTAGCCCAAAATACCGGTTGGGTACCCCCGGGCGGGCGGACCCGCCGGCGGGACGGTCCGAGCAAACGGGCCCCTTAAACCCTCTGCCGCCCTTTTTTTCCGCCCCctattc                                 >  2:689866/1‑236 (MQ=255)
                                                                                                                                                                                                       ggCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGGTCCGCCGCCCGTGCTTTAACCGCACCCAGGCGGGCAGTAACCTCCCCCCCATTCTTTCCCCCGATATCGAGATGGCCGCCCCCCGTTAAGCCAAAAGTACCGCTGGGGGACCACGGGCCGGGCGGACCCGCGGCCGGCGCGTTCAGAGCAAACGGCCCCCTTAAGCCCTCTTCCGCCCTTTTTTTCCGCCCCCTATcc                                 >  2:782493/1‑238 (MQ=255)
                                                                                                                                                                                                                 ccagccagCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCCCCGCAAAggc                      >  1:37076/1‑237 (MQ=255)
                                                                                                                                                                                                                           attgggggggTTTGCTCTCTTCCCGCCGCCCGTGCTTCAACCGTCCCCGGGCGGGAAATACCCCCCCCCGGATTCTTCCCCCAGATATGGAGTTCCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCgg             <  1:488303/232‑1 (MQ=255)
                                                                                                                                                                                                                               ggggggttttGCTCGCGTCCCCCCCCCCTTTCTTTACCCGTCCCCAGGGGGGCAGTACCCCCCCCCCGTTTCTTCCCCCAGATATCGAGATGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGa            <  2:643030/230‑1 (MQ=255)
                                                                                                                                                                                                                                       ttGCTCGCGTTCCCCCGCCCGTGCTTACACCGTACCCAGGGGGGCAGTACCCCCCCCCCGATTCTTCCCCCAGAAACCGGGATGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAg  <  2:531076/239‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
TTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAG  >  NC_000913/2227994‑2228459

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