Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,875,671 G→A 28.0% T4I (ACC→ATC)  dnaN ← DNA polymerase III, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,875,6710GA28.0% 694.1 / 303.9 386T4I (ACC→ATC) dnaNDNA polymerase III, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (145/132);  new base A (45/63);  total (190/196)
Fisher's exact test for biased strand distribution p-value = 6.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

Only 100 of 390 total aligned reads displayed.

CCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCT  >  CP009273/3875522‑3875819
                                                                                                                                                     |                                                                                                                                                    
  aTTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTa                                                                                                                                                    >  2:2847412/1‑150 (MQ=255)
     tCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAAt                                                                                                                                                 <  1:4946367/150‑1 (MQ=255)
     tCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAAt                                                                                                                                                 <  2:1637020/150‑1 (MQ=255)
     tCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAAt                                                                                                                                                 <  1:2825250/150‑1 (MQ=255)
     tCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAAt                                                                                                                                                 <  2:2360482/150‑1 (MQ=255)
       cATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAAttt                                                                                                                                               <  1:4465523/150‑1 (MQ=255)
         tctcGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCa                                                                                                                                             <  1:805859/150‑1 (MQ=255)
             gagaTCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAgg                                                                                                                                         >  2:2324142/1‑150 (MQ=255)
              agaTCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGt                                                                                                                                        <  1:5108541/150‑1 (MQ=255)
              agaTCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGt                                                                                                                                        >  1:5307836/1‑150 (MQ=255)
               gaTCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGtt                                                                                                                                       <  1:574276/150‑1 (MQ=255)
                aTCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGttt                                                                                                                                      <  1:2401232/150‑1 (MQ=255)
                aTCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGttt                                                                                                                                      >  1:3792840/1‑150 (MQ=255)
                  cAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTAc                                                                                                                                    <  1:1050639/150‑1 (MQ=255)
                   aGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACg                                                                                                                                   >  2:2084669/1‑150 (MQ=255)
                   aGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAAGCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCAGAGGTTTACg                                                                                                                                   >  2:3964442/1‑150 (MQ=255)
                     tACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGAt                                                                                                                                 <  1:5218876/150‑1 (MQ=255)
                     tACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGAt                                                                                                                                 <  1:5261758/150‑1 (MQ=255)
                     tACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCAGAGGTTTACGAt                                                                                                                                 >  1:1880926/1‑150 (MQ=255)
                        cGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGAc                                                                                                                              <  2:64024/150‑1 (MQ=255)
                         ggTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACa                                                                                                                             <  2:4227111/150‑1 (MQ=255)
                            cAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATg                                                                                                                          <  2:1872803/150‑1 (MQ=255)
                             aGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGt                                                                                                                         >  2:324382/1‑150 (MQ=255)
                                gaaaaCGTACCGTCAGCAACCAGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACACTGTTCt                                                                                                                      <  2:1269966/147‑1 (MQ=255)
                                gACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCt                                                                                                                      <  1:115344/150‑1 (MQ=255)
                                gACAACGTACCGTCAGCAACCTGTAACAGCAGATTACAGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTATCGGTTTTAATAAATGATCACTTTCTACGGTAAATTTCATATGTTTACGATGACAATGTTCt                                                                                                                      >  2:4950305/1‑150 (MQ=255)
                                 aCAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTg                                                                                                                     <  2:33425/150‑1 (MQ=255)
                                        aCCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCTCCTAACGGACCTCTCTCCTGTTGTAGCTGTTTTAATAAATGCTCACGTTCTACGGTAAATTACATAGGTTTACGATGACAATGTTCTGATTAAAt                                                                                                              <  1:2225498/150‑1 (MQ=255)
                                        aCCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTGTTAATAAATGCTCACGTTCTACGATAAATATCATAGGTTTACGATGACAAGGTTCTGATAAAAt                                                                                                              >  1:3587286/1‑150 (MQ=255)
                                          cGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAAttt                                                                                                            <  1:3258833/150‑1 (MQ=255)
                                           gTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTg                                                                                                           >  2:4684253/1‑150 (MQ=255)
                                            tCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGa                                                                                                          >  1:988948/1‑150 (MQ=255)
                                              aGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGaaa                                                                                                        >  1:3281592/1‑150 (MQ=255)
                                              aGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGaaa                                                                                                        <  2:3577694/150‑1 (MQ=255)
                                                       tAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATcttctt                                                                                               <  2:3673960/150‑1 (MQ=255)
                                                       tAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATcttctt                                                                                               >  1:700859/1‑150 (MQ=255)
                                                         aCAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTg                                                                                             <  1:2679869/150‑1 (MQ=255)
                                                         aCAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTg                                                                                             <  2:4754588/150‑1 (MQ=255)
                                                         aCAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTg                                                                                             <  1:4691914/150‑1 (MQ=255)
                                                              agATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATc                                                                                        <  2:4057045/150‑1 (MQ=255)
                                                              agATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATc                                                                                        <  1:1114781/150‑1 (MQ=255)
                                                               gATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCg                                                                                       >  2:1478548/1‑150 (MQ=255)
                                                               gATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCCTCTTTGATATCg                                                                                       <  2:3031038/150‑1 (MQ=255)
                                                                   aCCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCATTGATATCGTGGc                                                                                   <  2:556909/150‑1 (MQ=255)
                                                                      gagaATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGctct                                                                                >  2:2898599/1‑150 (MQ=255)
                                                                         aaTCGGTAGCGTAGGACGACCACCTAACGGGCCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCtt                                                                             <  1:415386/150‑1 (MQ=255)
                                                                         aaTCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCtt                                                                             >  1:5059552/1‑150 (MQ=255)
                                                                         aaTCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCtt                                                                             >  2:755854/1‑150 (MQ=255)
                                                                          ctCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCGGTTGTAGCGGTTTTAATAAATGCACACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTc                                                                            <  2:5026198/149‑1 (MQ=255)
                                                                          aTCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTc                                                                            <  2:1236421/150‑1 (MQ=255)
                                                                          aTCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTc                                                                            >  1:1544050/1‑150 (MQ=255)
                                                                             ggTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCTCg                                                                         >  2:893307/1‑150 (MQ=255)
                                                                                aGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCaa                                                                      <  1:5155243/150‑1 (MQ=255)
                                                                                aGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCaa                                                                      >  2:2584825/1‑150 (MQ=255)
                                                                                aGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCaa                                                                      >  1:1001516/1‑150 (MQ=255)
                                                                                   gTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATTGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTg                                                                   <  2:2455985/150‑1 (MQ=255)
                                                                                      ggACGACCACCTAACGGACCGCTCACCTCTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTAAAAGGTTTAAGATGACAATTTTCTGATTAAATTTGAAAAATCTTCTATTATATCGTGGCTCTCTTCACGCAACTGATc                                                                >  1:3208040/1‑150 (MQ=255)
                                                                                         cgacCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGAt                                                             <  1:310817/150‑1 (MQ=255)
                                                                                           accaccTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCt                                                           <  1:3675385/150‑1 (MQ=255)
                                                                                           accaccTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCt                                                           <  1:2089495/150‑1 (MQ=255)
                                                                                            ccaccTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCtt                                                          >  1:1441907/1‑150 (MQ=255)
                                                                                                  aaCGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACggca                                                    <  2:2150807/150‑1 (MQ=255)
                                                                                                    cGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACggcagg                                                  <  1:3041432/150‑1 (MQ=255)
                                                                                                      gACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACggcaggca                                                >  1:4666839/1‑150 (MQ=255)
                                                                                                       aCCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCAt                                               <  1:3759282/150‑1 (MQ=255)
                                                                                                       aCCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCAt                                               >  2:4813540/1‑150 (MQ=255)
                                                                                                       aCCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCAt                                               >  2:4303686/1‑150 (MQ=255)
                                                                                                        ccGCTCAACTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATg                                              <  2:3097211/150‑1 (MQ=255)
                                                                                                          gCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGaa                                            <  1:4001869/150‑1 (MQ=255)
                                                                                                           cTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAg                                           <  2:5077872/150‑1 (MQ=255)
                                                                                                           cTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAg                                           >  2:1728816/1‑150 (MQ=255)
                                                                                                           cTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAg                                           >  1:1198865/1‑150 (MQ=255)
                                                                                                             cACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCa                                         >  2:1386270/1‑150 (MQ=255)
                                                                                                              aCCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCAc                                        <  1:2029072/150‑1 (MQ=255)
                                                                                                               ccTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCAcc                                       <  1:4517404/150‑1 (MQ=255)
                                                                                                               ccTGTTGTAGCGGATTTAATAACTGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCAcc                                       <  1:3917592/150‑1 (MQ=255)
                                                                                                                cTGTTGTAGCTGTTTTAATAAATGCTCACGTTCTAAGATAAATTTCATAGGTTTACGATGACATTGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTATGGCAGGCATGAAGCACcg                                      <  1:2272909/150‑1 (MQ=255)
                                                                                                                 tgttgtAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACcgt                                     <  2:3804043/150‑1 (MQ=255)
                                                                                                                  gttgtAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACcgtc                                    >  2:1130755/1‑150 (MQ=255)
                                                                                                                    tgtAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACcgtcgt                                  >  1:4386464/1‑150 (MQ=255)
                                                                                                                     gtAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCgtg                                 >  1:1691098/1‑150 (MQ=255)
                                                                                                                     gtAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCgtg                                 <  1:5176056/150‑1 (MQ=255)
                                                                                                                      tAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCgtgt                                >  1:4087844/1‑150 (MQ=255)
                                                                                                                           gTTTTAATAAATGCTCACGTTCTACGATAAATTTCATAGGTTTACTATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCACGCATGAAGCACAGTCGTGTGGTCa                           >  2:94888/1‑150 (MQ=255)
                                                                                                                             tttAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACg                         >  1:74924/1‑150 (MQ=255)
                                                                                                                               taataaATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGc                       <  1:4249270/150‑1 (MQ=255)
                                                                                                                                aataaATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGcc                      >  2:2872360/1‑150 (MQ=255)
                                                                                                                                aataaATGCTCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGcc                      >  1:291315/1‑150 (MQ=255)
                                                                                                                                 ataaATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGcca                     >  2:2726047/1‑150 (MQ=255)
                                                                                                                                        cTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACg              >  1:4603198/1‑150 (MQ=255)
                                                                                                                                         tCACGTTCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGc             <  2:3480891/150‑1 (MQ=255)
                                                                                                                                            cGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGTCATGAAGCACCGTCGTGTGGTCACGGCCACCAAAAGAATc          >  2:4473956/1‑150 (MQ=255)
                                                                                                                                            cGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATc          >  1:1373535/1‑150 (MQ=255)
                                                                                                                                            cGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATc          >  1:2502534/1‑150 (MQ=255)
                                                                                                                                             gTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCg         >  1:1304051/1‑150 (MQ=255)
                                                                                                                                             gTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCg         >  2:345547/1‑150 (MQ=255)
                                                                                                                                               tCTACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCACACGCATCGcc       >  1:4862995/1‑150 (MQ=255)
                                                                                                                                                 tACGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCaa     >  2:5162851/1‑150 (MQ=255)
                                                                                                                                                  aCGATAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAAt    <  2:4605482/150‑1 (MQ=255)
                                                                                                                                                   cGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAAtc   <  2:869160/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCATTTCCATCTCGAGATCAGTACCGGTCAGCGACAACGTACCGTCAGCAACCTGTAACAGCAGATTACCGAGAATCGGTAGCGTAGGACGACCACCTAACGGACCGCTCACCTGTTGTAGCGGTTTTAATAAATGCTCACGTTCTACGGTAAATTTCATAGGTTTACGATGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCT  >  CP009273/3875522‑3875819

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