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

TATTGGCGATACGGTAACCATTGGTAGCTTCTCGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCT  >  NC_000913/489376‑489807
                                                                                                                                                                                                                             |                                                                                                                                                                                                                      
tATTGGCGATACGGTAACCATTGGTAGCTTCTCGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCAAGCGTCTGATCAACTGGTCGTTGTCTGACACTACTACGCGTCTGGGGATCCGTCTCGCGTTGCCCTAGTGCCCCAATCTGGAAAACGTGGGCAAAATGTTCCt                                                                                                                                                                                                                  >  1:250494/1‑228 (MQ=255)
      cGATACGGTAACCATTGGTAGCTTCTCGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGCGATCCGTCTCAGCGTGGCCTATGGCTCCGCTCTGGAAAAAGTGCGTCAAGTGTAACtcacgg                                                                                                                                                                                                             >  2:258152/1‑222 (MQ=255)
         tACGGTAACCATTGGTAGCTTCTCGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTCCTACACGTCTGGTGATCCGTCTCGGCGGGGCCTATGGCCCCGATCTGGAAAAAGGGCGTAAGGGGTTACTGAAggcgg                                                                                                                                                                                                         >  2:80146/1‑227 (MQ=255)
          aCGGTAACCATTGGTAGCTTCTCGGGGAAGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGCCGTTGACTGACATTCATACGCGTCGGGCGACCCGTCCCGGCTTGGCCATGGGCCCCGACCTGGACACAGTGGGTAAAGGGttcattaagtcggga                                                                                                                                                                                                        >  1:325797/1‑215 (MQ=255)
               aaCCATTGGTATCTTCTCGGGGACGGGAAGGAAGATCCGTATTCGTGCGACAACGATTAGCGATTTCGATCGAAAAGAAGTGATCCTCCCGAACAAAGCGTTGGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGa                                                                                                                                                                                                  <  2:237342/228‑1 (MQ=255)
                actaTTGGTAGCTTCGCGGGGACGTGAAGTAAGATCCGTATCCGGGCGACAACTATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAg                                                                                                                                                                                                 <  1:356639/225‑1 (MQ=255)
                aacATTGGTAGCTTCGCGGGGACGGTAAGTAAAATCCGTATTCGGGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAg                                                                                                                                                                                                 <  2:579039/226‑1 (MQ=255)
                  ccgtgggtgttttgtgggaACCGTAAGTAAGATCCGTTATGCTGCTGCGACGGTTTCCGATTTCGATTCCAAAGGAGTGATCATCCCGAACAAAGCGTTTGGTACCGCGCGTCTGATCCACGGGTCGTTGACTGACACTTCTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCa                                                                                                                                                                                               <  1:397222/213‑1 (MQ=255)
                     tgggagtctcgcGGGGAAGGTAAGTAAGAGCCGGATGCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCAccc                                                                                                                                                                                            <  2:131336/217‑1 (MQ=255)
                      ggTAGCTTCTCGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCa                                                                                                                                                                                           >  2:251821/1‑228 (MQ=255)
                       gTAGCTTCTCGGGGACGGTAAGTAAGGTCCGTAGGCGGGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATGCCGAACAAAGCGTTTGATACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCaa                                                                                                                                                                                          <  2:251894/228‑1 (MQ=255)
                             tctcGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGACCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTAGGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCACGGGGGa                                                                                                                                                                                   >  2:252363/1‑228 (MQ=255)
                                        gtaagtaaGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGCTGCCCATGGCCCCGATCTGGAAAACGTGCGGAAAGGATTCCTGAAGGCGGCGACCGGCAccccccgggggatgcaccaccc                                                                                                                                                                         >  1:251574/1‑208 (MQ=255)
                                                caTCCGTTTTCGTGCGAAAACGATTCCCGATTTCGATCGCAAAGAAGTGATCACCCCGAACAAAGCGTTTGTTACCGGGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCgg                                                                                                                                                                <  1:251821/226‑1 (MQ=255)
                                                 aTCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGTTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGTAGGCGGCGCCTGAGCACACAAGGGTGATGCACGAACCAATGCCGGaa                                                                                                                                                              >  2:159617/1‑228 (MQ=255)
                                                 aTCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGGGCGTAAAGGGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGaa                                                                                                                                                              >  2:467551/1‑228 (MQ=255)
                                                   ccGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGTAAAAGGTGCGTAAAGTGGTACTGATGGGGGCGACTGAGCACCCCAGGGTGTTGCATCAACCacaacgagcggt                                                                                                                                                             >  1:106202/1‑217 (MQ=255)
                                                   ccGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGAAAAGTGTTACTGAAGGCGGCGCCTGAGCACCCAAGGGTGATGCACGACCCACTGCCGGCAg                                                                                                                                                             >  2:388029/1‑227 (MQ=255)
                                                      tATTCGTGCGACAACGATTCCCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTg                                                                                                                                                                                                                <  1:480080/175‑1 (MQ=255)
                                                      tATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATAACCTGGCCGTTGACTGACACTACTACGCGTCTGGTGGTCCGTCTCGGCGTGGCCTATGGCTCCGGTCTGGGAAAAGGGCGTAAGGTGTTAGTCAAGGCGTCGACTGAGGAGCCACAGGCGATGCAAGAACCAATTCCGGaatgctc                                                                                                                                                         >  1:612164/1‑223 (MQ=255)
                                                        ttCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGCCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGGGTAAAGTGTTACTGAAGGCGGCGCCTGAGACCCCAAGGGGGATGCACGAACAACTGCGGGagatctcct                                                                                                                                                       >  1:462914/1‑220 (MQ=255)
                                                            tGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGCCTCCTTTa                                                                                                                                                    >  2:327219/1‑227 (MQ=255)
                                                            tGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACTGAAGGCGGCGACTGAGCACCCAAGGGGGATGCACGACCCAATGCCGGAGGCCTTCtttc                                                                                                                                                    >  2:3146/1‑226 (MQ=255)
                                                             gCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTAc                                                                                                                                                   >  2:75835/1‑227 (MQ=255)
                                                                aCAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGCCTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACgccc                                                                                                                                               >  2:234553/1‑225 (MQ=255)
                                                                aCAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACTGAAGGCGGCGACTGAGCACCCAAGGGGGATGCCCGAACAAATGCCGGAGGTCTTCTTTAAGGc                                                                                                                                                >  1:609099/1‑227 (MQ=255)
                                                                  aaggatgtgcGATTTCGATCGCAAAGAAGTGATCGTCCCGAACAAAGGGGTTGGTACCGAGCGTTTGATCAACTTGTCGTTGACTGACACTGCTACGCGTCGGGTGATCCGTATCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCAt                                                                                                                                              <  1:39276/218‑1 (MQ=255)
                                                                         aCCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACTGAAGGCGGCGACTGAGCCCCCAAGGGTGGTGCCCGAACCAATGCCGGAAGGCTTCTTTACGGAATTTGGTGc                                                                                                                                       >  1:556032/1‑227 (MQ=255)
                                                                          ccGATTTCGAACGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCa                                                                                                                                     <  2:251574/228‑1 (MQ=255)
                                                                               ttCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCCTCGATCCGGAAAAAGGGGGTAAATCGTTCCGGAAGCGGGGGCATGAGGCCCACAGGGGGCGGCCGAACCAAAGACGGGAGGTCTTCTTCACGCGATTTGGtcgcacgca                                                                                                                                 >  1:32125/1‑219 (MQ=255)
                                                                                 tgATCGAAAAGAAGTGATCAGCCCGAGCAAGGCTTTTTTTGCGGGGCGTCTGATCAACTGGTCGTTGAATGACACTACGACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGtt                                                                                                                             <  1:476011/227‑1 (MQ=255)
                                                                                 gcagcGCAAAGAGTTGATCACCCCGAAAAAAGCGTTTTGTACCGGGCGTTTGATTAACTGGGCGTTGACTGACACCTCTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGTCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGtt                                                                                                                             <  1:192128/224‑1 (MQ=255)
                                                                                  gaacGCAAAGGAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGtt                                                                                                                             <  2:462914/224‑1 (MQ=255)
                                                                                   aTCGCAAAGAAGGGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAACGCGGCGACTGAGCACCCAAggg                                                                                                                                                                                       >  2:534719/1‑171 (MQ=255)
                                                                                       gaaaaaaTTGATCACCCCGAAAAAAGGTTTTGGTTCCGGGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTTAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCa                                                                                                                       <  1:126707/227‑1 (MQ=255)
                                                                                           gaaggtatcatcCCGAACAAAGGTTTTGTTACCGGGCGTTTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAg                                                                                                                   <  1:603080/222‑1 (MQ=255)
                                                                                           gaagTGGTCCCCCCGAACAGAGGTTTTGTTACCGGGGGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAg                                                                                                                   <  1:634367/228‑1 (MQ=255)
                                                                                            aaggggtcctcCCCAAAAAGGCGTTTGTTTCCGGGGTTTCGATCAACTGGTCGGTGACTGACTCTACTACGCGTTTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGAGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGc                                                                                                                  <  1:356086/219‑1 (MQ=255)
                                                                                             tgTGATTCTCCCGAAGAAAGCTTTTTTTACCGGGGCTCTGAGCAACGGGTCGTTGGCTGACACTGCTACGGAGCTGGTGATCCGGGTCAGCGTGGCCTAGGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGcg                                                                                                                 <  1:258152/227‑2 (MQ=255)
                                                                                             agggATCATCCCGAACAAAGCGTTGGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGGGTCTGGTGATGCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCt                                                                                                                 <  2:224261/225‑1 (MQ=255)
                                                                                                      cccGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGGCTTCTTTACGGCATTTGGTGCCAGACCGTTGGCTCATGAGCTGGGGCTg                                                                                                          >  2:188378/1‑226 (MQ=255)
                                                                                                      cccGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGCGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACTGAAGGCGGCGACTGAGCACCCAAGGGGGATGCACGAACCAATGCCGGAAGGCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGccactg                                                                                                          >  2:243713/1‑221 (MQ=255)
                                                                                                        ctacacaaGCTTTTGTTCCCGAGCGTCGGATCAACTGGTGGTTGACTGACCCTGCTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCCCCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTCACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTa                                                                                                        <  1:80146/220‑1 (MQ=255)
                                                                                                         gAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACTGAAGGCGGCGCATGAGCACCCAAGGGGGATGCACGAACCAATGCCGGAAGGCTTCTTTCCGGCATTTGGTGCCACGCAGTGGGATAATGAACTGCGTGTGTAtgt                                                                                                     >  1:276863/1‑228 (MQ=255)
                                                                                                              aaGGGTTTGTTACCGGGCGTCTGTTCAACTGGTCGTTGACGGACACTACGACGCGTCTGGTGATCCGTCGCGGGGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCCGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTg                                                                                                <  2:283333/228‑1 (MQ=255)
                                                                                                              aaGCTTTTGTTACCGGGCGTCTGATCAACTGGTCGTTGACTGACCCTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTg                                                                                                <  1:467551/228‑1 (MQ=255)
                                                                                                              aaGCGTTTGTTACCGGGCGTCTGAGCAACGGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTg                                                                                                <  2:504409/228‑1 (MQ=255)
                                                                                                               aGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAATTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGGGATGCACGAACCAATGCCGGAAGGCTTCTTTACGGCATTTGGTGCCAGCCCGTTTGACAATGAGcgtgctccgtattggcgtga                                                                                               >  1:214629/1‑209 (MQ=255)
                                                                                                               aGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCCGCGACTGAGCACCCACGGGTGAAGCACGAACCAATGCCGGAAGCCTTCTTTACGGCCTTTGGTGCCACCCCGTTGGAACATCAgcggcggtcttgggtgcgtgg                                                                                               >  1:602693/1‑213 (MQ=255)
                                                                                                                    ttGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCCCCCAAGGGTGATGCACGAACCAATTCCGGAAGCCGTCTTTAGCGCATTTGGTGCCAGCCCGGTGGTACATTAGGCGCGCCGTTATGTGCGTGAAGTGc                                                                                          >  1:151768/1‑228 (MQ=255)
                                                                                                                      gTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGt                                                                                        <  2:76490/228‑1 (MQ=255)
                                                                                                                             aGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGTACTGCGTGCCCgga                                                                                 >  2:513924/1‑226 (MQ=255)
                                                                                                                               cGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCATATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGGGATCGACGAACCAATGTCGGAAGCCTTCTTTACGGCATTCGGGGCCCGCACGTTGGCTCATGAGCTCAGTCGGCATTGGCGGTAACTGCGTGAcaggaag                                                                               >  1:187595/1‑222 (MQ=255)
                                                                                                                                    gATCAACTGGGTGTTGACTGACACTACTCCGCGTCTGGTGATCCGTCTCGGCGTGGACTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATGTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTag                                                                          <  2:572909/228‑2 (MQ=255)
                                                                                                                                          ctgtgcGTTGCCTGACAATAATAAGCGTCTGGTGGTCCGGCCCGGCGTGGCCTATGGCTACGATCTGGAAAAAGTGCGAAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGAATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCg                                                                     <  2:106202/222‑1 (MQ=255)
                                                                                                                                               cGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGGGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAAATGCGTGACCGTAGTCGTACTGTCCGTGAGc                                                               >  2:281974/1‑228 (MQ=255)
                                                                                                                                                  tgagtggcacacttactgtgcTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGTCGGCGACTGAGCCCCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTGGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                                             <  2:550712/207‑1 (MQ=255)
                                                                                                                                                       gacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGa                                                                                                                                                               >  2:358779/1‑125 (MQ=255)
                                                                                                                                                        acacTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGaa                                                                                                                                                              <  1:358779/125‑1 (MQ=255)
                                                                                                                                                        acacTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGCCTTCTTTACGGCATTTGGGGCCAGCAGTTTGGATCATGAACTGGGTCTTTATGGCGGTGAACTGCGTGACGGGAGTCGTCGTGTCgtttggcggacccgga                                                      >  1:495167/1‑213 (MQ=255)
                                                                                                                                                         cattactccGCGTCTGGTGATCCGTCTGGGCGGGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTAc                                                     <  1:234553/225‑1 (MQ=255)
                                                                                                                                                           gtactacGCGTCTGGGGATCCGTCTCGGCGTGGCATATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTa                                                   <  2:187595/227‑1 (MQ=255)
                                                                                                                                                              ctcggcgTTTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGGTCTGGAAAAAGGGCGGAAAGTGTTACGGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCg                                                <  1:252363/224‑1 (MQ=255)
                                                                                                                                                              ctacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTg                                                                                                                                           >  2:118544/1‑138 (MQ=255)
                                                                                                                                                              ctacGCGTCTGGTGATCCGTCGCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGTACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCg                                                <  2:32125/228‑1 (MQ=255)
                                                                                                                                                               tacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTgg                                                                                                                                          <  1:118544/138‑1 (MQ=255)
                                                                                                                                                                              ccGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGGGACCGTAGTCGTACTGTCGATGAGCTGAACCGGACTATCGATAACATGTGCCGTGa                                >  2:472512/1‑228 (MQ=255)
                                                                                                                                                                                gactcGGCGGTGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAAGGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaa                              <  1:75835/226‑1 (MQ=255)
                                                                                                                                                                                 ttgcgggggggCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGGTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa                             <  2:609099/219‑1 (MQ=255)
                                                                                                                                                                                 tcgcggcggGGCCTATGGGTGCGATCGGGAAAAAGTTCGTAAAGTGGTACTGAAGGCGGCGTCTGAGCACGCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCGTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGGGaaa                              <  2:87456/224‑1 (MQ=255)
                                                                                                                                                                                   tcGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGGGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGGGACCGTAGTCGTACTGTGGATGAGCTGAACCGTa                                                      >  2:52881/1‑201 (MQ=255)
                                                                                                                                                                                       gcgtgCCTAGGGCTCCGAACTGGAAAAAGGTGGTAAAGGGTTACTGAAGGGGGCGACTGGGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACAt                       <  2:571952/224‑1 (MQ=255)
                                                                                                                                                                                           gCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCaa                    >  2:120392/1‑227 (MQ=255)
                                                                                                                                                                                             ctaatgCTCCGATCTTGAAAAAGTGCTTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGAATTTGGTGCCTGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACAt                 <  2:276863/223‑1 (MQ=255)
                                                                                                                                                                                                      cGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTaa        >  2:435790/1‑228 (MQ=255)
                                                                                                                                                                                                      cGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTaa        >  2:219019/1‑228 (MQ=255)
                                                                                                                                                                                                      cGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACAATCGTTCAGCTGTGCCGTGAAACCGCCATCACCATTGCCTTTaa        >  2:644917/1‑228 (MQ=255)
                                                                                                                                                                                                      cGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGAACTATCGATCAGCTTTCCCTGGAAAACGAGATCACAATTGCCTTTaa        >  2:371195/1‑228 (MQ=255)
                                                                                                                                                                                                       gATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAg                                                                <  1:657509/171‑1 (MQ=255)
                                                                                                                                                                                                            ggAAAAAGTGGCTAAAGGTTTACTGGAGGGCGGAATCTGACCGCAAAGGGTGGATGCCGGAGCAATTCCGGAAGTCTGCTTTACGGGATTTGGTGCCAGGACGTTGGATCATGAGCTGCGTCTGTATGTGCGGGAACTGCGTGACCGTGGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCt  <  1:132936/228‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                      
TATTGGCGATACGGTAACCATTGGTAGCTTCTCGGGGACGGTAAGTAAGATCCGTATTCGTGCGACAACGATTACCGATTTCGATCGCAAAGAAGTGATCATCCCGAACAAAGCGTTTGTTACCGAGCGTCTGATCAACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATCAACATTGCCTTTAACCAGCT  >  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: