Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 482,397 G→A 6.0% P669P (CCC→CCT acrB ← multidrug efflux pump RND permease AcrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913482,3970GA 6.0% 202.2 / 2.4 66P669P (CCC→CCTacrBmultidrug efflux pump RND permease AcrB
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base A (2/2);  total (33/33)
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.94e-01

TACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACG  >  NC_000913/482256‑482546
                                                                                                                                             |                                                                                                                                                     
tACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGtt                                                                                                                                               >  2:1751455/1‑150 (MQ=255)
    cTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAg                                                                                                                                           >  2:1825725/1‑150 (MQ=255)
         cAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGaa                                                                                                                                      >  1:746205/1‑150 (MQ=255)
             atatCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAAcc                                                                                                                                  >  1:1062115/1‑150 (MQ=255)
             atatCAGGGTGCTTCGCGGCTTCTGCAAGCAACTGGTTACGCGGCTGAGTGAGTTTTTCGTGACAAATGCCAGCCTGGTCATTCAGCTCAAAGTCAATGCGGGTTTCACTACCCAGTTCCACGATTGCGGGCAGGTTAATGGCGAAAAcc                                                                                                                                  <  2:1062906/150‑1 (MQ=255)
                 cAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCg                                                                                                                              <  2:1617995/150‑1 (MQ=255)
                 cAGGGTGCTTCGCTGCATCTGCAAGCAACTGGTTACGCGCCTGAGTCAGATTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCg                                                                                                                              >  2:1238815/1‑150 (MQ=255)
                  tgtttGCTTCTCTGGTTTTGCAAGCAACTTGTTACGGGTCTGAGTCAGTTTTTCGTGAGGAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCTGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGc                                                                                                                             <  2:1464942/146‑1 (MQ=255)
                       gCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCtt                                                                                                                        >  1:496278/1‑150 (MQ=255)
                           cGCTGCTTCTGCAAGCAACTTGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGAt                                                                                                                    >  1:35745/1‑150 (MQ=255)
                           cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGAt                                                                                                                    <  1:1825725/150‑1 (MQ=255)
                             ctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGAttt                                                                                                                  <  1:367748/150‑1 (MQ=255)
                                     gcaagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGaaa                                                                                                          >  2:1119267/1‑150 (MQ=255)
                                     gcaagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGaaa                                                                                                          >  1:656584/1‑150 (MQ=255)
                                     gcaaccaaCTGGTTACGCGACTGAGTCAGTTTTTAGTGACCAAGGCCAGCCTGGTCAATCAGATCAAAATCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGaaa                                                                                                          >  2:1742477/1‑150 (MQ=255)
                                      caagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAg                                                                                                         >  2:1629693/1‑150 (MQ=255)
                                      caagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAg                                                                                                         <  2:205437/150‑1 (MQ=255)
                                        agcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAgcg                                                                                                       >  1:1039745/1‑150 (MQ=255)
                                                gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGc                                                                                               <  2:915127/150‑1 (MQ=255)
                                                  tACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCAGGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTTCGAGAAAGCGCGTGTTGCAc                                                                                             <  2:1081876/150‑1 (MQ=255)
                                                    cgcgCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGc                                                                                           >  1:1709908/1‑150 (MQ=255)
                                                      cgcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCAt                                                                                         >  2:508064/1‑150 (MQ=255)
                                                       gcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATg                                                                                        <  1:54359/150‑1 (MQ=255)
                                                              tCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCg                                                                                 <  1:1795120/150‑1 (MQ=255)
                                                                aGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCt                                                                               >  2:1466444/1‑150 (MQ=255)
                                                                      tCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACt                                                                         >  1:825131/1‑150 (MQ=255)
                                                                       cGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAACGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACtt                                                                        >  2:435770/1‑150 (MQ=255)
                                                                       cGTGACAAAGGCAAGCCTGGTCAATCAGATCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACtt                                                                        >  2:863111/1‑150 (MQ=255)
                                                                          gACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGt                                                                     <  2:35745/150‑1 (MQ=255)
                                                                           tccAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGtt                                                                    >  2:1354678/2‑150 (MQ=255)
                                                                            ccAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGttt                                                                   >  2:509420/1‑150 (MQ=255)
                                                                             cAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGtttt                                                                  >  2:571114/1‑150 (MQ=255)
                                                                                 gCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTttcttc                                                              >  2:229318/1‑150 (MQ=255)
                                                                                  ccAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATg                                                                                        <  1:1629796/123‑1 (MQ=255)
                                                                                  ccAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATg                                                                                        >  2:1629796/1‑123 (MQ=255)
                                                                                  ccAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCg                                                             >  1:524861/1‑150 (MQ=255)
                                                                                   cAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGc                                                            <  2:1116404/150‑1 (MQ=255)
                                                                                        tGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGa                                                       <  1:124605/150‑1 (MQ=255)
                                                                                           tCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGc                                                                                                                                         >  2:686071/1‑65 (MQ=255)
                                                                                           tCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGc                                                                                                                                         <  1:686071/65‑1 (MQ=255)
                                                                                              aTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCg                                                 <  2:1478316/150‑1 (MQ=255)
                                                                                              aTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCg                                                 <  2:1066972/150‑1 (MQ=255)
                                                                                                cAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGc                                               <  2:254735/150‑1 (MQ=255)
                                                                                                cAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGc                                               <  1:1629693/150‑1 (MQ=255)
                                                                                                 aGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGcc                                              <  1:509420/150‑1 (MQ=255)
                                                                                                    tcaaagtcaaagCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAg                                           <  2:100013/150‑1 (MQ=255)
                                                                                                              agCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAgg                                 <  2:656584/150‑1 (MQ=255)
                                                                                                              agCCGGTTGCAGTACCCAGTTCCACGATTGCAGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAgg                                 >  2:1936215/1‑150 (MQ=255)
                                                                                                              agCCGGTTGCAGTACCCAGTTCCACGATTGCAGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGACCTTCAAgg                                 >  2:1929310/1‑150 (MQ=255)
                                                                                                               gCCGGTTGCAGTACCCAGTTCCACGATTGCAGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGa                                <  1:1929310/150‑1 (MQ=255)
                                                                                                               gCCGGTTGCAGTACCCAGTTCCACGATTGCAGGCAGGGTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGa                                <  1:1936215/150‑1 (MQ=255)
                                                                                                                 cGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGAGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGaaa                              >  2:114897/1‑150 (MQ=255)
                                                                                                                       cAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGTACGATCGGCACAGTCCTTCAAGTAaactaaag                        <  2:499320/150‑3 (MQ=255)
                                                                                                                       cAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAaacgaacg                        <  1:912704/150‑1 (MQ=255)
                                                                                                                        aGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGc                       <  1:1982394/150‑1 (MQ=255)
                                                                                                                        aGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTAAAGGAAACGAACGc                       <  1:1119267/150‑1 (MQ=255)
                                                                                                                        aGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGAGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGc                       >  2:346457/1‑150 (MQ=255)
                                                                                                                         gTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCa                      >  2:443256/1‑150 (MQ=255)
                                                                                                                            cccAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATa                   >  1:421742/1‑150 (MQ=255)
                                                                                                                            cccAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATa                   >  1:433831/1‑150 (MQ=255)
                                                                                                                             ccAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATAc                  <  2:1709908/150‑1 (MQ=255)
                                                                                                                              cAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATAcc                 >  1:750863/1‑150 (MQ=255)
                                                                                                                                gTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCgg               <  1:1742477/150‑1 (MQ=255)
                                                                                                                                     aCGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGGTATTc          <  1:1466444/150‑1 (MQ=255)
                                                                                                                                     aCGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGGTATTc          <  2:1062115/150‑1 (MQ=255)
                                                                                                                                        aTTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGa       <  2:433831/150‑1 (MQ=255)
                                                                                                                                             gggCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACg  >  2:1094030/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
TACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACG  >  NC_000913/482256‑482546

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: