Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 489,597 T→G 14.2% V1021G (GTG→GGG)  mscK → mechanosensitive channel protein, intermediate conductance, K+ regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913489,5970TG14.2% 163.2 / 2.4 80V1021G (GTG→GGG) mscKmechanosensitive channel protein, intermediate conductance, K+ regulated
Reads supporting (aligned to +/- strand):  ref base T (28/38);  new base G (11/2);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGGCGATACGGTAACCATTGGTAGCTTCTCGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCT  >  NC_000913/489376‑489807
                                                                                                                                                                                                                             |                                                                                                                                                                                                                      
tATTGGCGATACGGTAACCATTGGTAGCTTCTCGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCAAGCGTCTGATCAACTGGTCGTTGTCTGACACTACTACGCGTCTGGGGATCCGTCTCGCGTTGCCCTAGTGCCCCAATCTGGAAAACGTGGGCAAAATGTTCCt                                                                                                                                                                                                                  >  1:250494/1‑228 (MQ=255)
      cGATACGGTAACCATTGGTAGCTTCTCGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGCGATCCGTCTCAGCGTGGCCTATGGCTCCGCTCTGGAAAAAGTGCGTCAAGTGTAACtcacgg                                                                                                                                                                                                             >  2:258152/1‑222 (MQ=255)
         tACGGTAACCATTGGTAGCTTCTCGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTCCTACACGTCTGGTGATCCGTCTCGGCGGGGCCTATGGCCCCGATCTGGAAAAAGGGCGTAAGGGGTTACT‑GAAggcgg                                                                                                                                                                                                         >  2:80146/1‑227 (MQ=255)
          aCGGTAACCATTGGTAGCTTCTCGGGGAAGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGCCGTTGACTGACATTCATACGCGTCGGGCGACCCGTCCCGGCTTGGCCATGGGCCCCGACCTGGACACAGTGGGTAAAGGGttcattaagtcggga                                                                                                                                                                                                        >  1:325797/1‑215 (MQ=255)
               aaCCATTGGTATCTTCTCGGGGACGGGAAGGAAGATCCGTATTCGTGCGACAACGATTAGCGATTTCGATCGAAAAGAAGTGATCCTCCCGAACAAAGCGTTGGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGa                                                                                                                                                                                                  <  2:237342/228‑1 (MQ=255)
                actaTTGGTAGCTTCGCGGGGACGTGAAGTAAGATCCGTATCCGGGCGACAACTATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAg                                                                                                                                                                                                 <  1:356639/225‑1 (MQ=255)
                aacATTGGTAGCTTCGCGGGGACGGTAAGTAAAATCCGTATTCGGGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAg                                                                                                                                                                                                 <  2:579039/226‑1 (MQ=255)
                  ccgtgggtgttttgtgggaACCGTAAGTAAGATCCGTTATGCTGCTGCGACGGTTTCCGATTTCGATTCCAAAGGAGTGATCATCCCGAACAAAGCGTTTGGTACCGCGCGTCTGATCCACGGGTCGTTGACTGACACTTCTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCa                                                                                                                                                                                               <  1:397222/213‑1 (MQ=255)
                     tgggagtctcgcGGGGAAGGTAAGTAAGAGCCGGATGCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCAccc                                                                                                                                                                                            <  2:131336/217‑1 (MQ=255)
                      ggTAGCTTCTCGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCa                                                                                                                                                                                           >  2:251821/1‑228 (MQ=255)
                       gTAGCTTCTCGGGGACGGTAAGTAAGGTCCGTAGGCGGGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATGCCGAACAAAGCGTTTGATACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCaa                                                                                                                                                                                          <  2:251894/228‑1 (MQ=255)
                             tctcGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGACCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTAGGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCACGGGG‑Ga                                                                                                                                                                                   >  2:252363/1‑228 (MQ=255)
                                        gtaagtaaGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGCTGCCCATGGCCCCGATCTGGAAAACGTGCGGAAAGGATTCCT‑GAAGGCGGCGACCG‑GCAccccccgggggatgcaccaccc                                                                                                                                                                         >  1:251574/1‑208 (MQ=255)
                                                caTCCGTTTTCGTGCGAAAACGATTCCCGATTTCGATCGCAAAGAAGTGATCACCCCGAACAAAGCGTTTGTTACCGGGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCgg                                                                                                                                                                <  1:251821/226‑1 (MQ=255)
                                                 aTCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGTTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GTAGGCGGCGCCTGAGCACACAAGGGT‑GATGCACGAACCAATG‑CCGGaa                                                                                                                                                              >  2:159617/1‑228 (MQ=255)
                                                 aTCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGGGCGTAAAGGGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGaa                                                                                                                                                              >  2:467551/1‑228 (MQ=255)
                                                   ccGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGTAAAAGGTGCGTAAAGTGGTACT‑GATGGGGGCGACTGAGCACCCCAGGGT‑GTTGCATCAACCacaacgagcggt                                                                                                                                                             >  1:106202/1‑217 (MQ=255)
                                                   ccGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGAAAAGTGTTACT‑GAAGGCGGCGCCTGAGCACCCAAGGGT‑GATGCACGACCCACTG‑CCGGCAg                                                                                                                                                             >  2:388029/1‑227 (MQ=255)
                                                      tATTCGTGCGACAACGATTCCCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑g                                                                                                                                                                                                                <  1:480080/175‑1 (MQ=255)
                                                      tATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATAACCTGGCCGTTGACTGACACTACTACGCGTCTGGTGGTCCGTCTCGGCGTGGCCTATGGCTCCGGTCTGGGAAAAGGGCGTAAGGTGTTAGT‑CAAGGCGTCGACTGAGGAGCCACAGGC‑GATGCAAGAACCAATT‑CCGGaatgctc                                                                                                                                                         >  1:612164/1‑223 (MQ=255)
                                                        ttCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGCCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGGGTAAAGTGTTACT‑GAAGGCGGCGCCTGAGACCCCAAGGGG‑GATGCACGAACAACTG‑CGGGagatctcct                                                                                                                                                       >  1:462914/1‑220 (MQ=255)
                                                            tGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGCCTCCTTTa                                                                                                                                                    >  2:327219/1‑227 (MQ=255)
                                                            tGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGG‑GATGCACGACCCAATG‑CCGGAGGCCTTCtttc                                                                                                                                                    >  2:3146/1‑226 (MQ=255)
                                                             gCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTAc                                                                                                                                                   >  2:75835/1‑227 (MQ=255)
                                                                aCAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGCCTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACgccc                                                                                                                                               >  2:234553/1‑225 (MQ=255)
                                                                aCAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGG‑GATGCCCGAACAAATG‑CCGGAGGTCTTCTTTAAGGc                                                                                                                                                >  1:609099/1‑227 (MQ=255)
                                                                  aaggatgtgcGATTTCGATCGCAAAGAAGTGATCGTCCCGAACAAAGGGGTTGGTACCGAGCGTTTGATCAACTTGTCGTTGACTGACACTGCTACGCGTCGGGTGATCCGTATCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCAt                                                                                                                                              <  1:39276/218‑1 (MQ=255)
                                                                         aCCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACT‑GAAGGCGGCGACTGAGCCCCCAAGGGT‑GGTGCCCGAACCAATG‑CCGGAAGGCTTCTTTACGGAATTTGGTGc                                                                                                                                       >  1:556032/1‑227 (MQ=255)
                                                                          ccGATTTCGAACGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCa                                                                                                                                     <  2:251574/228‑1 (MQ=255)
                                                                               ttCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCCTCGATCCGGAAAAAGGGGGTAAATCGTTCCG‑GAAGCGGGGGCATGAGGCCCACAGGGG‑GCGGC‑CGAACCAAAGACGGGAGGTCTTCTTCACGCGATTTGGtcgcacgca                                                                                                                                 >  1:32125/1‑219 (MQ=255)
                                                                                 tgATCGAAAAGAAGTGATCAGCCCGAGCAAGGCTTTTTTTGCGGGGCGTCTGATCAACTGGTCGTTGAATGACACTACGACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGtt                                                                                                                             <  1:476011/227‑1 (MQ=255)
                                                                                 gcagcGCAAAGAGTTGATCACCCCGAAAAAAGCGTTTTGTACCGGGCGTTTGATTAACTGGGCGTTGACTGACACCTCTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGTCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGtt                                                                                                                             <  1:192128/224‑1 (MQ=255)
                                                                                  gaacGCAAAGGAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGtt                                                                                                                             <  2:462914/224‑1 (MQ=255)
                                                                                   aTCGCAAAGAAGGGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAACGCGGCGACTGAGCACCCAAggg                                                                                                                                                                                       >  2:534719/1‑171 (MQ=255)
                                                                                       gaaaaaaTTGATCACCCCGAAAAAAGGTTTTGGTTCCGGGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTTAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCa                                                                                                                       <  1:126707/227‑1 (MQ=255)
                                                                                           gaaggtatcatcCCGAACAAAGGTTTTGTTACCGGGCGTTTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAg                                                                                                                   <  1:603080/222‑1 (MQ=255)
                                                                                           gaagTGGTCCCCCCGAACAGAGGTTTTGTTACCGGGGGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAg                                                                                                                   <  1:634367/228‑1 (MQ=255)
                                                                                            aaggggtcctcCCCAAAAAGGCGTTTGTTTCCGGGGTTTCGATCAACTGGTCGGTGACTGACTCTACTACGCGTTTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGAGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGc                                                                                                                  <  1:356086/219‑1 (MQ=255)
                                                                                             tgTGATTCTCCCGAAGAAAGCTTTTTTTACCGGGGCTCTGAGCAACGGGTCGTTGGCTGACACTGCTACGGAGCTGGTGATCCGGGTCAGCGTGGCCTAGGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGcg                                                                                                                 <  1:258152/227‑2 (MQ=255)
                                                                                             agggATCATCCCGAACAAAGCGTTGGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGGGTCTGGTGATGCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCt                                                                                                                 <  2:224261/225‑1 (MQ=255)
                                                                                                      cccGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGGCTTCTTTACGGCATTTGGTGCCA‑GACCGTTGGCTCATGAGCTGGGGCTg                                                                                                          >  2:188378/1‑226 (MQ=255)
                                                                                                      cccGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGCGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGG‑GATGCACGAACCAATG‑CCGGAAGGCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGccactg                                                                                                          >  2:243713/1‑221 (MQ=255)
                                                                                                        ctacacaaGCTTTTGTTCCCGAGCGTCGGATCAACTGGTGGTTGACTGACCCTGCTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCCCCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTCACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTa                                                                                                        <  1:80146/220‑1 (MQ=255)
                                                                                                         gAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACT‑GAAGGCGGCGCATGAGCACCCAAGGGG‑GATGCACGAACCAATG‑CCGGAAGGCTTCTTTCCGGCATTTGGTGCCACGCA‑GTGGGATAATGAACTGCGTGTGTAtgt                                                                                                     >  1:276863/1‑228 (MQ=255)
                                                                                                              aaGGGTTTGTTACCGGGCGTCTGTTCAACTGGTCGTTGACGGACACTACGACGCGTCTGGTGATCCGTCGCGGGGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCCGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTg                                                                                                <  2:283333/228‑1 (MQ=255)
                                                                                                              aaGCTTTTGTTACCGGGCGTCTGATCAACTGGTCGTTGACTGACCCTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTg                                                                                                <  1:467551/228‑1 (MQ=255)
                                                                                                              aaGCGTTTGTTACCGGGCGTCTGAGCAACGGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTg                                                                                                <  2:504409/228‑1 (MQ=255)
                                                                                                               aGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAATTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGG‑GATGCACGAACCAATG‑CCGGAAGGCTTCTTTACGGCATTTGGTGCCA‑GCCCGTTTGACAATGAGcgtgctccgtattggcgtga                                                                                               >  1:214629/1‑209 (MQ=255)
                                                                                                               aGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCCGCGACTGAGCACCCACGGGT‑GAAGCACGAACCAATG‑CCGGAAGCCTTCTTTACGGCCTTTGGTGCCA‑CCCCGTTGGAACATCAgcggcggtcttgggtgcgtgg                                                                                               >  1:602693/1‑213 (MQ=255)
                                                                                                                    ttGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCCCCCAAGGGT‑GATGCACGAACCAATT‑CCGGAAGCCGTCTTTAGCGCATTTGGTGCCA‑GCCCGGTGGTACATTAGGCGCGCCGTTATGTGCGTGAAGTGc                                                                                          >  1:151768/1‑228 (MQ=255)
                                                                                                                      gTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGt                                                                                        <  2:76490/228‑1 (MQ=255)
                                                                                                                             aGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGTACTGCGTGCCCgga                                                                                 >  2:513924/1‑226 (MQ=255)
                                                                                                                               cGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCATATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGG‑GATCGACGAACCAATG‑TCGGAAGCCTTCTTTACGGCATTCGGGGCCC‑GCACGTTGGCTCATGAGCTCAGTCGGCATTGGCGGTAACTGCGTGAcaggaag                                                                               >  1:187595/1‑222 (MQ=255)
                                                                                                                                    gATCAACTGGGTGTTGACTGACACTACTCCGCGTCTGGTGATCCGTCTCGGCGTGGACTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATGTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTag                                                                          <  2:572909/228‑2 (MQ=255)
                                                                                                                                          ctgtgcGTTGCCTGACAATAATAAGCGTCTGGTGGTCCGGCCCGGCGTGGCCTATGGCTACGATCTGGAAAAAGTGCGAAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGAATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCg                                                                     <  2:106202/222‑1 (MQ=255)
                                                                                                                                               cGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAAATGCGTGACCGTAGTCGTACTGTCCGTGAGc                                                               >  2:281974/1‑228 (MQ=255)
                                                                                                                                                  tgagtggcacacttactgtgcTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGTCGGCGACTGAGCCCCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTGGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                                             <  2:550712/207‑1 (MQ=255)
                                                                                                                                                       gacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGa                                                                                                                                                               >  2:358779/1‑125 (MQ=255)
                                                                                                                                                        acacTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGaa                                                                                                                                                              <  1:358779/125‑1 (MQ=255)
                                                                                                                                                        acacTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGCCTTCTTTACGGCATTTGGGGCCA‑GCAGTTTGGATCATGAACTGGGTCTTTATGGCGGTGAACTGCGTGACGGGAGTCGTCGTGTCgtttggcggacccgga                                                      >  1:495167/1‑213 (MQ=255)
                                                                                                                                                         cattactccGCGTCTGGTGATCCGTCTGGGCGGGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTAc                                                     <  1:234553/225‑1 (MQ=255)
                                                                                                                                                           gtactacGCGTCTGGGGATCCGTCTCGGCGTGGCATATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTa                                                   <  2:187595/227‑1 (MQ=255)
                                                                                                                                                              ctcggcgTTTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGGTCTGGAAAAAGGGCGGAAAGTGTTACG‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCg                                                <  1:252363/224‑1 (MQ=255)
                                                                                                                                                              ctacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTg                                                                                                                                           >  2:118544/1‑138 (MQ=255)
                                                                                                                                                              ctacGCGTCTGGTGATCCGTCGCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGTACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCg                                                <  2:32125/228‑1 (MQ=255)
                                                                                                                                                               tacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTgg                                                                                                                                          <  1:118544/138‑1 (MQ=255)
                                                                                                                                                                              ccGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGGGACCGTAGTCGTACTGTCGATGAGCTGAACCGGACTATCGATAACATGTGCCGTGa                                >  2:472512/1‑228 (MQ=255)
                                                                                                                                                                                gactcGGCGGTGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAAGG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaa                              <  1:75835/226‑1 (MQ=255)
                                                                                                                                                                                 ttgcgggggggCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGGTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa                             <  2:609099/219‑1 (MQ=255)
                                                                                                                                                                                 tcgcggcggGGCCTATGGGTGCGATCGGGAAAAAGTTCGTAAAGTGGTACT‑GAAGGCGGCGTCTGAGCACGCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCGTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGGGaaa                              <  2:87456/224‑1 (MQ=255)
                                                                                                                                                                                   tcGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGGGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGGGACCGTAGTCGTACTGTGGATGAGCTGAACCGTa                                                      >  2:52881/1‑201 (MQ=255)
                                                                                                                                                                                       gcgtgCCTAGGGCTCCGAACTGGAAAAAGGTGGTAAAGGGTTACT‑GAAGGGGGCGACTGGGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACAt                       <  2:571952/224‑1 (MQ=255)
                                                                                                                                                                                           gCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCaa                    >  2:120392/1‑227 (MQ=255)
                                                                                                                                                                                             ctaatgCTCCGATCTTGAAAAAGTGCTTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGAATTTGGTGCCT‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACAt                 <  2:276863/223‑1 (MQ=255)
                                                                                                                                                                                                      cGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTaa        >  2:435790/1‑228 (MQ=255)
                                                                                                                                                                                                      cGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTaa        >  2:219019/1‑228 (MQ=255)
                                                                                                                                                                                                      cGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACAATCGTTCAGCTGTGCCGTGAAACCGCCATCACCATTGCCTTTaa        >  2:644917/1‑228 (MQ=255)
                                                                                                                                                                                                      cGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGAACTATCGATCAGCTTTCCCTGGAAAACGAGATCACAATTGCCTTTaa        >  2:371195/1‑228 (MQ=255)
                                                                                                                                                                                                       gATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAg                                                                <  1:657509/171‑1 (MQ=255)
                                                                                                                                                                                                            ggAAAAAGTGGCTAAAGGTTTACTGGAGGGCGGAATCTGA‑CCGCAAAGGGTGGATGC‑CGGAGCAATT‑CCGGAAGTCTGCTTTACGGGATTTGGTGCCA‑GGACGTTGGATCATGAGCTGCGTCTGTATGTGCGGGAACTGCGTGACCGTGGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCt  <  1:132936/228‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                      
TATTGGCGATACGGTAACCATTGGTAGCTTCTCGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACT‑GAAGGCGGCGACTGAGCACCCAAGGGT‑GATGCACGAACCAATG‑CCGGAAGTCTTCTTTACGGCATTTGGTGCCA‑GCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCT  >  NC_000913/489376‑489807

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