Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 134,775 T→C 19.4% intergenic (+25/+13) yacL → / ← speD UPF0231 family protein/S‑adenosylmethionine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913134,7750TC19.4% 102.8 / 5.1 73intergenic (+25/+13)yacL/speDUPF0231 family protein/S‑adenosylmethionine decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (36/19);  new base C (2/14);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 3.59e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTCACATCCACCATGTCATACAGCGCCTTCA  >  NC_000913/134550‑134996
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              
tCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTTGCGGCTTACCGCAATTTCGTGCAGCAGAAGGAAATAAATTTGGCGGAGC‑CGGGGGGCTCCg                                                                                                                                                                                                                             >  1:720133/1‑229 (MQ=255)
tCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGGGCTCCg                                                                                                                                                                                                                             >  1:153145/1‑229 (MQ=255)
tCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCATGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCCATTTGTTGCACAGGAAGTAAAAAAATTTTGCGGGGC‑CTGGGGGCTCCg                                                                                                                                                                                                                             >  2:108222/1‑229 (MQ=255)
     gcacggggcaggggaGGAATCCCCCCTGTGGATGGTGGGGGAAGAGGTGATGGTTCGGTCAAATAAACTGAAATTCGCGGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGGGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCag                                                                                                                                                                                                                         <  2:79894/223‑1 (MQ=255)
          gggcagggcaTGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CGGGGGGCTCCGCCAGAGCCGt                                                                                                                                                                                                                   >  1:446867/1‑229 (MQ=255)
          gggcagggcaTGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCCCATTTTGGTGCAGCAGAGGGAAAAAAACTTGGGGGGGC‑CTGGGGGTCCCGcccagaccct                                                                                                                                                                                                                   >  2:364882/1‑221 (MQ=255)
               gggcaTGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTAAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCTGGCAGCAGAAGCAAAAAAATCGTGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTacaa                                                                                                                                                                                                              >  1:693839/1‑226 (MQ=255)
                     gAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCCCCAGAGCCGTTAAACAGCTg                                                                                                                                                                                                         >  1:367664/1‑228 (MQ=255)
                          cacCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGTCGGAGC‑CTGGGAGCTCCCCCAGAGCCGTTAAACAGCTGGCaaat                                                                                                                                                                                                   >  1:392551/1‑229 (MQ=255)
                           acCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGATCCGCCAGAGCCGTTAAACAGCTGGCCTAtt                                                                                                                                                                                                  >  1:647617/1‑229 (MQ=255)
                            cccTGTGGATGGCCGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGACGAAGTAAATAAATCGGGCGGAGC‑CTGGGAGCTCCGCCAAGGCCGTTAAAAACATGGCAAAtta                                                                                                                                                                                                 >  2:281126/1‑228 (MQ=255)
                               tgtgGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGGGCCGTTAAACAGCTGGCATATTgcgc                                                                                                                                                                                              >  1:234421/1‑229 (MQ=255)
                                gtgGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCCGAGC‑CTGGGGGCTCCGCCCGAGCCGTTAAACAGCTGGCCAATTtcgc                                                                                                                                                                                              >  1:643183/1‑228 (MQ=255)
                                 tggatggaCGGGGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGTCGGAGC‑CTGGGGGGTCCGCCAGAGCCGTTAAAAATCGGGCATATTGCGcc                                                                                                                                                                                             >  1:267104/1‑228 (MQ=255)
                                 ggggtggaCGGGGAAGGGGTGGTGTGTCGCGCCAATCAACTGAAATGCGGTGGCGATGAAATGGAAGAGGGTATGAACTAATCCGACGAAGAAAGCCTTGCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCCCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGcc                                                                                                                                                                                             <  1:376387/224‑1 (MQ=255)
                                 gggatgggCGGTGAAGGGGTGGTGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGGGGGGATTAACTACTACGACGAAGAAAGCCTGTCGCTAGGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTgcgc                                                                                                                                                                                              <  1:364882/226‑1 (MQ=255)
                                           gTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTAACGCAATTTCGTGCAGCCGAAGTAAATAAATCTGGCGGGGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCGGGCACATTGGGCCCGTAATAAAtc                                                                                                                                                                                  >  1:640269/1‑229 (MQ=255)
                                           gTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGAGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCATAAGTAAATAAATCTGGGGGAGC‑CTGGGGCCCCCGCCAGAGCCCTAAAACACATGGCAATCTGGGCCCGtaactaaca                                                                                                                                                                                  >  2:399512/1‑223 (MQ=255)
                                                   ggggtGGTTCGCGCCAATCAACTGGAATTCGGTGGCGATGAAATGGAAGAGGGGATGACCTCCTACGACGAAGAAAGCCTGTCGCTAGGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCAAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCa                                                                                                                                                                           <  1:544859/224‑1 (MQ=255)
                                                   ggggtGGGTCGCACCAATCACGGGGAAGTCGGTGGCGATGAAATGGAAGAGGGGAGGAACTACTACGACGAAAAAGGCCTGTCGGTATGGGGTGTTGAGGATTTTAT‑GA‑AG‑GTCGTGGGGGCTTACCGCAATTTCGTGCAGCAGAATTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCa                                                                                                                                                                           <  2:153145/224‑1 (MQ=255)
                                                      aTGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCCGAGCCGTTAAAAAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCAttt                                                                                                                                                                        >  1:384224/1‑228 (MQ=255)
                                                      aTGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAAACGCTGGAATATTGCGCCCGGAATAAATCTCGCGCAatt                                                                                                                                                                        >  1:150945/1‑228 (MQ=255)
                                                      aTGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGCCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTCGGCCCGTAATAAATCTCGCGCAttt                                                                                                                                                                        >  1:416039/1‑228 (MQ=255)
                                                      aTGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGGAAACAAATCTGGCGGAGC‑CCGGGCGCTCCGCCAGAGCCGTTAAAAAACTGGCCTATTTGCACCGTAATAAACCTCCCGCCtttt                                                                                                                                                                       >  2:745622/1‑228 (MQ=255)
                                                            cgcgCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAGGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGCGCCGTTAAACACCGGGCATATTCCGCCCGTACTAAATCTCGCGCATTTCTTTCca                                                                                                                                                                 >  1:380034/1‑229 (MQ=255)
                                                            cgcgCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCATAGCCGTTAAACAGCTGGCATATTGCGCCCGGAATAAATCTCGCGCATTTCTTTccc                                                                                                                                                                 >  1:532121/1‑228 (MQ=255)
                                                               gcaaagcAACGGGAATTGGCTGGCAGGGAAATGGAAGAGGGGATGAACTACGACGACGAAGAAAGCCTGTCCTTATGCGGCGTTGAGGATTTTCG‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTGCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCCCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCAcag                                                                                                                                                              <  2:522622/223‑1 (MQ=255)
                                                                                  cgggCGGTGAAAGGGAAGGGGGGATGAACTCCTACGACGAAGAAAGCCTGGCGCTATGCGGCGTTGGGGATTTTCT‑GC‑AG‑GGCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTcgcg                                                                                                                                                                             <  2:485839/193‑1 (MQ=255)
                                                                                  cTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTccctgcctcttttacacacctccg                                                                                                                                            >  1:485839/1‑206 (MQ=255)
                                                                                  cTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTGATTTCCTgg                                                                                                                                           >  1:657344/1‑229 (MQ=255)
                                                                                                 ggggggggATGAATTCCTAGGACAAAGAAAGCTTGTCGCTATGCGGCTGTGAGGATGTTCG‑GG‑AG‑GTCGGGGCGGCTTCCCGAAATTCCGTGCAGCAGAAGTAAATAAATCGGGCGGAGC‑CGGGGAGCCCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGtt                                                                                                                            <  2:450247/225‑1 (MQ=255)
                                                                                                               tcttacGAGCAAGAAAGCCTGCCCCTATGGGGCTTTGGGGATTTTCG‑GC‑GG‑GTCGTGGCGGTTTCCCGCATTTTCGGGCAGAAGAAGAAAATAAAGCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGtt                                                                                                               <  2:367664/225‑1 (MQ=255)
                                                                                                                actacGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTCATTTCCTGGCGCTCGCTGTCGGTTAGGTCTTCCGGtt                                                                                                               >  1:600107/1‑227 (MQ=255)
                                                                                                                actacGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCGCGCTGCCGGTTAAGTCCTCCGGGTTg                                                                                                             >  2:321633/1‑229 (MQ=255)
                                                                                                                 ctacGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCATGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATAATGCGCTCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCACCGGGAATATCCTGGCGCTCCCTGTCGGTTAAGGCTTCCCGTTTgg                                                                                                            >  1:683673/1‑229 (MQ=255)
                                                                                                                  tctgtcgaAGAAGGTCTGTTCCTATGCGGGGTTGAG‑‑TTTTTT‑GCAAGAGTGGTGGCGTCTCCGCGAATTTTGGTCCATAAGACGTAATTACTGCGGGCGGAGC‑TTGGGATCTCCGCATGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCAGTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTGGCTGTCGGTTAAGTATTCCGGTTTGgt                                                                                                           <  2:131971/224‑1 (MQ=255)
                                                                                                                   acgagaaagaaagctggtgcttATGTGGCGATGAGGATTTTCT‑GC‑GG‑GTGTGGGCTGCTTACGGCATGTTCGAGAAGCAGAAGTAAATAAGTCTGGCGTGCC‑CTGGGAGCACGCCCAGAGCCGTTAAACAGCTGTCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGgt                                                                                                           <  2:108958/207‑1 (MQ=255)
                                                                                                                      acgaAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCCCAGCCCAGGGGGAATTTCCTTGGGCCCGCCGTCGGTTAAGTCTTCCCGGTTggggggg                                                                                                       >  2:666221/1‑224 (MQ=255)
                                                                                                                           gaaAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGGAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACAt                                                                                                  >  1:423893/1‑229 (MQ=255)
                                                                                                                            gaggctctctgccttGCGGCGTTGAGTATTTTCG‑GC‑GG‑GTGGTGGCGGCTTACCGCAATTTCGTGCAGCCGAAGTAAAAAAATCGGGGGGCGC‑CTGGGACCTCCGCACGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATg                                                                                                 <  2:657343/215‑1 (MQ=255)
                                                                                                                              agcttttggCTAGGGGGCGTGGAGGTTTTGTG‑GC‑GG‑GTCGTGGCGGCTTCCCGAAATTTCGGGCAGCGGAAGTAAATAAACCGGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTa                                                                                               <  2:392551/221‑1 (MQ=255)
                                                                                                                                   gTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCCCAGCGCAGCGGTAAGTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGCTTGGGGTGGagacggtagggc                                                                                           >  1:603190/1‑217 (MQ=255)
                                                                                                                                   gTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGGCTTCCGGTTTGGGGTGGAACATGTAGGGc                                                                                           >  1:195196/1‑228 (MQ=255)
                                                                                                                                     cctttcgggggttttgagtttttGT‑GG‑AG‑GTTGTGGGCGCTGTCCGCCTCTTCGTGAAGCGGAAATAAAAAAGTCTGGCGGAGC‑CTGGGCGCCCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTa                                                                                        <  2:347462/210‑1 (MQ=255)
                                                                                                                                            tggcgTTTGGGATTTTGC‑GC‑AG‑GGCGTGGCGGGTTGCCGCAATTTCGTGCAGCAAGAGAAAAAAAATCCGGGGGGGC‑CTGGGACCCGCGCCAGGGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGa                                                                                 <  2:301029/228‑1 (MQ=255)
                                                                                                                                            gggctTTGAGGTTTTTTT‑GC‑AG‑GTCGGGGCGGGTTCCCGGCATTGAGTCAAGGCATAGAAAATAAATCGGGCGGGGC‑TGGGAAGATCCGCCGGAGGCGGTAAACAGCTGGCATATTACGCCCGTAATAAATCTCGCGCAGTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGa                                                                                 <  2:635774/228‑1 (MQ=255)
                                                                                                                                            cggtttgggtgtgtttCT‑GC‑AG‑GGCCGGGGAGTTTCCCGAAATTCCGGGAAGCAGAATTAACAAAATGTGGTGGAGC‑ATGGGAGCCCGGCAGGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAAGCTCGCGCATGTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTGTGGTGTGGAACATGTAGTGCTTAAGGTCGa                                                                                 <  2:401209/216‑1 (MQ=255)
                                                                                                                                            cggcgTTGA‑GATTTTCTGGA‑AG‑GTCGTGGCGGCTTACCGGAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGACCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGa                                                                                 <  1:666222/229‑1 (MQ=255)
                                                                                                                                             ggcgTTGAGGATTTTCT‑GC‑AG‑GTCGTGGGGGCTTACCGCAATTTCGTGCAGCAAAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGaa                                                                                <  1:496109/229‑1 (MQ=255)
                                                                                                                                               cgTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGGTAAGTCTTCCGGTTTGGTGTGGAACATGTAAGGCTTAAGGTCGacat                                                                              >  1:186094/1‑226 (MQ=255)
                                                                                                                                                        ttttCT‑GG‑AG‑GTGGTGGCGGCTTCCCGCAATTTCGGGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCAGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGc                                                                      <  2:186094/228‑1 (MQ=255)
                                                                                                                                                         tttCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGGAGAAGTAAATAAATCTGGCGGAGC‑CGGGGAGCTCCGCCAGGGCCCTTAAACAGCTGGCATATTGCGCCCGGAATAACTCTCTGGCATTTCTTTCCGCAGCGCAGCGGGAATTTCCTGTCGCCTGAGGTCGGTTAAGTCGTCCGGTTTGGTGGGGAACCTTATGTGCTTAAGGTCGAACACCTTCAGCaa                                                                    >  2:540287/1‑229 (MQ=255)
                                                                                                                                                         tttCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGGTTGGTGGGGAACATGTGATGCTTAAGGTGGAACCCTTTacacaa                                                                    >  2:499560/1‑224 (MQ=255)
                                                                                                                                                          tttttttggtgtgcggggcgtcTCACGGTAATTTCGTGCAGCAGAAGTAAAAAAAGCGGGGGGCGG‑CTGGGAGCCCCGCACGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAg                                                                       <  2:688642/212‑1 (MQ=255)
                                                                                                                                                               cgtgggggcgggggggcTTCCCCAAATTTCGTTCAGCAGAAGAAAAGAACT‑TGGGGGGGCTCTGGAGGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGTCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACAt                                                                 <  2:647616/215‑1 (MQ=255)
                                                                                                                                                                          ggggcgttacggatatgccGTGCAGGAGAAGAAAAGAATTGGGGGGGAGC‑CGGGGAGCCGCGCAGGAGCCGTTAAACAGCTGGCATAGTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGa                                                     <  2:183634/211‑1 (MQ=255)
                                                                                                                                                                          ggcggcTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAGATCTTTAAGCAACATCTTG‑GTATgg                                                      >  1:602252/1‑228 (MQ=255)
                                                                                                                                                                          ggcgcgttcccGCAATTTGGGTCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGCGCCCCGCCAGGGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGa                                                     <  1:72759/220‑1 (MQ=255)
                                                                                                                                                                                ttACCGCAATTTCGTGCAGCAGAAGGAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCCTATTGGGACCGGAATAAATCTCGCGGATTTCTTTTCCCAGCGCAACGGATATTTCCTGCGCCCCGCCGTGGGTTTAGTCTTCCGGCTTGGGGGGGGAAAGTGATGGCTGTAGCTCGGACTCTTTTATGAAAGACTTG‑GGATGGAAGatat                                               >  2:497355/1‑229 (MQ=255)
                                                                                                                                                                                         ttgcggccagaagAAGTAAATAAATCTGGCGGAGC‑TGGGAACCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGa                                       <  1:108222/221‑1 (MQ=255)
                                                                                                                                                                                              tGCAGCAGAAGTAAATAAATCTGTCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATATCGCGCGTTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAATCTTCCGGATTGGTGTGGAACATGTAGGGCTTATGGTCGAACTCTTTAAGCAACATCTTT‑GTTTGGTAGGTGTTTTCCTGCTAGAc                                  >  1:728832/1‑228 (MQ=255)
                                                                                                                                                                                                  gcagAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCACCGCCAGAGCCGTTAAACAGCTGGCATAGTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTca                             <  1:87094/229‑1 (MQ=255)
                                                                                                                                                                                                      aaGTAAATAATTGTCGCCGAGC‑TGGGGATCCCCGCCCGAGCCGTTAAACAGCTTGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCACCATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTCACATc                         <  1:702844/229‑1 (MQ=255)
                                                                                                                                                                                                        gtaaatatgttggggggggC‑CTGGGGGCCCCCCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTCACATcca                       <  2:195196/211‑1 (MQ=255)
                                                                                                                                                                                                         taaataatGCGGGCGGTGC‑CTGGGATCTCCCCCAGAGCCTTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGACCTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTCACATcca                       <  1:497355/228‑1 (MQ=255)
                                                                                                                                                                                                         taaataaaTCGGGCGGATC‑CTGGGAGCTCCGCCGGGGCCGTTAAACAGCGGGCATATTGCGCCCGTAATAAAGCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTCACATcca                       <  1:399512/228‑1 (MQ=255)
                                                                                                                                                                                                         taaataaaGTTGGGGGAGC‑CGGGGAACTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCACGTAATAAAGCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTCACATcca                       <  1:609046/228‑1 (MQ=255)
                                                                                                                                                                                                         gagacaggtCTGGCGGTGC‑CTGGGTGCTCCGCCAGGGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTCACATcca                       <  1:409309/220‑1 (MQ=255)
                                                                                                                                                                                                          aaataattTGGGGGGGGC‑CTGGGCGCCCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTCACATccacc                     <  2:446867/229‑1 (MQ=255)
                                                                                                                                                                                                           aataatgtgggCGGGGC‑TCGGGAGCCCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTCACATccacc                     <  2:423893/219‑1 (MQ=255)
                                                                                                                                                                                                                    ggCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCAGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGGAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGGTTGGACCATGTAGGGCTTAAGGTCGAACTCTTTAAGCAAC‑TATTGTGGATGGAAGAAATTTTCCCGTTCGACGtctccacacaccacctcatcaa           >  2:494845/1‑208 (MQ=255)
                                                                                                                                                                                                                     ggggggC‑CTGGGGGCCCCCCCAGGGCCCTTAATCAGCGGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTCACATCCACCATGTCATACAg          <  2:640268/224‑1 (MQ=255)
                                                                                                                                                                                                                           c‑cTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGCAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGGAGATATTTTCCTGTAAGACGTTCACATCTATCAttccactctggcctct    >  2:409309/1‑214 (MQ=255)
                                                                                                                                                                                                                              tggggtgtCCGCGAGAGCCGTTAAACAGCTGGCACATTGCCCCCGTAATAAATCCCGCGGGTTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTCACATCCACCATGTCATACAGCGCCTTCa  <  2:532121/222‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              
TCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCAACTGGAATTCGCTGGCGATGAAATGGAAGAGGGGATGAACTACTACGACGAAGAAAGCCTGTCGCTATGCGGCGTTGAGGATTTTCT‑GC‑AG‑GTCGTGGCGGCTTACCGCAATTTCGTGCAGCAGAAGTAAATAAATCTGGCGGAGC‑CTGGGAGCTCCGCCAGAGCCGTTAAACAGCTGGCATATTGCGCCCGTAATAAATCTCGCGCATTTCTTTCCACAGCGCAGCGGTAATTTCCTGGCGCTCGCTGTCGGTTAAGTCTTCCGGTTTGGTGTGGAACATGTAGTGCTTAAGGTCGAACTCTTTAAGCAACATCTTG‑GTATGGAAGATATTTTCCTGATAGACGTTCACATCCACCATGTCATACAGCGCCTTCA  >  NC_000913/134550‑134996

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: ‑