Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 470,035 A→G 7.7% S389G (AGC→GGC)  mdlA → putative multidrug ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913470,0350AG 7.7% 151.0 / 3.0 52S389G (AGC→GGC) mdlAputative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (28/18);  new base G (2/2);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

TAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTT  >  NC_000913/469764‑470289
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                              
tAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGc                                                                                                                                                                                                                                                              >  1:32016/1‑274 (MQ=255)
  ccaatgtagaagttGGTAGGGCTGCGTACAGCCGTACTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAg                                                                                                                                                                                                                                                           <  1:160436/262‑1 (MQ=255)
           aaCGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATAtt                                                                                                                                                                                                                                                  <  2:32016/275‑1 (MQ=255)
             cgggggagggcggcGTACAGCCGTATTCGGGGGATGCTGGCGGAAGGGCCGGTGGGGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCg                                                                                                                                                                                                                                                <  2:204505/263‑1 (MQ=255)
             cGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGCTGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCATGAAACCCGGTCCGATGCTGGGTATCTGCGGGCAGACTGCTTCCGGCAACAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACTTCCGCGAGGGGGATATTc                                                                                                                                                                                                                                                 >  1:84649/1‑274 (MQ=255)
                ggTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGGGGGGGagtgtcgct                                                                                                                                                                                                                                              >  2:34289/1‑266 (MQ=255)
                     tgctgcGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCCGACTGACCATCATGCGCTGGAAAACGCCAATTTCGCCCTGAAAACCGGTCACATGCTGGGTATCTCCCGGCCGACGGGGTCCGGCCAAAGTACCCTCTTGTCGCCCATTCAGCGTCCTTTCGACGTTATCGAGGGGGTATTTCGGTTTCCTg                                                                                                                                                                                                                                        >  2:309259/1‑275 (MQ=255)
                       ctgcGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGCCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCCCTCATTCCGCGTCATTTCGCCGTCAGGGCGGGGGACTTTCGCGTTCATGat                                                                                                                                                                                                                                      >  2:325450/1‑275 (MQ=255)
                          cGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGCTATTCGCTTTCATGATAtt                                                                                                                                                                                                                                   >  1:71317/1‑275 (MQ=255)
                          cGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATAGTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCCGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACACTGTTGTCGCCCATTACGCGGCATTTCGACGTCAGCGAGGGGGGTATTCGCTTTCAGGAtttt                                                                                                                                                                                                                                   >  1:251596/1‑275 (MQ=255)
                               aGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCCAACGTACCCTGTTGTCGCTCCTTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTTCCTTTCATGATATTCctt                                                                                                                                                                                                                               >  2:134778/1‑273 (MQ=255)
                                           gcgATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGGTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGGGGGGGAGATTCGCTTTCATGATATTCCTCCGACGAAGTTAc                                                                                                                                                                                                                   >  2:262593/1‑274 (MQ=255)
                                                         aGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATTTTCCTCTTGCGAAGTTACAACTCGATAgtcgcc                                                                                                                                                                                                    >  2:484650/1‑270 (MQ=255)
                                                         aGCGACGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAAACCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCAGGCAAAAGTACCATGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGTGAGGGGGCTATTCGCTTTCATGATATTCCTCTGACGAAGTTACACCGCGAGAGCt                                                                                                                                                                                                       >  1:54616/1‑272 (MQ=255)
                                                                                        gTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGa                                                                                                                                                                         >  2:100328/1‑271 (MQ=255)
                                                                                                   gCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCtttttt                                                                                                                                                          >  1:160490/1‑275 (MQ=255)
                                                                                                          gAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGaca                                                                                                                                                    <  1:100328/274‑1 (MQ=255)
                                                                                                             cTGGATGTAAATATTCCCCCGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTAAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACtg                                                                                                                                                 <  2:498384/274‑1 (MQ=255)
                                                                                                             cTGGAGTTAAATACTCACCAGTTCACGTATCCGCAGACTGACCACCCTGCGCTGGAAAACGTCTATTGCGCCCTGAAGCCCGTTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGTGAGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGTTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACtg                                                                                                                                                 <  1:313837/274‑1 (MQ=255)
                                                                                                              tGGGTGTAAAGATTCCCCAGTTCACGTATCCGCAGGCTGACCAACCTGCGCTGGAAAACGGCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACtgtg                                                                                                                                               <  2:71317/275‑1 (MQ=255)
                                                                                                                          ttCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGCTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCTTGATATTCCTCTGACGACGTTACCACTCGATAGCTGGCGTAGCCTCCTGGCGGTCGTTCGCCCGACGCCATTCCTTTTTTCCGACCCTGTGGCGAATAACAt                                                                                                                                    >  2:165541/1‑274 (MQ=255)
                                                                                                                            cACCATGTCACGTATCCGCTGACTGACCATCCTGGGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCATATGCTGGGTATCTGCGTGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATcgc                                                                                                                                 <  2:54616/275‑1 (MQ=255)
                                                                                                                                        tATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGccc                                                                                                                     <  1:262593/275‑1 (MQ=255)
                                                                                                                                                 aCTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCCCGCTGGGTTGCCCGAATGCCAc                                                                                                            >  2:230151/1‑275 (MQ=255)
                                                                                                                                                  gtGACCATCCTGCGCTGGAAAACTTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCAcc                                                                                                           <  2:84649/274‑1 (MQ=255)
                                                                                                                                                  cTGACCATCCGGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGGCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCAGGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCAcc                                                                                                           <  2:113877/275‑1 (MQ=255)
                                                                                                                                                     aCCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGCCACTCTGGCGAATAACATCGCGCTGGGTTGCCCGAGTGCCAccacg                                                                                                        >  2:317265/1‑272 (MQ=255)
                                                                                                                                                     aCCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCCTGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACCCTGGGGCGAATAAAACCCGCCTGGGTCTCCCGACTGcccccccg                                                                                                        >  2:275943/1‑273 (MQ=255)
                                                                                                                                                      ccATCCCGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATGCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGc                                                                                                       <  2:166781/275‑1 (MQ=255)
                                                                                                                                                           cTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTTTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCagggg                                                                                                  >  1:401691/1‑271 (MQ=255)
                                                                                                                                                           cTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTTGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAagag                                                                                                  >  2:201126/1‑275 (MQ=255)
                                                                                                                                                                             aTTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACTTCAGCGAGGGGGATATTCGCATTCATGATCTTCCTCTGCCGAAGTTACAACTCGATAGCTGGCGTAGCCGACTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCGTACACTGTGGCGCATAACATCGCGCGTGGTTGCCCGACTCCCACCCACAAAGAGATTGAGCCTGACGCGCgg                                                                                >  1:474163/1‑275 (MQ=255)
                                                                                                                                                                             aTTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGCCGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGCCGCCATTCCTTTTTTTTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGAATGTCGCGCgg                                                                                >  2:134328/1‑275 (MQ=255)
                                                                                                                                                                                                    aTGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCTAATAACATCCCGCTGGTCTGCCCGAATGCCCCCCCACCAGAGCTTGAGACGGTCGCGGCGTTCGCCCGCGTACATGACGaaat                                                         >  1:68962/1‑275 (MQ=255)
                                                                                                                                                                                                          ggTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCAt                                                                                         >  2:456088/1‑237 (MQ=255)
                                                                                                                                                                                                              tCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCGGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCtc                                               >  2:144990/1‑274 (MQ=255)
                                                                                                                                                                                                              tCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTcg                                                                                     <  1:456089/237‑1 (MQ=255)
                                                                                                                                                                                                              tCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTa                                               >  1:128418/1‑275 (MQ=255)
                                                                                                                                                                                                              tCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTa                                               >  1:142616/1‑275 (MQ=255)
                                                                                                                                                                                                               cTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTCGTTAGCCAGACGCAATTCCTTTTTTCTGACACTGTGGCGAATAACATCGGCCTGGATTGCCGGAATTCCACCGAGCCAGAGATAAAGAATGTCGCGCGGTTGGGAAGCGTACTTGACGACGTTTTGCGTCTa                                               >  1:178485/1‑274 (MQ=255)
                                                                                                                                                                                                                   gggcggcattgTTCCTGCAAAAGTACCCTTTTGTCGCTGCTTCAGCGTCATTTCGACGTGCGCGAGGGGGATATTCGCTGTCAGGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGc                                           <  2:496160/264‑1 (MQ=255)
                                                                                                                                                                                                                   gggCCGTGTGGGTACGGCAGAAGTACCCTGTTGTCGCTCATTCGGCGCCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCAGGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGc                                           <  2:343810/274‑1 (MQ=255)
                                                                                                                                                                                                                   gggCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCGGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGc                                           <  1:144990/274‑1 (MQ=255)
                                                                                                                                                                                                                    ggCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCaa                                         <  1:364972/275‑1 (MQ=255)
                                                                                                                                                                                                                    ggCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCaa                                         >  1:229704/1‑275 (MQ=255)
                                                                                                                                                                                                                        gatgggTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCAGTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGt                                      <  1:134778/270‑1 (MQ=255)
                                                                                                                                                                                                                                ccGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATAc                              <  2:142616/274‑1 (MQ=255)
                                                                                                                                                                                                                                   gCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCGGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGc                                                                                                               >  2:201606/1‑190 (MQ=255)
                                                                                                                                                                                                                                       aaGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCGGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCAcc                                                                                                           <  1:201606/190‑1 (MQ=255)
                                                                                                                                                                                                                                              cTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAgcg               <  2:128418/275‑1 (MQ=255)
                                                                                                                                                                                                                                                      gCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGCAGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGGGTGa       >  2:449542/1‑275 (MQ=255)
                                                                                                                                                                                                                                                           ttCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACCCTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATCTTTTGCGTCTCCCGCCAGGGTACGACACAGCGGTGGGCGAGCGCTGTGGGCTGCtt  >  1:424896/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                              
TAACATTGTGGAACGTGGTAGTGCTGCGTACAGCCGTATTCGCGCGATGCTGGCGGAAGCGCCGGTGGTGAACGATGGTAGCGAACCCGTGCCGGAAGGGCGTGGCGAACTGGATGTAAATATTCACCAGTTCACGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCCAGCGTACATGACGATATTTTGCGTCTACCGCAAGGTTACGATACAGAGGTGGGCGAGCGCGGTGTGATGCTT  >  NC_000913/469764‑470289

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: