Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,667,443 T→G 13.7% T106P (ACC→CCC)  mdtJ ← multidrug efflux system transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,667,4430TG13.7% 699.3 / 117.7 289T106P (ACC→CCC) mdtJmultidrug efflux system transporter
Reads supporting (aligned to +/- strand):  ref base T (117/132);  new base G (26/14);  total (143/146)
Fisher's exact test for biased strand distribution p-value = 4.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

CCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAA  >  CP009273/1667294‑1667590
                                                                                                                                                     |                                                                                                                                                   
 cAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTa                                                                                                                                                    <  1:756354/150‑1 (MQ=255)
  aaaTATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTAc                                                                                                                                                   <  1:2753978/150‑1 (MQ=255)
    atatTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCt                                                                                                                                                 >  1:904099/1‑150 (MQ=255)
    atatTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTGTACGCGCTTTACGGGTACCt                                                                                                                                                 >  1:2007396/1‑150 (MQ=255)
       ttttGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGAt                                                                                                                                              <  2:914373/150‑1 (MQ=255)
       ttttGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGAt                                                                                                                                              <  2:2706661/150‑1 (MQ=255)
         ttGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGAttt                                                                                                                                            <  1:2589080/150‑1 (MQ=255)
           gCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTa                                                                                                                                          <  1:2777496/150‑1 (MQ=255)
            cgacgaAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTAt                                                                                                                                         >  2:1140838/1‑150 (MQ=255)
               cgaAAGCCGTCATAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATACAAACTGCGCCATGGAACACCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCTGAATTTATCaa                                                                                                                                      >  2:302992/1‑150 (MQ=255)
                gaAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCTGATTTTATCAac                                                                                                                                     >  1:4273367/1‑150 (MQ=255)
                    gCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACaa                                                                                                                                 >  2:2048400/1‑150 (MQ=255)
                     ccGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAAt                                                                                                                                <  1:4546621/150‑1 (MQ=255)
                      cGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATc                                                                                                                               <  2:4058975/150‑1 (MQ=255)
                      cGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATc                                                                                                                               >  1:324620/1‑150 (MQ=255)
                        tCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATccc                                                                                                                             <  2:2440984/150‑1 (MQ=255)
                          aGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCgg                                                                                                                           >  2:2559809/1‑150 (MQ=255)
                          aGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCgg                                                                                                                           >  2:64243/1‑150 (MQ=255)
                          aGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCgg                                                                                                                           >  1:1280475/1‑150 (MQ=255)
                           gAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGc                                                                                                                          >  2:2488317/1‑150 (MQ=255)
                            aaaaTTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCg                                                                                                                         <  2:1413596/150‑1 (MQ=255)
                             aaaTTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGa                                                                                                                        <  1:3187057/150‑1 (MQ=255)
                                tttCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCa                                                                                                                     <  1:301164/150‑1 (MQ=255)
                                     aaaaGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGggtg                                                                                                                <  1:1516090/150‑1 (MQ=255)
                                     aaaaGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGggtg                                                                                                                <  1:2926035/150‑1 (MQ=255)
                                      aaaGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGggtgg                                                                                                               >  2:4567041/1‑150 (MQ=255)
                                      aaaGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCTGATTTTATCAACACAATCCCGGCGACCAGggtgg                                                                                                               <  2:3148213/150‑1 (MQ=255)
                                        aGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGtt                                                                                                             >  2:3449755/1‑150 (MQ=255)
                                        aGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGtt                                                                                                             >  2:1206926/1‑150 (MQ=255)
                                         gACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTa                                                                                                            >  2:545825/1‑150 (MQ=255)
                                            gTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAAcc                                                                                                         >  1:3864306/1‑150 (MQ=255)
                                             ttAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAAccc                                                                                                        <  2:1634947/150‑1 (MQ=255)
                                                 cAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCa                                                                                                    >  2:4535296/1‑150 (MQ=255)
                                                 cAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCa                                                                                                    >  1:2266494/1‑150 (MQ=255)
                                                 cAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCa                                                                                                    >  1:2266370/1‑150 (MQ=255)
                                                  aaCGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCaa                                                                                                   >  1:490534/1‑150 (MQ=255)
                                                     gATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAAttt                                                                                                <  1:40270/150‑1 (MQ=255)
                                                        ttCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTtca                                                                                             <  2:4485192/150‑1 (MQ=255)
                                                         tCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTtcat                                                                                            <  1:507660/150‑1 (MQ=255)
                                                         tCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTtcat                                                                                            >  1:3681833/1‑150 (MQ=255)
                                                           cAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTtcatca                                                                                          <  2:3243342/150‑1 (MQ=255)
                                                           cAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTtcatca                                                                                          >  1:1701564/1‑150 (MQ=255)
                                                           cAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTtcatca                                                                                          >  1:40615/1‑150 (MQ=255)
                                                            aGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAg                                                                                         >  1:2907129/1‑150 (MQ=255)
                                                              cACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCg                                                                                       >  2:1506455/1‑150 (MQ=255)
                                                               aCGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGa                                                                                      >  1:4171653/1‑150 (MQ=255)
                                                               aCGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGa                                                                                      <  1:624258/150‑1 (MQ=255)
                                                                  aTTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATaa                                                                                   >  2:1789416/1‑150 (MQ=255)
                                                                    tGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAAc                                                                                 >  2:398650/1‑150 (MQ=255)
                                                                     gCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACt                                                                                <  1:1368267/150‑1 (MQ=255)
                                                                     gCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGGTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACt                                                                                <  2:2289997/150‑1 (MQ=255)
                                                                       cAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACttt                                                                              <  1:1651903/150‑1 (MQ=255)
                                                                       cAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACttt                                                                              <  2:1885385/150‑1 (MQ=255)
                                                                          tGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTtcgt                                                                           <  1:1152364/150‑1 (MQ=255)
                                                                             cagccagGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGa                                                                        >  1:306335/1‑150 (MQ=255)
                                                                               gccagGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAAt                                                                      <  1:2026050/150‑1 (MQ=255)
                                                                                ccagGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATa                                                                     <  2:1542533/150‑1 (MQ=255)
                                                                                 cagGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATaa                                                                    >  1:4119381/1‑150 (MQ=255)
                                                                                 cagGCGGCGTGAACCCATTCAAACTGAGCCATGGTTCTCCTCGAGTTCAGGTTTTCGAGCTTTACGGGTTCATGATTTTTTCAACCCTTTCCCGGCGACCAGGGTGGTTATCCCGGCAATTTTCTTCAGCTATATACTTTCTTCGTATaa                                                                    <  2:828958/150‑1 (MQ=255)
                                                                                   ggcggcGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACa                                                                  <  1:3896465/150‑1 (MQ=255)
                                                                                    gcggcgTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACaa                                                                 <  1:1871660/150‑1 (MQ=255)
                                                                                     cggcgTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACaaa                                                                >  2:3568017/1‑150 (MQ=255)
                                                                                      ggcgTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACaaaa                                                               >  2:483515/1‑150 (MQ=255)
                                                                                      ggcgTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACaaaa                                                               >  2:4168153/1‑150 (MQ=255)
                                                                                          tGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCt                                                           <  2:3958882/150‑1 (MQ=255)
                                                                                          tGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCt                                                           >  2:1925408/1‑150 (MQ=255)
                                                                                          tGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCt                                                           >  2:4074734/1‑150 (MQ=255)
                                                                                            aaCCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTaa                                                         <  2:1051386/150‑1 (MQ=255)
                                                                                                   tCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGt                                                  <  2:2800214/150‑1 (MQ=255)
                                                                                                    cAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGta                                                 >  1:3992674/1‑150 (MQ=255)
                                                                                                    cAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGta                                                 <  1:3492728/150‑1 (MQ=255)
                                                                                                    cAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGta                                                 >  1:1203733/1‑150 (MQ=255)
                                                                                                         tGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTaataaa                                            <  2:1575072/150‑1 (MQ=255)
                                                                                                              cATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATaaataaaa                                       <  2:2619335/150‑1 (MQ=255)
                                                                                                              cATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATaaataaaa                                       >  1:4339962/1‑150 (MQ=255)
                                                                                                               aTGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATaaataaaat                                      >  1:1023678/1‑150 (MQ=255)
                                                                                                                 ggTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATAc                                    >  1:2400606/1‑150 (MQ=255)
                                                                                                                  gTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATAcc                                   <  1:3310249/150‑1 (MQ=255)
                                                                                                                   ttCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCg                                  <  1:2636602/150‑1 (MQ=255)
                                                                                                                      aCCTCCAGTTCAGGTTTACGCGCTTTACGGGGACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATa                               <  1:2514864/150‑1 (MQ=255)
                                                                                                                        cTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATAcc                             >  1:3548050/1‑150 (MQ=255)
                                                                                                                        cTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATAcc                             >  2:3872044/1‑150 (MQ=255)
                                                                                                                             gTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTccc                        <  2:1937818/150‑1 (MQ=255)
                                                                                                                             gTTCAGGTTTACGCGCTTTACGGGTACCTGATTNTATCAACACNATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTccc                        <  1:2577363/150‑1 (MQ=255)
                                                                                                                              ttCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCca                       <  1:4616786/150‑1 (MQ=255)
                                                                                                                               tCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCcac                      <  2:3633497/150‑1 (MQ=255)
                                                                                                                                 aGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAg                    >  2:2338603/1‑150 (MQ=255)
                                                                                                                                  ggTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAgc                   >  2:573972/1‑150 (MQ=255)
                                                                                                                                     ttACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCa                >  2:2147265/1‑150 (MQ=255)
                                                                                                                                      tACGCGCTTTACGGGGACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCAt               <  1:3796908/150‑1 (MQ=255)
                                                                                                                                        cgcgCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATaa             <  2:1095328/150‑1 (MQ=255)
                                                                                                                                        cgcgCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATaa             <  2:3645973/150‑1 (MQ=255)
                                                                                                                                         gcgcTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAg            <  2:2952628/150‑1 (MQ=255)
                                                                                                                                            cTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTa         <  2:2751806/150‑1 (MQ=255)
                                                                                                                                            cTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTa         >  1:1053948/1‑150 (MQ=255)
                                                                                                                                             tttACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTAc        <  2:4639418/150‑1 (MQ=255)
                                                                                                                                             tttACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTAc        <  2:807227/150‑1 (MQ=255)
                                                                                                                                                  gggTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTa   <  2:153470/150‑1 (MQ=255)
                                                                                                                                                   ggTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTaa  >  2:1130877/1‑150 (MQ=255)
                                                                                                                                                   ggTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTaa  >  1:2553461/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCAAATATTTTGCGACGAAAGCCGTCAGAAAATTTCAAAAAGACGTTAGCAACGATTTCCAGCACGATTGCCAATGCCAGCCAGGCGGCGTGAACCCATTCAAACTGCGCCATGGTTCACCTCCAGTTCAGGTTTACGCGCTTTACGGGTACCTGATTTTATCAACACAATCCCGGCGACCAGGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAA  >  CP009273/1667294‑1667590

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