Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,228,254 A→C 22.4% pseudogene (354/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,2540AC22.4% 87.0 / 4.7 67pseudogene (354/1287 nt)mdtQpseudogene, putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base A (36/12);  new base C (5/13);  total (42/25)
Fisher's exact test for biased strand distribution p-value = 1.22e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTA  >  NC_000913/2228022‑2228489
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAgg                                                                                                                                                                                                                                                   <  1:1304874/239‑1 (MQ=255)
gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAgg                                                                                                                                                                                                                                                   <  1:187512/239‑1 (MQ=255)
  tgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGc                                                                                                                                                                                                                                                  <  2:43019/238‑1 (MQ=255)
     cGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATTTAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCGCTTTTTCTATTTGCTGCAAGACCGGGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCAGGCggcg                                                                                                                                                                                                                                             <  2:242440/239‑3 (MQ=255)
                 gCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCCCCCCCCGGGGCTTAACCCGTCCCCAGCCGGGAGTAACccccc                                                                                                                                                                                                                                 >  1:248805/1‑237 (MQ=255)
                      aaTATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGGTCCCCCCCCCCGGCCTTAACCCTACCCCGGGGGGGAGTAAACTCCTCCCCAt                                                                                                                                                                                                                          >  1:544568/1‑239 (MQ=255)
                           cGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCCCCCCCCCCGCCTTTACCCTCCCCCCGCggggcgttaaccccccccccttcttt                                                                                                                                                                                                                       >  2:638993/1‑216 (MQ=255)
                           cGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGCTTTCTCGCGTTCCCCCGCCCGCGCTTTTACCCTTCCCCGGCGGGCCGCAAACTCCCCCCGGTTCttt                                                                                                                                                                                                                    >  1:105349/1‑239 (MQ=255)
                           cGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTAAGCTGCCAGCCAGCAATTGGCGGGCTTTCTCGCGGTCCCCCGCCCCTGCTTTTAACGCACACCGGGTGGCAGTAAACCCCCCCCCACTCttt                                                                                                                                                                                                                    >  2:779984/1‑239 (MQ=255)
                             gTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGTTCGCGCTCCGCCGCCCCTGCTTTAAACGTACCCAAGCGGGCCGTAACCTCCGCCCCAGTCtttt                                                                                                                                                                                                                   >  1:393417/1‑237 (MQ=255)
                             gTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCt                                                                                                                                                                                                                                    <  1:1097723/223‑1 (MQ=255)
                               ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTcccc                                                                                                                                                                                                                <  1:1067846/239‑1 (MQ=255)
                               ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTccc                                                                                                                                                                                                                 <  1:1269034/238‑1 (MQ=255)
                               ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTccc                                                                                                                                                                                                                 <  2:316400/238‑1 (MQ=255)
                               tgccACACCTTCAACTGAAGAAGTAATCCCGGTCTGATAGAGCTGGCGATCGGTCGCGATAATGGTGGTCTGCTCTTTTGCTATTTGCTTCAAGACCGTGTTTACCGACGCCTGGGTTTGCCACACCCAGTTCAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTGcccc                                                                                                                                                                                                                <  2:218677/237‑1 (MQ=255)
                                  cacaCCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGTCGGGGTTGCTCGCGGTCCCCCGCCCCGGCTTTAACCTCACCCACGGGGGGAGCTACATCCGCCCGGTTCTTTCCCCa                                                                                                                                                                                                               >  2:161552/1‑237 (MQ=255)
                                    cacCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCCCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGGAAACTCCGCCCGATTCTTTCCCCAGa                                                                                                                                                                                                             >  1:1117300/1‑237 (MQ=255)
                                       cTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTTCGCCGCCCGCGCCTTTAACCCACCCAGGGGGGGCGTAAACCCCCCCCgcgttctgttctttataca                                                                                                                                                                                                           >  2:1097723/1‑220 (MQ=255)
                                          cAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACTCTGCCAGCCAGCCATTGGCGGGGTTTGTCCCGGTCCGGCGCCCCGTCCTTTAACCGGCCCCGGCGGGGGACAAACCCCCGCCGCTTTCTTTccccctcactgc                                                                                                                                                                                                        >  2:751771/1‑229 (MQ=255)
                                          cAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCCGGCGGGCAGTAACCCCCGCCCCGTTTTTTCCCCCGATATCgaga                                                                                                                                                                                                    >  1:1261303/1‑239 (MQ=255)
                                          cAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCCGCGGTTCCGCCCCCCCTGCTTTTACCGCACCCCCGGGGGGCACGACCCCCCCCCCGATCTTTCCCCCAGaattgtga                                                                                                                                                                                                      >  2:588789/1‑232 (MQ=255)
                                           aaCTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAGTGGAGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGATCGGCGTTCCCCGCCCCCTCTTTTACCATACCCCGGCGGGGCGTTACCCCCCGCCCCATCTTTTccccgcaatcgcgat                                                                                                                                                                                                     >  2:510482/1‑228 (MQ=255)
                                                 agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCCCCCCCCGTGCTTTAACCGTACCCAGGCGGGCCGTAACCTCCGCCCGATTCTTTCCCCCGATATCGAGATGCCAg                                                                                                                                                                                             >  1:556770/1‑238 (MQ=255)
                                                     gTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCGCCCCCCGTGCTTTAACCGCACCCATGCGGGCCGTAACCTCCGCCCCCTTCTTTTCCCCGAGATCGTGATGCCAGCccg                                                                                                                                                                                         >  1:1159706/1‑238 (MQ=255)
                                                     gTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCCGTTCGCCGCCCCGGCTTTAACCCCACCCCGGCGGGGAGTAACCTCCCCCCCCGTCTTTTCCCACATATCCAGATGCCAGcccc                                                                                                                                                                                         >  2:1233532/1‑237 (MQ=255)
                                                             gTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCCCCGCCCGGGCTTTTACCGTACCCAGGCGGGCAGTAACCTCCCCCCGATTTTTTCCCCAGATATCGAGATGCCCGCCCGCCGTTTAg                                                                                                                                                                                >  1:615416/1‑238 (MQ=255)
                                                              ttCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTGTCGCGCTTCGCCCCCCGTGCTTTTACCCGTCCCCAGCGGGGCAGGCACTCCGCCCGGCTCTTTTCCCCATATTCGAGATGGTCGCCCGCCGTTatac                                                                                                                                                                               >  2:628952/1‑235 (MQ=255)
                                                              ttCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGGTCCGCCGCCCGGGCTTTTACCCTACCCCGGGGGGGCATAAACTCCCCCCCGTTCTTTTCCCCCATATCGGGGTTGCAGCCCCCCGTTACGc                                                                                                                                                                               >  2:1309124/1‑238 (MQ=255)
                                                              ttCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCCCGTTCCGCCCCCCGGGCTTTAACCCTACCCACGGGGGCCGGAACCCTCCCCCGGTTCTTTTCCCCCATATTGTGAATGCAAGCCCCCCGTTAAGc                                                                                                                                                                               >  1:4625/1‑239 (MQ=255)
                                                                       aGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCCCCGCCCGTGCTTTTACCGTACCCAGGCGGGCAGTAACCTCCCCCCGATTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCCAAAATAc                                                                                                                                                                      >  1:86242/1‑238 (MQ=255)
                                                                             cGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCCTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTTCCCCGATATCGAGATGCCAGCCCCCCGTTAAGCCAAAAGTACCGTTGgt                                                                                                                                                              >  1:1286388/1‑239 (MQ=255)
                                                                              gATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGgtg                                                                                                                                                             <  1:241165/239‑1 (MQ=255)
                                                                                        gcaagtGGTGTTCGCCGCTTTTTCTATTGCCCGAACGCCCGTGTTAGACCCCGCCTGGGTTGCCCCCCCCCAGTACGGGCGGGCTGCCCGGCCAGCCAGCATTTGGCGGGTTTGCGCGCGTTCCGCCCCCCGTGCTTTACCCGTACCCAGGCGGGCGGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACggg                                                                                                                                                  <  2:1002013/233‑1 (MQ=255)
                                                                                        gATAATGTTGTCCTCCTTTTTTTCTTTTTGCTGCAAGACCGTTTTTAACGCCGCCGGGGTTGGCCCCCCCCAGTACGGGCGGGCTACGCTGCCAGCCAGCAATGGGGGGGTTTGCTCGCTTTCCGCCGCCGGTGCTTTAACCGTCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACggg                                                                                                                                                  <  2:1261303/238‑1 (MQ=255)
                                                                                              ggTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACCCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGTTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCg                                                                                                                                          <  1:1130764/239‑1 (MQ=255)
                                                                                                  ttCTGCTTTTTTTCTATTGTCGGCAAGCCCGTGTTTACCGCCGCCGGGGTTCGCCACCCCCAGTACAGGCGGGCTACGCTCCCGCCCAGCAATGGGGGGGTTTGCTCGCGTCCCGCCGCCCGTTCTTTACCCGTCCCCAGGGGGGCAGTACCCTCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTAcc                                                                                                                                      <  2:615416/239‑1 (MQ=255)
                                                                                                    cTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCCGTAACCTCCGCCCGATTCTTTTCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCt                                                                                                                                     >  1:860516/1‑238 (MQ=255)
                                                                                                            tttgcgtttGCTGAAAGCCGTTTTTACCCGCCGCCGGGGTTGGCCACGCCCACACACGGGGCGTTACGCTGCCAGCAGGCAGTGGTGGGTTTGGCTCGTGTCCCGCGGCCGATGCTTGACCAGGACCATGGCGGGAGTGACCCTCCGCCCACTCTTTTCCCCGAATACCGTGAGTCCTGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGa                                                                                                                            <  2:763338/233‑1 (MQ=255)
                                                                                                            tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGCGCTTTTACCGTACCCAGGCGGGGCGTAACCCCCGCCCGCTTCTTTCCCCCGATATCCAGATGCCCACCCGCCGGTAAGCCCAAAGTAACCCTGGTGTACCACCGGGCGGCCGTACCCTCGGccacg                                                                                                                            >  2:335867/1‑236 (MQ=255)
                                                                                                                 aTTGGCTGAAAGACCGTTTTACACGCCGCCTGGGTGCGCCACCCCCAGTACGGGCGGGCTACCCTGCCCGCCAGCAATGGGCGGGTTTGCTCGCGTCCCGCCCCCCGTGCTTTAACCGTCCCCGGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGtt                                                                                                                       <  2:849228/239‑1 (MQ=255)
                                                                                                                    tgctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTTCCAGCCCGCAATTGGCGGGGTTTCTCGCCTTCCGCCCCCCCGGCTTTAACCCTACCCCGGCGGGCCGGAACCTCCCCCCCCTTCTTTTCCCCGATATCTCGATGCCCGCCCCCCCCTAAACCCAAAGTAACGCTTGGGTGCCCCCGGCCCGCCCTACCCGCGCCCCGCTCGTTCac                                                                                                                    >  2:1185461/1‑238 (MQ=255)
                                                                                                                    tgctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGGTCCGCCGCCCGTGCTTTTAACGTACCCAGGCGGGCCGGAACCCCCCCCCCCTTCTTTTCCCCCAAATCCAGATGCCCGCCCCCCGGTTAACCAAAAGTACCCCTTGTGGACCCCCGGCCGGGGGGACCTCCgggccgggcgtttag                                                                                                                    >  2:1207028/1‑226 (MQ=255)
                                                                                                                            gACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTGTCCCGTTCCCCCGCCCCTGCTTTAACCCTAACCCGGCGGGGCGTAACCTCCCCCCGCTTCTTTTCCCAGATATCCAGATACCCGCCCCCCGTTAAGCCCAAAAATCCCTTGGTGGACCCCCGGCCGGTCGGAACTGCGGCCGGAGCGGTCAGAGCAAACg                                                                                                           >  2:372653/1‑239 (MQ=255)
                                                                                                                              ccGTGTTTAACGCCGCCGGGTTTTGCCACCCCCAGTACGGGCGGGCTCCGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGc                                                                                                         <  1:682587/239‑1 (MQ=255)
                                                                                                                                 tgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCCCCCCCCCCGGTTTTAACCTACCCCAGCGGGGAGTAACCTCCCCCCCCTTCTTTTCCCCGATAACCAGATTCCCGGCCCCCCCTTAGCCCAAAATACCCGTGGGGGACCCAGGGCCCGGCCTACccgccgccgcggtgttcctgggaaacggccc                                                                                                        >  2:6426/1‑216 (MQ=255)
                                                                                                                                 tgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGGCCCCGGTCCCCCGCCCCGGCCTTATACGCACCCCGGCGGGCCGGAACCCCCCCCCCGATCTTTCTCCAGATATCCCGATGGCCCCCCCCCTCTAAACCAAAAGTACCCTTGGTGTGCCACAGGCCGGGCGTCACCGTCGCCGGGTAGTTCCAAGCCAAACGGcc                                                                                                        >  2:180379/1‑238 (MQ=255)
                                                                                                                                         cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGGTCCCCCCCCCCTGGTTTAACCCTACCCCGGCCGGCAGGAACCTCCCCCCCCTTCTTTTCCCCCATATCCTGATGGCCACCCCCCCGTTAACCCAAAATACCGTTTGGTTACCCACGGCCGGTCGGTCCCCGCGCCGGGTCGTTCACGGCACACAGcccccctaaa                                                                                                >  2:671070/1‑232 (MQ=255)
                                                                                                                                              ccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTGTCCCCGTCCGCCCCCCCGGGCTTAACCCGAACCCCGGCGGGAGTGACCTCCCCCCGCTTTCTTTCCCCCAATACCAGGTGCCCGCCCCCCCTTAAAGCCCAAATACCCTTGGTTGTACACCGGCCCGTCGTACCCTCGGCCTGGCTGCTCAGAACAAACACCCCCAATAAACCTTc                                                                                           >  2:1303864/1‑238 (MQ=255)
                                                                                                                                              ccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCCCGCGGTCCGCCGCCCGGGCTTTTACCCTACCCAGGCGGGGAGTAAACCCCGCCCGGTTCTTTCCCCCGATATCGAGATGCCAGCCCCCCCGTAAGCCAAAAATAACCTTTGGGTACCACGGGGCGGGCGTACCCCCGGCCGGAGCGTTCAGAGGCAACGGCCCCCCTTAACCTTCt                                                                                          >  2:505640/1‑239 (MQ=255)
                                                                                                                                              ccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCCCCGCCCCCGCTTTTACCCGACCCACGCGGGGCGTAAACTCCCCCCCGTTCTTTTCCCCGATATCGAGATGCCCACCCCCCGCTAAGCCCAAAATAACGTTGGTGTAACACCGGGCCGGCGTACCCGCGGCCCGGACGGTCAGAGCAAACGGGCCCCTTTAACCTTCt                                                                                          >  2:496678/1‑239 (MQ=255)
                                                                                                                                                   ttttGCCCCCCCCATAACAGGCGGTCTACGCTGCCAGCCAGCAATTGGGGGGTTTCCTGGCGTCCCGCCCCCCGTGCTTTAACCGTCCCCGGGGGGGAAGTACCCTCCCCCCGATTCTTCCCCCAGATATCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCg                                                                                      <  2:544568/237‑1 (MQ=255)
                                                                                                                                                    gtttcccccccccgtgcggggCGGGCTCCGCGCCCCGCCAGCATTTGGCGGGTTTGCCCGCTTCCCCCCGCCCGTTTTTTACCCGTCCCCAGGCGGGCGTTACCCCCCGCCCGATTTTTCCCCCAGATAGCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCg                                                                                      <  2:1315979/219‑1 (MQ=255)
                                                                                                                                                        ccACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGTTCCGCCGCCCGTGCTTTAACCCTACCCAGGCGGGCAGTAAACTCCCCCCGATTCTTTCCCCCGATATCGAGATGCCCACCCGCCCTTAAACCCAAAGTACCCTTTGGGTACCCCCGGGCCGGCGGACCTGCGGCCGGGTCGTTCCGAGCAACAGGCCCCACTAAGCCTTCTGCCGCCAttt                                                                                >  2:230716/1‑239 (MQ=255)
                                                                                                                                                             cccAGTACAGGCGGTCTCCGCTGCAGGCCAGAAATTGGCGGGTTTCCTCGTGTTCCCCCGCCCGGTTTTAAACCGTCCCCAGGGGGGCAGTACCCTCCCCCCGATTTTTCCCCCAGATAGCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGc                                                                           <  2:1267132/239‑1 (MQ=255)
                                                                                                                                                                        gCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTGTCGCGGTTCGCCGCCCGTGCGTTTACCGTCCCCAGGCCGGGAGTAACCTCCCCCCGCTTTTTTTCCCACGGATCGAGATGCCCGCCCCCCGTTAAACCCAAAATAACGTTGGGGGTCCACGGGCCGGGCGTGCGCGCGGCCCGGGCGGTCACGGCCAACGGGCCCCCTAACAACTCTCTCGACCTTTTTTGCTGCCCCatt                                                                  >  2:168018/1‑237 (MQ=255)
                                                                                                                                                                        gCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTTCGCCGCCCCTGCTTTTACCGTACCCCCGCGGGGAGTAACCTCCGCCCCATTCTTTTCCCCGGTATCCCGAAGCCAACCCGCCGTTAAGCCAAAAGGACCGCTGGTGTACCACGGGCCGGGCGTGCCTGCGGCCCGGTCCGGCTGAGCCAACTGCCCCCCTAAGCCTCCTCCCGACCTTTTTTTCCGCTccctc                                                                  >  2:1023423/1‑235 (MQ=255)
                                                                                                                                                                             cTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTTCGCCCCCCGTGCTTTTACCCCACCCAGGCGGGGAGTAACCCCCGCCCCCTTTTTTTCCCCAATATCGCGAAGCCACCCCGCCCCTAAGCCCAAAAGACCGTTGGTGTACCACGGGCCGGGCGCACCCGCGGCCGGGTCGTTTCGCGCAAAAGGGCCCCCTAAACCCTCTGCCGGCAtttttttgcccgccacatcccc                                                             >  2:1224294/1‑222 (MQ=255)
                                                                                                                                                                                        gccagcATTGGGGGGGTTGCCTCGCTTTCCGCCGCCCTTGCTTTACCCGTCCCCAGGCGGGCATTACCCCCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGAt                                                 <  1:1233532/239‑1 (MQ=255)
                                                                                                                                                                                                   tgggggTTTGCCCCGCGTTCCGCCCCCCGTGCTTAAACCGTCCCCAGGGGGGAAGAACCCCCCGCCCAATTTTTCCCCCAAAAATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGa                                          <  1:163008/234‑1 (MQ=255)
                                                                                                                                                                                                                      gTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTACGCCTTCTGCCGACATTTTTTGCCGCCCCATCTCCGCCGCAAAGTCCACCTGCGGACCATCCTGACTGGcg                           >  2:1094475/1‑235 (MQ=255)
                                                                                                                                                                                                                         ccgcccCCCGTGCTTAAACCGTACCCAGGCGGGAAGTAACCCCCCCCCGATTCTTCCCCCAGATACCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCAGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGcc                     <  1:335867/239‑1 (MQ=255)
                                                                                                                                                                                                                                        cTTTAACCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCAAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCt         <  1:505640/239‑1 (MQ=255)
                                                                                                                                                                                                                                         ttttAACCGTCCCCGGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCt         <  1:230716/238‑1 (MQ=255)
                                                                                                                                                                                                                                                 gTACCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCACAACTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTACAGAGCAACCGGCCCCATTAAGCCTTCCGCCGCCCTTTTTTGCCGCCCCATATCCGCCGAAAAGTCGATCTGCGGACCACCCTGAGCGGCACCTGCCTTCGCCGCGCCTTCaactc  >  2:166032/1‑237 (MQ=255)
                                                                                                                                                                                                                                                 gTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta  >  2:934381/1‑238 (MQ=255)
                                                                                                                                                                                                                                                  tACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta  >  2:741399/1‑237 (MQ=255)
                                                                                                                                                                                                                                                  tACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta  >  1:715508/1‑237 (MQ=255)
                                                                                                                                                                                                                                                   aCCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta  >  2:77098/1‑236 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
GCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTA  >  NC_000913/2228022‑2228489

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