Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099741,896,1840AG70.6% 2348.1 / inf 3354intergenic (‑3962/‑4635)RPPX_08470/RPPX_08480tryptophan synthase subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (419/359);  new base G (910/956);  total (1701/1653)
Fisher's exact test for biased strand distribution p-value = 1.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 3379 total aligned reads displayed.

GGCCACCACTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGAGCCAGCGGCTGGTTGACGTGCACGGTGAAGGTCGGTTTGACGTTCTCGGCTACCGAAGTCTGGTCGGCC  >  CP009974/1896034‑1896334
                                                                                                                                                      |                                                                                                                                                       
 gCCACGACTTCGTCGAGGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCGGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCAt                                                                                                                                                        <  1:766720/151‑1 (MQ=38)
    accacTTCGTCGGTGGTGTCGGTGACTTGAACGGTAGCTGCACCGCCCAACTCAAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGATAGACGTCATCGCCTTGCGCATCGt                                                                                                                                                     >  1:3728906/1‑151 (MQ=255)
    accacTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCGTCCACTGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCCTGTGcgtcgt                                                                                                                                                     <  1:2761241/151‑1 (MQ=9)
      tacTTCGTCGGTGGTGTCGGTCACTTGGACCGAAGCCGCACCGCCCAACTCAAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCATCGTgg                                                                                                                                                   <  1:68345/150‑1 (MQ=37)
      gacTTCGTCGAGGGTGTCGGTGACTTGAACCGAAGCGTTGCCGCCCAATTCCAGGTTCTCGAACGCAGCGCCGGTGGCGTCTACTGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTgg                                                                                                                                                   <  1:2947780/150‑1 (MQ=12)
      gacTTCGTCGAGGGTGTCGGTGACTTGAACCGAAGCGTTGCCGCCCAATTCCAGGTTCTCGAACGCAGCGCCGGTGGCGTCTACTGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTgg                                                                                                                                                   >  1:4448767/2‑151 (MQ=12)
      cacTTCGTCGGTGGTGTCGGTTACTTGGACCGAAGCAGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTgg                                                                                                                                                   >  2:1292299/1‑151 (MQ=25)
        gttCGTCGGGGGTGTCGGTCACTTTAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGATAGACGTCATCGCCTTGCGCGTCATGGgt                                                                                                                                                 <  2:150235/150‑1 (MQ=255)
        cTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCGTCTACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGgt                                                                                                                                                 <  1:1724542/151‑1 (MQ=16)
         ttCGTCGGTGGTGTCGGTCACTTGGACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGgtg                                                                                                                                                <  1:4530741/151‑1 (MQ=255)
           cgtcgGTGGTGTCGGTTACTTGGACCGAAGCAGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTa                                                                                                                                              <  1:4722712/151‑1 (MQ=25)
             tcgagagtgtCGGTGACTTGGACCGAAGCGTTGCCGCCCAATTCCAGGTTCTCGAACGCAGCGCCGGTGGCGTCTACTGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACg                                                                                                                                            <  2:2658663/145‑1 (MQ=12)
              cgGTGGTGTCGGTCACTTGGACCGAAGCCGCACCGCCCAACTCAAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACgg                                                                                                                                           >  2:1607214/1‑151 (MQ=255)
                 gggtgTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCGGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGgcgc                                                                                                                                        <  1:3991196/150‑1 (MQ=255)
                  ggtgTCGGTGACTTGAACCGAAGCGTTGCCGCCCAATTCCAGGTTCTCGAACGCAGCGCCGGTGGCGTCTACTCCCCAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTCCGTCGTGGCAGGACGgcgcg                                                                                                                                       >  1:4823832/1‑151 (MQ=12)
                   gtgtcggtgacatccactttcgcgGCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGgcgcgc                                                                                                                                      <  1:4363641‑M1/127‑1 (MQ=255)
                    tgtcggtgacatccactttcgcgGCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCATCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGATAGACATCATCGCCTTGCGCGTCGTGGGTGTACGGCGCGCt                                                                                                                                     <  2:2986725‑M1/128‑1 (MQ=255)
                    tgtCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGATAGACGTCATCGCCTTGCGCGTCATGGGTGTACGGCGCGCt                                                                                                                                     >  1:3032235/1‑151 (MQ=38)
                     gtCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCGTCTACCGCCGAGTTGAAGCCCAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTCGCGCCGCGCGGGCGACCGGCGCCCtc                                                                                                                                    >  2:4650664/1‑150 (MQ=255)
                       cGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCGGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCACTGGt                                                                                                                                  <  2:309904/151‑1 (MQ=255)
                        ggTCACTTGAACCGAAGCCGCACCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCGTCCACTGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGtt                                                                                                                                 >  2:5064628/1‑151 (MQ=9)
                         gtgaCTTGAACGGTAGCTGCACCGCCCAACTCAAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGATAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGtct                                                                                                                                >  1:2851054/4‑149 (MQ=21)
                          tgaCTTGAACCGAAGCGTTGCCGCCCAATTCCAGGTTCTCGAACGCAGCGCCGGTGGCGTCTACTGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTc                                                                                                                               <  2:823818/149‑1 (MQ=9)
                           cACTTGAACAGTAGCCGCACCGCCCAGTTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCg                                                                                                                              >  1:887216/1‑151 (MQ=16)
                             cTTGAACCGGAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGATAGACGTCATCGCCTTGCGCGTCATGGGTGTACGGCGCGCTGGTTTCGcc                                                                                                                            <  2:4238043/151‑1 (MQ=37)
                              ttGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGATAGACGTCATCGCCTTGCGCGTCATGGGTGTACGGCGCGCTGGTTTCGcca                                                                                                                           >  2:2026597/1‑150 (MQ=37)
                               tggaCCGAAGCCGCACCGCCCAACTCAAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTTTCGCCgg                                                                                                                          <  2:1734018/148‑1 (MQ=25)
                                 aacagtaGCCGCACCGCCCAGTTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGcc                                                                                                                        <  1:1814606/145‑1 (MQ=16)
                                   ctgtcgcgGCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCtt                                                                                                                      >  2:3428322‑M1/9‑151 (MQ=255)
                                    tttcgcgGCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCATCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGATAGACATCATCGCCTTGCGCGTCGTGGGTGTACGGCGCGCTGGTTTCACCGGCCTTg                                                                                                                     >  2:1820548‑M1/8‑151 (MQ=255)
                                         gggCACCGCCCAACTCGAGGTGCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGt                                                                                                                <  2:3738740/149‑1 (MQ=17)
                                         ccGCACCGCCCAACTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGt                                                                                                                <  1:2157295/151‑1 (MQ=16)
                                          gGCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTg                                                                                                               <  2:276977‑M1/150‑1 (MQ=255)
                                          gGCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTg                                                                                                               <  1:128469‑M1/150‑1 (MQ=255)
                                           gCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGa                                                                                                              >  2:653196/1‑151 (MQ=17)
                                            tgccGCCCAATTCCAGGTTCTCGAACGCAGCGCCGGTGGCGTCTACTGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGAc                                                                                                             <  1:2331708/149‑1 (MQ=14)
                                            cACCGCCCAGTTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTCTCGCCGGCCTTGATGGTGAc                                                                                                             >  2:139681/1‑151 (MQ=18)
                                             cccGCCCAACTCAAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGAcc                                                                                                            <  2:4427995/150‑1 (MQ=21)
                                             aCCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGATAGACGTCATCGCCTTGCGCGTCATGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGAcc                                                                                                            >  1:4163495/1‑151 (MQ=21)
                                              ccGCCCAACTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCt                                                                                                           <  2:3831012/151‑1 (MQ=18)
                                               cGCCCAACTCAAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTg                                                                                                          >  2:2862000/1‑151 (MQ=18)
                                                  ccAACTCGAGGTTCTCGAACGCAGCACCAGTGGCATCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGATAGACATCATCGCCTTGCGCGTCGTGGGTGTACGGCGCGCTGGTTTCACCGGCCTTGATGGTGACCTGGGc                                                                                                       >  1:3326356/1‑151 (MQ=18)
                                                  ccAACTCAAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGATAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTCTCGCCGGCCTTGATGGTGACCTgagc                                                                                                       >  1:1614447/1‑148 (MQ=18)
                                                     aCTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCAtt                                                                                                    <  2:1919173/151‑1 (MQ=17)
                                                     aCTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCAtt                                                                                                    <  2:2579785/151‑1 (MQ=17)
                                                       tgcaGGTTCTCGAAGGTCGCACCAGTGGCGTCCACTGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCAttgt                                                                                                  >  2:903628/4‑151 (MQ=2)
                                                       tgcaGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCGttgt                                                                                                  >  1:5361274/4‑151 (MQ=9)
                                                       tCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCAttgt                                                                                                  >  2:3634518/1‑151 (MQ=17)
                                                         caGGTTCTCGAAGGTCGCACCAGTGGCGTCCACTGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCGttgttg                                                                                                >  1:3546967/2‑151 (MQ=2)
                                                          aGGTTCTCGAAGGTCGCACCAGTGGCGTCCACTGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCATTGTTGc                                                                                               >  2:1187095/1‑151 (MQ=2)
                                                          aGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGATAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTCTCGCCGGCCTTGATGGTGACCTGAGCGTTGTTGc                                                                                               >  1:976460/1‑151 (MQ=18)
                                                           ggTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGTCCTGAGCATGGTTGCt                                                                                              >  2:3943201/1‑151 (MQ=14)
                                                             ttggcGATCGCAGCACCAGTGGCATCCGCCGTGGATTTGATGGCCTGGGTGATCTGGCCCGCATGTTGATAGACATCTTCGTGTTGCGCGTCGTGGGTGTACGGCGCGCTGGTTTCACCGGCCTTGCTGGTGACCTGGGCATTGTTGCtga                                                                                            <  2:3192163/147‑3 (MQ=9)
                                                             ttCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGATAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTCTCGCCGGCCTTGATGGTGACCTGAGCGTTGTTGCTcc                                                                                            >  2:1063901/1‑150 (MQ=17)
                                                                 cGAACGCAGCACCAGTGGCATCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGATAGACATCATCGCCTTGCGCGTCGTGGGTGTACGGCGCGCTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGt                                                                                        <  1:3231797/151‑1 (MQ=17)
                                                                  gAACGCAGCGCCGGTGGCGTCTACTGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCATTGTTGCTCAGGGTc                                                                                       <  1:4078378/151‑1 (MQ=11)
                                                                  gAACGCAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTc                                                                                       <  1:4892072/151‑1 (MQ=14)
                                                                   aaCGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGATAGACGTCATCGCCTTGCGCGTCATGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCa                                                                                      <  1:5085282/151‑1 (MQ=18)
                                                                    aCGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGATAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTCTCGCCGGCCTTGATGGTGACCTGAGCGTTGTTGCTCAGGGTCAc                                                                                     <  2:556032/151‑1 (MQ=17)
                                                                     cGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCATTGTTGCTCAGGGTCACg                                                                                    <  1:2638886/151‑1 (MQ=11)
                                                                      gttgcgcCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCGTTGTTGCTCAGGGTCACGa                                                                                   <  1:343423/148‑1 (MQ=9)
                                                                      gtaGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAGGGTCACGa                                                                                   >  1:606856/3‑151 (MQ=2)
                                                                       cAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGAc                                                                                  <  1:3826150/151‑1 (MQ=14)
                                                                         gcaccAGTGGCGTCCACTGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCa                                                                                <  1:1489107/148‑1 (MQ=11)
                                                                          cgcCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAGGGTCACGACCaa                                                                               >  2:5267593/1‑150 (MQ=11)
                                                                          cgcCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAg                                                                               <  2:4025238/151‑1 (MQ=11)
                                                                           gccggTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGAGCGTTGTTGCTCAGGGTCACGACCAgg                                                                              <  1:4628245/147‑1 (MQ=14)
                                                                             cAGTGGCGTCTACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCACTGGTTTCGCCAGCCTTGATGGTGACCTGGGCGTTGTTGCTCAGGGTCACGACCAAGTc                                                                            >  1:425512/1‑151 (MQ=14)
                                                                              aGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCATTGTTGCTCAGGGTCACGATCAGGTCg                                                                           <  1:204228/151‑1 (MQ=11)
                                                                               gTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGt                                                                          <  2:4942545/151‑1 (MQ=11)
                                                                               gTGGCATCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGATAGACATCATCGCCTTGCGCGTCGTGGGTGTACGGCGCGCTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGt                                                                          <  1:3907608/151‑1 (MQ=17)
                                                                                tGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAGGGTCACGACCAAGTCGTg                                                                         >  2:4771819/1‑151 (MQ=255)
                                                                                 ggCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCGTTGTTGCTCAGGGTCACGACTAGGTCGTGa                                                                        <  1:193410/151‑1 (MQ=9)
                                                                                  gCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGAg                                                                       <  2:2477240/151‑1 (MQ=11)
                                                                                  gCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGAg                                                                       >  2:4866786/1‑151 (MQ=11)
                                                                                  gCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGATAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTCTCGCCGGCCTTGATGGTGACCTGAGCGTTGTTGCTCAGGGTCACGACCAGGTCGTggg                                                                       <  1:5256282/151‑3 (MQ=14)
                                                                                    gTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGATAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTCTCGCCGGCCTTGATGGTGACCTGAGCGTTGTTGCTCAGGGTCACGACCAGGTCGTGGGcc                                                                     <  1:1386749/151‑1 (MQ=11)
                                                                                       taCCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGGGCCAGc                                                                  <  1:1000780/150‑1 (MQ=14)
                                                                                           gccgAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAGGGTCACGACCAAGTCGTGGGCCAGCGGCt                                                              >  2:739823/1‑151 (MQ=17)
                                                                                                      gCCCAGGCTGATCTGGCCAGCATCCTGATAGACGTCATCGCCTTGCGCGTCATGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAAGTCGTGGGCCAGCGGCTGGTTGACGTGc                                                  <  1:1087058/151‑1 (MQ=11)
                                                                                                          aGGCTGATCTGGCCAGCATCCTGATAGACGTCATCGCCTTGCGCGTCATGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAAGTCGTGGGCCAGCGGCTGGTTGACGTGCACgg                                              >  1:4968575/1‑151 (MQ=11)
                                                                                                            gCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAGGGTCACGACTAGGTCGTGAGCCAGGGCGGTGTTGATGTGCACGGTg                                            <  1:5180672/151‑1 (MQ=11)
                                                                                                               gatctggCCCGCATCCTGATAGACATCATCGCCTTGCGCGTCGTGGGTGTACGGCGCGCTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGGGCCAGCGGCTGGTTGACGTGCACGGTGAAg                                         >  1:2244769/7‑151 (MQ=11)
                                                                                                                 tctggCCCGCATCGTTGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGGGCCAGCGGCTGGTTGACGTGTACGGTGAAGGt                                       <  1:1180239/147‑1 (MQ=14)
                                                                                                                  ctggCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCACTGGTTTCGCCAGCCTTGATGGTGACCTGGGCGTTGTTGCTCAGGGTCACGACCAAGTCGTGGGCCAGCGGCTGGTTGACGTGCACGGTGAAGGTc                                      <  2:4141384/148‑1 (MQ=11)
                                                                                                                    ggccagCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAGGGTCACGACTAGGTCGTGAGCCAGGGCGGTGTTGATGTGCACGGTGAACGTCgg                                    <  1:3628421/146‑1 (MQ=11)
                                                                                                                    ggCCCGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGAGCGTTGTTGCTCAGGGTCACGACCAGGTCGTGAGCCAGGGCGGTGTTGATGTGCACGGTGAACGTCgg                                    <  2:4511142/150‑1 (MQ=11)
                                                                                                                     gCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTCGATGGTGACCTGGGCGTTGTTGCTCAGGGTCACGACCAAGTCGTGGGCCAGCGGCTGGTTGACGTGCACGGTGAAGGGCGGt                                   >  2:2402353/1‑151 (MQ=11)
                                                                                                                     gCCCGCATCCTGATAGACATCATCGCCTTGCGCGTCGTGGGTGTACGGCGCGCTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGGGCCAGCGGCTGGTTGACGTGCACGGTGAAGGTCGGt                                   <  2:2313819/151‑1 (MQ=11)
                                                                                                                       cagCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCGTTGTTGCTCAGGGTCACGACTAGGTCGTGAGCCAGGGCGGTGTTGATGTGCACGGTGAACGTCGGttt                                 >  1:2448583/3‑151 (MQ=11)
                                                                                                                        agCATCCTGATAGACGTCATCGCCTTGCGCGTCATGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAAGTCGTGGGCCAGCGGCTGGTTGACGTGCACGGTGAAGGTCGGTTTg                                >  2:5171666/2‑151 (MQ=11)
                                                                                                                              cTGATAGACGTCATCGCCTTGCGCGTCATGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAAGTCGTGGGCCAGCGGCTGGTTGACGTGCACGGTGAAGGTCGGTTTGACGTtc                          >  2:5381831/1‑151 (MQ=11)
                                                                                                                               ttgtgggcGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCAGGCCTTGATGGTGACCTGAGCGTTGTTGCTCAGGGTCACGACCAGGTCGTGAGCCAGGGCGGTGTTGATGTGCACGGTGAACGTCGGTTTGCCGTtct                         <  2:211499/144‑1 (MQ=11)
                                                                                                                                tgTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGGGCCAGCGGCTGGTTGACGTGTACGGTGAAGGTCGGTTTGACGTtctc                        >  2:3318746/2‑151 (MQ=255)
                                                                                                                                gatAGACGTCAGCGCCTTGCGCGTCGTGGGTGTACGGCGCGCTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGGGCCAGCGGCTGGTTGACGTGCACGGTGAAGGTCGGTTTGACGTtctc                        <  2:3539799/149‑1 (MQ=11)
                                                                                                                                       gTCATCGCCTTGCGCGTCGTGGGTGTACGGCGCGCTGGTCTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGGGCCAGCGGCTGGTTGACGTGCACAGTGAAGGTTGGTTTGACGTTCTCGGCTAcc                 <  1:5436406/151‑1 (MQ=11)
                                                                                                                                          ttCGCCTTGCGCGTCGTGGGTGTACGGCGCGCTGGTCTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGGGCCAGCGGCTGGTTGACGTGCACAGTGAAGGTTGGTTTGACGTTCTCGGCTACCGaa              <  1:3134876/150‑1 (MQ=11)
                                                                                                                                            cGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCGTTGTTGCTCAGGGTCACGACTAGGTCGTGAGCCAGGGCGGTGTTGATGTGCACGGTGAACGTCGGTTTGACGTTCTCGGCCACCGAAGt            <  2:2215189/151‑1 (MQ=11)
                                                                                                                                              ccTTGCGCATCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGGGCCAGCGGCTGGTTGACGTGTACGGTGAAGGTCGGTTTGACGTTCTCGGCTACCGAAGTCt          <  2:4323375/151‑1 (MQ=255)
                                                                                                                                                   cgcATCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAGGGTCACGACCAAGTCGTGGGCCAGCGGCTGGTTGACGTGCACGGTGAAGGTCGGTTTGACGTTCTCGGCTACCGAGGTCTGGTCg     >  1:2535669/1‑151 (MQ=1)
                                                                                                                                                      |                                                                                                                                                       
GGCCACCACTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGAGCCAGCGGCTGGTTGACGTGCACGGTGAAGGTCGGTTTGACGTTCTCGGCTACCGAAGTCTGGTCGGCC  >  CP009974/1896034‑1896334

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