Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,228,229 G→C 17.0% 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.0% 109.2 / 3.3 51pseudogene (381/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base G (16/25);  new base C (7/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 6.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

TTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCA  >  NC_000913/2227996‑2228454
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             
ttcattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGTTAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCCGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGAGGGGTTTCTcgccgtc                                                                                                                                                                                                                                          >  2:1366354/1‑235 (MQ=255)
  cattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTccg                                                                                                                                                                                                                                        <  1:1268961/239‑1 (MQ=255)
   attttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTccgc                                                                                                                                                                                                                                       <  2:490828/239‑1 (MQ=255)
       tcCCCTCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGAGATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCgccgcc                                                                                                                                                                                                                                   <  1:740733/239‑1 (MQ=255)
       tcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCgccgcc                                                                                                                                                                                                                                   <  1:1050857/239‑1 (MQ=255)
       tcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGCGCTACCCTGCCAGCCCGCAATTGGCGGGTTTTGTCTCGTTTccccgc                                                                                                                                                                                                                                    >  2:524944/1‑238 (MQ=255)
       gccccGCGACATCGGTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCAAACCTTAAAATGAGGAAGTAGTCCCGTTCTTATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCGATCTGCTGCAAGACGGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCgccgcc                                                                                                                                                                                                                                   <  2:731897/238‑1 (MQ=255)
              gACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGTGGGTTTGCTCGCGTTCCGCCGCCCGTGCtt                                                                                                                                                                                                                           <  2:651489/239‑1 (MQ=255)
                 aTCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTaa                                                                                                                                                                                                                        <  2:1040693/239‑1 (MQ=255)
                 aTCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTaa                                                                                                                                                                                                                        <  1:718385/239‑1 (MQ=255)
                      tGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTa                                                                                                                                                                                                                   <  2:1272767/239‑1 (MQ=255)
                          ctgctgCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTAccc                                                                                                                                                                                                                <  2:699243/238‑1 (MQ=255)
                            gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGCCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAgg                                                                                                                                                                                                             <  2:275186/239‑1 (MQ=255)
                            gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAgg                                                                                                                                                                                                             <  1:744827/239‑1 (MQ=255)
                             ctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTGTTCGCTTTCCCCCCCCCCGCCTTTACCCGTcccccgc                                                                                                                                                                                                             >  2:849964/1‑236 (MQ=255)
                             ctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCCTAccccagc                                                                                                                                                                                                            >  1:658867/1‑235 (MQ=255)
                              tgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGc                                                                                                                                                                                                            <  1:1394819/238‑1 (MQ=255)
                              tgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTCTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGc                                                                                                                                                                                                            <  1:754416/238‑1 (MQ=255)
                               gcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTTCGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCgg                                                                                                                                                                                                          <  1:459330/239‑1 (MQ=255)
                                 cGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCggg                                                                                                                                                                                                         <  1:812763/238‑1 (MQ=255)
                                                         gTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCCCCCGTGCTTTAACCGTACCCAAGCGGGCAGTAAACTCCGCCCGGTTCtttt                                                                                                                                                                                >  1:33122/1‑237 (MQ=255)
                                                           ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTccc                                                                                                                                                                              <  1:293422/238‑1 (MQ=255)
                                                            tCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCGCCGCCCGTGCTTTTACCGTACCCAGGCGGGCAGTAACCCCCCCCCGGTTCTTTTCCCa                                                                                                                                                                            >  1:454244/1‑239 (MQ=255)
                                                                    ttCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCCCCCGGGCTTTAACCGTACCCCGGCGGGGCGTAACCTCCCCCCGATTCTTTcccccagttttca                                                                                                                                                                    >  1:1065543/1‑230 (MQ=255)
                                                                             agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGGAACCTCCGCCCGATTCTTTCCCCCGATATCGAGATGCCAg                                                                                                                                                         >  1:1159355/1‑238 (MQ=255)
                                                                             agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCCCCCCCCGTGCTTTTACCCGACCCCGGCGGGCCGTAACCCCCCCCCCGTTCTTTTCCCCCGTATCCCGATGcccg                                                                                                                                                         >  2:300563/1‑236 (MQ=255)
                                                                             agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTGTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCCGTAACCTCCGCCCCATTCTTTCCCCCGATATCgcgcagccaa                                                                                                                                                           >  1:1060747/1‑229 (MQ=255)
                                                                                          ttCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCCTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGGAAACTCCGCCCGATTCTTTTCCCAGATATCGCGATGCCAGCCCCCCGCTAAGcc                                                                                                                                         >  1:892628/1‑239 (MQ=255)
                                                                                               atatAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTTCCCCGCCCGTGCCTTTACCCTACCCCAGCGGGCCATAACCTCCCCCCCGTTCTTTCCCCCAAGATCCCGAAGCCCACCCCCCGTTTAGccccaaa                                                                                                                                    >  2:1247715/1‑234 (MQ=255)
                                                                                                         cGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTTTCGCCGTCCGCCACCCCGGCTTTAACCCGACCCCGGCGGGGAGTGAACTCCGCCCGGTTCTTTTCTCAGATAATGAGGCGCGACCCCGCCCCTTAGCCCCAAGAACCGCTggg                                                                                                                          >  2:76914/1‑238 (MQ=255)
                                                                                                         cGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCCGGCTTTTACCCTAACCAGGCGGGCCGTAACCTCCGCCCGGTTCTTTTCCCCGATATCGAGATGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTTgt                                                                                                                          >  1:509571/1‑239 (MQ=255)
                                                                                                               ggcgcaatatgggTTTTCTGCTTTTTTTCTTTTTGCTGAAGGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACGGGGGGGCTACGCTGCCACCCAGCAATGGGCGGTTTTGCTCGCGTCCCGCCGCCCGTGCTTTACCCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTAcc                                                                                                                     <  2:1159355/227‑1 (MQ=255)
                                                                                                                              ttCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCCCCGTTTGCCGCCCCCGCCTTTACCCTACCCCCGCGGGCCGTAAACCCCCCCCCCCTCTTTTCCCAGCTATCCACATGCCCAGCCCCCCGTAACCCCAAAAAACCCTTTGGTTGCCCCGGGCCGGGCGtca                                                                                                     >  2:10902/1‑236 (MQ=255)
                                                                                                                                        tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGAGGTCCGCCCCCCCGGCTTTTAACCTTACCCAGCGGGGGATTAACCCCGCCCCCTTCTTTCCCCCAATATCGAGATGCCACCCCCCCGTTAAGCCCAAAATACCCTTGGTGTACCACGGGCCCGGCGGACCTGAGGcccgg                                                                                         >  2:704923/1‑236 (MQ=255)
                                                                                                                                        tttCGATTTGCGGCAAGACCGTGTTTACCGCCGCTGGGGTTTGCCCCTCCCAGTACGGGCGGGCTTCGCGGCCAGCCGGAAATGGGCGGGTTTGCTCGCTTCCCGCCGCCCGTGCTTTAACCGTACCCGGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGa                                                                                        <  2:1065543/239‑1 (MQ=255)
                                                                                                                                                          ccGTTTTTACCCCCCCCTGGGTTTCCACCCCCCTGTGCGGGCGGCCTCCCCCCCCAGCCAGCAATTGGCGGGTTTGCCCCCTTTCCGCCCCCCGTCCTTAACCCTTCCCCGGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATACCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGc                                                                     <  2:33122/239‑1 (MQ=255)
                                                                                                                                                             tgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCCCCCCCTTCCCCGCCCCCGGTTTTAACCGACCCCGGCGGGGCGTAACCCCCCGCCCGTTTTTTTCCCCAGTAACCAGAATGCCACCCCCCCGTTAACCCAAAAGACCCGTGGGTGTACCCGGGCCGGTCGGCACCGCGCCCCGAACGTTCCGAGCACACGGccc                                                                   >  2:679777/1‑238 (MQ=255)
                                                                                                                                                             tgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCCACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCCCCGCCCGTGGTTTAAACGTACCCAAGCGGGCAGGTAACTCCCCCCGATTCTTTTCCCCGATATCGAGGAGCCAGCCCGCCCGTAAGCCCAAAATGCCGCTGGGGTACCAACGGCCGGTCGTACCTGCGGGCGGGTCGTTCCGGGCGAAACGGcc                                                                    >  2:416611/1‑238 (MQ=255)
                                                                                                                                                                         ccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCCGTCCTCCCCCCCGGCTTTAACCCCACCCCGGCGGGCCGTAACATCCGCCCCCTTCTTTCCCCACGTATCCAGATGCCAGCCCGCCCGTAAACCCAAAATACCCCTTGTGGACCACGGGGCCGGCGGACCCGCGGCCGGGATCGTCCAAGCAACCGCCCCCCTTACGCCTTc                                                       >  2:410459/1‑238 (MQ=255)
                                                                                                                                                                                       tCCCATTACAGGCGGGCTCCGCTGCCAGCAAGCAATGGGCGGGTTTGCCGGCGTTCCGCCGCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTAACCCCCGCCCTATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTg                                       <  1:1015008/239‑1 (MQ=255)
                                                                                                                                                                                        cccGGTACAGGCGGGCTCCGCTGCCAGCCAGCAATTGGCGGGTTGGCTCGCGTCCCGCCGCCCGTGCTTTAACCGTCCCCAGGGGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGc                                      <  2:509571/239‑1 (MQ=255)
                                                                                                                                                                                           aGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCg                                                                                                                                                                                        <  2:212451/100‑1 (MQ=255)
                                                                                                                                                                                              aCAGGCGGGCTACGCTGCCAGCCACCAATTGGCGGGTTTCCTCGCGTCCCGCGGCCGGTGCTTTAACCGTCCCCGGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATACCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGc                                   <  1:1247715/236‑1 (MQ=255)
                                                                                                                                                                                                  gCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCCCCCGTTCCCCCCCCCGCGCTTTTAACCTACCCCGGCGGGCAGGAACCTCCCCCCCATTTTTTTCCCCGATATCGGGAAGGCCCCCCGCCGTTTACCCCAAAAAACCCGTGGTGCACCACCGGGCGGTCGGTACCCTCGTCCGGGAGCTACGGCACAAACGGCCCCCTTAAGCCCTCTGTCGCAAttttttacacgcccttt                             >  2:1000316/1‑228 (MQ=255)
                                                                                                                                                                                                  gCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTGTCCCGGTCCGCCGCCCGTGCTTTAACCGCAACCAGGCGGGGCGGAACCTCCGCCCCGTTCTTTCCCCCGATATCGAGATGCCACCCCCCCGTTTAGGCCAAAATACCCTTGGGGTACCACGGGCCGGGCGTGCCCGGGGCCGCGGCGTTCAGAGCAAAACGCCCCACTAAGCCCTCCTCCGGCCTTCTTTTCCggccccac                             >  2:831934/1‑231 (MQ=255)
                                                                                                                                                                                                  gCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCCCCGCCCCTGCCTTAACCGTACCCACGCGGGCAGTAAACCCCGCCCCATTCTTTCCCCCGCTATCGAGAGGCCAGCCCCCCGGTAAGACCAAAAGAACGTTGGTGTACCACGGGGCCGCCGGAACCGCGGCCCGATCGTTCCTAGCAAACCGCCCCATTAACCCTCCTGCCGCCCATTTTTTTccctcctac                              >  2:22901/1‑231 (MQ=255)
                                                                                                                                                                                                     ggCTACGCTGCCAGCCAGCAATTGGCGGGTTTTGTCCCGTTTCGCCCCCCGCGCCTTAACCGTACCCCGGCGGGGAGTAACCTCCCCCCCATTCTTTccc                                                                                                                                                                              >  1:212451/1‑100 (MQ=255)
                                                                                                                                                                                                     ggCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGTTCCCCCGCCCCGGCTTTTAACCTACCCCGGCGGGGAGGAACCTCCCCCCCACTTTTTTCCCCCATAACGAGAAGCCCGCCCCCCGGTTAACCCCAAACAACCGTTGTGTACCACCGGCCCGTCGGACCTGCGGGCCGGACGGTCCAGGCCAAACGCCCCCATTACCCTTTCTCCGGCCATTTTTTGCCGcccaaactc                           >  2:334400/1‑230 (MQ=255)
                                                                                                                                                                                                                    cagcAATTGGCGGGTTTGCGCGCGTCCCGCCCCCCGTGTTTTAACCGTACCCGGGCGGGCATTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCt          <  1:853612/239‑1 (MQ=255)
                                                                                                                                                                                                                         atgggCGGGTTTGCTCGCGTCCCGCCGCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGa     <  1:300563/236‑1 (MQ=255)
                                                                                                                                                                                                                         aggggCGTGTTTGCCCCCTTCCCGCCCCCCGTACTTTACCCGTCCCCAGGCGGGCACTAACCTCCCCCCTATTCTTCCCCAGTATACCGAGCGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCGTATCCGCCGAAAAGTCGATCTGCGGa     <  2:645117/236‑1 (MQ=255)
                                                                                                                                                                                                                             gCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGCTCTGCGGACCa  >  1:1086643/1‑238 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             
TTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCA  >  NC_000913/2227996‑2228454

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: