Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,228,229 G→C 17.6% 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,2290GC17.6% 116.6 / 3.4 55pseudogene (381/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base G (18/25);  new base C (8/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 3.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

TTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAG  >  NC_000913/2227996‑2228461
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
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)
                                                                            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)
                                                                                               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)
                                                                                                                                                             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/2227996‑2228461

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