Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099741,896,1210CT67.8% 2272.5 / inf 4004intergenic (‑3899/‑4698)RPPX_08470/RPPX_08480tryptophan synthase subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (502/511);  new base T (1076/1061);  total (2034/1970)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 4018 total aligned reads displayed.

GGTCAGCGTGATGGTGTAAGTGATCTCGCCGCCTTCGGTGACCGAAGGGGTCGCGGTCAGCTTGGCCACCACTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGAGCCAGC  >  CP009974/1895971‑1896271
                                                                                                                                                      |                                                                                                                                                      
    aGCGTGATGGTGTAGGTGATCTCGCCGCCTTCGGTCACCGACGGGGTGGCAGTCAGCTTGGCCACGACTTCATCGAGGGTGTCGGTGACATCCACTTTCGCGGCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAAccg                                                                                                                                                    >  1:5829431/1‑151 (MQ=2)
       gTGATGGTGTAAGTGATCTCGCCGCCTTCGGTGACCGACGGGGTCGCGGTCAGCTTGGCCACTACTTCGTCGGTGGTGTCGGTCACTTGGACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTAccgccg                                                                                                                                                 >  1:3798527/1‑151 (MQ=17)
           tGGTGTAAGTGATCTCGCCGCCTTCGGTGACCGAAGGGGTCGCGGTCAGCTTGGCCACTACTTCGTCGGTGGTGTCGGTCACTTGGACCGAAGCCGCACCGCCCAACTCAAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAAtt                                                                                                                                             <  1:4569781/151‑1 (MQ=21)
            ggTGTAGGTGGTCTCGCCGCCGTCGGTCACCGAAGGGGGCGCGGTCAGCTTCGCCACGACTTCGTCGAGGGTGTCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAACGCAGCACCAGTGGCGTCTACTGCCGAGTTg                                                                                                                                            <  2:4177363/151‑1 (MQ=9)
             gtgtAAGTGATCTCGCCGCCTTCGGTGACCGACGGGGTCGCGGTCAGCTTGGCCACGACTTCGTCGGTGGTGTCGGTCACTTGAACAGTAGCCGCACCGCCCAGTTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGa                                                                                                                                           >  2:2214523/1‑151 (MQ=17)
              tgtAAGTGATCTCGCCGCCTTCGGTGACCGACGGGGTCGCGGTCAGCTTGGCCACGACTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCGTCTACCGCCGAGTTGAt                                                                                                                                          >  2:3400980/1‑151 (MQ=14)
                taggTGATCTCGCCGCCTTCGGTCACCGAAGGGGTCGCGGTCAGCTTCGCCACGACTTCGTCGAGGGTGTCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAACGCAGCACCAGTGGCGTCTACTGCCGAGTTGATGc                                                                                                                                        >  1:6238310/4‑151 (MQ=17)
                 aaGTGATCTCGCCGCCTTCGGTGACCGAAGGGGTCGCGGTCAGCTTGGCCACCACTTCGTCGGTGGTGTCGGTCACTTGAACAGTAGCCGCACCGCCCAGTTGCAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATAcc                                                                                                                                       <  2:1030563/151‑1 (MQ=18)
                   gTGATCTCGCCGCCTTCGGTGACCGAAGGGGTCGCGGTCAGCTTGGCCACCACTTCGTCGGTGGTGTCGGTCACTTGAACAGTAGCCGCACCGCCCAGTTGCAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGAtaccaa                                                                                                                                     >  1:634371/1‑146 (MQ=18)
                      aTCTCGCCGCCTTCGGTGACCGACGGGGTCGCGGTCAGCTTGGCCACGACTTCGTCGGTGGTGTCGGTCACTTGAACAGTAGCCGCACCGCCCAGTTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGc                                                                                                                                  >  1:155718/1‑151 (MQ=17)
                       gcgagccgccGGCGGTCACCGACGGGGGGGCAGGCAGCGTGGCCACGACTTCATCGAGGGTGTCGGTGACAGCCACGGTCGCGGCACCGCCCAACTCGAGGTTCGCGAACGCAGCGCCAGTGGGGTCAACCGCCGAGTTGATGCGCAGACt                                                                                                                                 <  2:3931401/147‑1 (MQ=255)
                        ctcGCCGCCTTCGGTGACGGACGGGGTGGCAGTCAGCTTGGCCACCACTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCGTCCACTGCCGAATTGATACCCAGGCTg                                                                                                                                >  2:6084509/1‑151 (MQ=14)
                        ctcGCCGCCTTCGGTGACCGAAGGGGGGCCGGTCAGCGTGGCCACGAGGTCGTCGCGGGTGTCGGTGACATCCAATTTCGCGGCACCGCCCACCTCGAGGTTCTCGAACGCAGCACCAGTGGCATCCATCGCCGCGGTGATGCCCACGCTg                                                                                                                                >  2:4445254/1‑151 (MQ=255)
                        ctcGCCGCATTGGGTGACAGACAAGGTGGCGATCAGCTTGGCCACGACTTCATCGAGGGTGTCGGTGACATCCACTTTCGCGTCACCGCCCATCTCGAGGTTCTCGAACGCAGCACCAGTGGCATCCACCGCCGAGTTGATGCCCAGGCTg                                                                                                                                <  1:1002190/151‑1 (MQ=16)
                          caccgccTTCGGTGACCGACGGGGTGGCAGTCAGCTTGGCCACGACTTCATCGAGGGTGTCGGTGACATGCACTTTCGCGGCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGCGTTGATGCCCAGACTGAt                                                                                                                              >  1:5697684/3‑151 (MQ=16)
                           gccgccTTCGGTCACCGAAGGGGTGGCGGTCAGCTTGGCCACGACTTCGTCGAGGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGAtc                                                                                                                             >  2:1135751/1‑150 (MQ=25)
                              gccTTCGGTGACCGAAGGGGTCGCGGTCAGCTTGGCCACTACTTCGTCGGTGGTGTCGGTCACTTGAACAGTAGCCGCACCGCCCAACTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGAtctgg                                                                                                                          >  1:5575375/1‑147 (MQ=17)
                                cTTCGGTGACCGAAGGGGTGGCGGTCAGCTCGGCCGCGACTTCGTCGAGGGTGTCGGTCACTTGAACCGAGGCCGCACCGCCCACCTCGAGGTTCTCGAACGCAGCGCCGGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGcc                                                                                                                        <  2:6206222/151‑1 (MQ=25)
                                 ttCGGTCACCGAAGGGGTCGCGGTCAGCTTCGCCACGACTTCGTCGAGGGTGTCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAACGCAGCACCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGAtctggcca                                                                                                                       >  1:2494529/1‑144 (MQ=9)
                                  tCGGTCACCGACGGGGTGGCAGTCAGCTTGGCCACGACTTCATCGAGGGTGTCGGTGACATCCACTTTCGCGGCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGACGAGTTGATGCCCAGACTGATCTGGCCCg                                                                                                                      <  2:1504688/151‑1 (MQ=12)
                                   cGGTCACCGACGGGGTGGCAGTCAGCTTGGCCACGACTTCATCGAGGGTGTCGGTGACATCCACTTTCGCGGCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGc                                                                                                                     >  1:3670421/1‑151 (MQ=9)
                                    ggTGACGGACGGGGTGGCAGTCAGCTTGGCCACCACTTCGTCGGTGGTGTCGGTGACTTGAACGGTAGCTGCACCGCCCAACTCAAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCa                                                                                                                    >  2:2549589/1‑151 (MQ=21)
                                    ggTGACCGAAGGGGTCGCGGTCAGCTTGGCCACCACTTCGTCGGTGGCGTCGGTCACTTGGACAGTAGCCGCACCGCCCAGTTGCAGGTTCTCGAACGTAGCCCCGGTGGCGTCGACCGCCGAATTGATACCAAAGCTGATCTGGCACGCa                                                                                                                    >  2:4799034/1‑151 (MQ=12)
                                      tGACGGACGGGGTGGCAGTCAGCTTGGCCACCACTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCGTCCACTGCCGAATTGATACCAAGGCTGATCTGGCCAGCATc                                                                                                                  <  1:1428577/151‑1 (MQ=9)
                                      tGACCGACGGGGTCGCGGTCAGCTTGTCCACGACTTCGTCGGTGGTGTCGGTCACTTGAACAGTAGCCGCACCGCCCAGTTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCATGCTGATCTGGCCCGCATc                                                                                                                  >  1:2983230/1‑151 (MQ=17)
                                         ccGACGGGGTCGCAGTCAGCTTGGCCACGACTTCGTCGAGGGTGTCGGTGACTTGAACCGAAGCGTTGCCGCCCAATTCCAGGTTCTCGAACGCAGCGCCGGTGGCGTCTACTGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTg                                                                                                               >  1:5363107/1‑151 (MQ=12)
                                         ccGAAGGGGTCGCGGTCAGCTTGGCCACTACTTCGTCGGTGGTGTCGGTCACTTGAACAGTAGCCGCACCGCCCAACTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATcgtt                                                                                                               >  2:5173933/1‑148 (MQ=12)
                                           gAAGGGGTCGCGGTCAGCTTGGCCACTACTTCGTCGGTGGTGTCGGTCACTTGGACCGAAGCCGCACCGCCCAACTCAAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATcgttgt                                                                                                             >  2:4660726/1‑146 (MQ=21)
                                             cggggTCGCGGTCAGCTTGGCCACGACTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCGTCTACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCGTTGTAg                                                                                                           <  2:1724905/150‑1 (MQ=9)
                                                 gtggcgGTCAGCTTGGCCACGACTTCGTCGAGGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGATAGACGt                                                                                                       >  1:2387846/4‑151 (MQ=25)
                                                 gTCGCGGTCAGCTTGGCCACTACTTCGTCGGTGGTGTCGGTCACTTGAACAGTAGCCGCACCGCCCAACTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGt                                                                                                       <  1:74378/151‑1 (MQ=16)
                                                 gTCGCGGTCAGCTTGGCCACGACTTCGTCGGTGGTGTCGGTCACTTGAACAGTAGCCGCACCGCCCAGTTGCAGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGt                                                                                                       >  2:208329/1‑151 (MQ=12)
                                                  tCGCGGTCAGCTTGGCCACTACTTCGTCGGTGGTGTCGGTCACTTGGACCGAAGCCGCACCGCCCAACTCAAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTc                                                                                                      <  1:6038275/151‑1 (MQ=21)
                                                    gcgGTCAGCTTGGCCACTACTTCGTCGGTGGTGTCGGTCACTTGGACCGAAGCCGCACCGCCCAACTCAAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCAt                                                                                                    >  1:5929579/1‑151 (MQ=21)
                                                    gcagTCAGCTTGGCCACGACTTCGTCGAGGGTGTCGGTGACTTGAACCGAAGCGTTGCCGCCCAATTCCAGGTTCTCGAACGCAGCGCCGGTGGCGTCTACTGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCAt                                                                                                    <  1:2051458/148‑1 (MQ=16)
                                                    gcagTCAGCTTGGCCACGACTTCATCGAGGGTGTCGGTGACATCCACTTTCGCGGCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCAt                                                                                                    <  1:4696222/148‑1 (MQ=255)
                                                     cgGTCAGCTTGGCCACGACTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCATCTACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATc                                                                                                   >  2:3011526/1‑151 (MQ=12)
                                                     cgGTCAGCTTGGCCACGACTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCATCTACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATc                                                                                                   >  2:694243/1‑151 (MQ=9)
                                                     cgGTCAGCTTGGCCACCACTTCGTCGGTGGTGTCGGTTACTTGGACCGAAGCAGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATc                                                                                                   <  2:1604255/151‑1 (MQ=255)
                                                     cgGTCAGCTTGGCCACCACTTCGTCGGTGGTGTCGGTTACTTGGACCGAAGCAGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATc                                                                                                   <  1:130609/151‑1 (MQ=255)
                                                      ggTCAGCTTGGCCACTACTTCGTCGGTGGTGTCGGTCACTTGGACCGAAGCCGCACCGCCCAGCTCGAGGTTCTCGACCGTAGCGCCGGTGGCGCCTACGGCCGTGCTGATGCCCAGGCTGTTCTGGCTCGCCTCGGTGTAGACCTCAtag                                                                                                  >  2:6364331/1‑149 (MQ=12)
                                                      ggTCAGCTTGGCCACCACTTCGTCGGTGGTGTCGGTTACTTGGACCGAAGCAGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCg                                                                                                  >  1:4159256/1‑151 (MQ=37)
                                                          aGCTTGGCCACTACTTCGTCGGTGGTGTCGGTCACTTGAACAGTAGCCGCACCGCCCAACTGCAGGCTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCtt                                                                                              >  2:800383/1‑151 (MQ=18)
                                                          aGCTTGGCCACGACTTCGTCGAGGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGATAGACGTCATCGCCtt                                                                                              >  1:5139841/1‑151 (MQ=37)
                                                          aGCTTGGCCACGACTTCGTCGAGGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGATAGACGTCATCGCCtt                                                                                              >  2:5974000/1‑151 (MQ=37)
                                                           gCTTGGCCACGACTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCGTCTACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTg                                                                                             >  1:1926663/1‑151 (MQ=12)
                                                             ttGGCCACGACTTCGTCGAGGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCGGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTgcg                                                                                           >  1:4054669/1‑151 (MQ=38)
                                                              tGGCCACCACTTCGTCGGTGGTGTCGGTTACTTGGACCGAAGCAGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTgtgc                                                                                          >  1:5242929/1‑151 (MQ=37)
                                                               ggCCACTACTTCGTCGGTGGTGTCGGTCACTTGGACCGAAGCCGCACCGCCCAACTCAAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCa                                                                                         >  1:2854815/1‑151 (MQ=25)
                                                                 ccaccaCTTCGTCGGTGGTGTCGGTGACTTGAACGGTAGCTGCACCGCCCAACTCAAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGATAGACGTCATCGCCTTGCGCATc                                                                                       >  1:2652503/1‑151 (MQ=255)
                                                                  cacgacTTCGTCGAGGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCGGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCg                                                                                      <  1:308657/147‑1 (MQ=38)
                                                                       cTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCGTCTACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGgt                                                                                 <  2:544185/151‑1 (MQ=16)
                                                                         tcgtcgGTGGTGTCGGTCACTTGAACAGTAGCCGCACCGCCCAGTTGCAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCGTCGTGGgtgt                                                                               <  1:841689/151‑1 (MQ=16)
                                                                            tcgagggtgTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGATAGACGTCATCGCCTTGCGCGTCATGGGTGTACg                                                                            <  1:2174909/146‑1 (MQ=37)
                                                                             cgGTGGTGTCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCGTCTACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCGTTGTAGACGCCATCGCCTTGCGCCGCGTGGGGGTCCgg                                                                           >  2:5581353/1‑151 (MQ=16)
                                                                              ggtggtGTCGGTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAAGGTCGCACCAGTGGCGTCTACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGgc                                                                          >  2:3565359/1‑151 (MQ=255)
                                                                              gagggtgTCGGTGACTTGAACCGAAGCGTTGCCGCCCAATTCCAGGTTCTCGAACGCAGCGCCGGTGGCGTCTACTGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGgc                                                                          <  2:431330/148‑1 (MQ=12)
                                                                                tggtgTCGGTCACTTGGACCGAAGCCGCACCGCCCAACTCAAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGgcgc                                                                        <  1:5041492/151‑1 (MQ=255)
                                                                                 ggtgTCGGTCACTTGAACAGTAGCCGCACCGCCCAGTTGCAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCGTCGTGGGTGTACGgcgca                                                                       >  2:3120231/1‑150 (MQ=255)
                                                                                        gTCACTTGAACCGAAGCGTTGCCGCCCAGTTGCAGGTTCTCGAACGCAGCACCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGttt                                                                <  1:3198751/151‑1 (MQ=12)
                                                                                             ttGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGATAGACGTCATCGCCTTGCGCGTCATGGGTGTACGGCGCGCTGGTTTCGcca                                                           >  1:6133679/1‑150 (MQ=37)
                                                                                                  ccGAAGCAGCACCGCCCAACTCGAGGTTCTCGAGCGCAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCGCGCATCCTGGTAGTCGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCtt                                                      <  2:2909106/151‑1 (MQ=18)
                                                                                                    ttcgcggCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGa                                                    <  2:2433792/145‑1 (MQ=17)
                                                                                                    gAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCGGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGa                                                    <  2:2329283/151‑1 (MQ=255)
                                                                                                       gcggCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATgg                                                 >  1:3570717/4‑151 (MQ=17)
                                                                                                        cgttgccGCCCAATTCCAGGTTCTCGAACGCAGCGCCGGTGGCGTCTACTGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGgcgcg                                                                       >  2:3669918/6‑128 (MQ=9)
                                                                                                        cggCACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGt                                                <  1:430247/149‑1 (MQ=17)
                                                                                                        cagCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGt                                                >  2:3776984/3‑151 (MQ=18)
                                                                                                         cGCACCGCCCAACTCAAGGTTCTCGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTg                                               >  2:1821044/1‑151 (MQ=21)
                                                                                                           cACCGCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGAc                                             >  2:1979240/1‑151 (MQ=17)
                                                                                                             ccGCCCAATTCCAGGTTCTCGAACGCAGCGCCGGTGGCGTCTACTGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCt                                           <  2:4552535/151‑1 (MQ=14)
                                                                                                               gCCCAACTCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTga                                         >  2:3246691/1‑150 (MQ=17)
                                                                                                                 ccAACTCAAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGATAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTCTCGCCGGCCTTGATGGTGACCTGAGc                                       <  2:903508/151‑1 (MQ=21)
                                                                                                                     cTCGAGGTTCTCGAACGCAGCGCCGGTGGCGTCCACCGCCGAGTTGCTGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCGttg                                   <  2:5720782/151‑1 (MQ=25)
                                                                                                                     cTCGAGGTTCTCGAACGCAGCACCAGTGGCATCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGATAGACATCATCGCCTTGCGCGTCGTGGGTGTACGGCGCGCTGGTTTCACCGGCCTTGATGGTGACCTg                                          <  1:93214/144‑1 (MQ=18)
                                                                                                                      tcaaGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGATAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTCTCGCCGGCCTTGATGGTGACCTGAGCGttgt                                  >  2:2799644/4‑151 (MQ=18)
                                                                                                                      tCGAGGTTCTCGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACATGAGCActgt                                  >  2:5830020/1‑151 (MQ=14)
                                                                                                                      tCGAGGTTCTCGAACGCAGCACCAGTGGCATCCACCGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGATAGACATCATCGCCTTGCGCGTCGTGGGTGTACGGCGCGCTGGTTTCACCGGCCTTGATGGTGACCTGGGCAttgt                                  >  2:6481410/1‑151 (MQ=18)
                                                                                                                       gcaGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCgttgtt                                 <  1:2267587/149‑1 (MQ=9)
                                                                                                                        caGGTTCTCGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCGttgttg                                >  1:4808029/2‑151 (MQ=9)
                                                                                                                           gTTCTCGAACGCAGCGCCGGTGGCGTCTACTGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCATTGTTGCTc                             >  1:2766723/1‑151 (MQ=14)
                                                                                                                           gTTCTCGAACGCAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTc                             <  1:4435703/151‑1 (MQ=17)
                                                                                                                            ttCTCGAACGCAGCGCCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCATTGTTGCTCa                            >  1:1798797/1‑151 (MQ=14)
                                                                                                                            ttCTCGAACGCAGCACCAGTGGCGTCTACTGCCGAGTTGATGCCCAGGCTGATCTGGCCAGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCACCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCa                            >  1:3450286/1‑151 (MQ=11)
                                                                                                                              ctcGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAgg                          >  1:1120864/1‑151 (MQ=14)
                                                                                                                              ctcGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCGTTGTTGCTCAgg                          <  1:5470300/151‑1 (MQ=9)
                                                                                                                               tcGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAggg                         <  2:4927579/151‑1 (MQ=14)
                                                                                                                               tcGAACGTTGCGCCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCGTTGTTGCTCAggg                         >  2:127971/1‑151 (MQ=9)
                                                                                                                               tcGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGCTTCGCCGGCCTTGATGGTGACCTGGGGCTTGTTGCTCAggg                         >  2:5071199/1‑151 (MQ=17)
                                                                                                                               tcGAACGTAGCGCCGGTGGCGTCTACCGCCGAATTGATACCACGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAggg                         <  1:3351121/151‑1 (MQ=14)
                                                                                                                               tcGAACGCAGCGCCGGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAggg                         >  2:4053296/1‑151 (MQ=255)
                                                                                                                                cGAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCATTGTTGCTCAGGGt                        >  1:64439/1‑151 (MQ=14)
                                                                                                                                 gaacgttgcgcCGGTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAGGGTc                       >  1:1680330/8‑151 (MQ=9)
                                                                                                                                 gAACGCAGCACCAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCATTGTTGCTCAGGGTc                       <  1:493684/151‑1 (MQ=14)
                                                                                                                                        gcaccAGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCATTGTTGCTCAGGGTCACGATCa                <  2:5295147/148‑1 (MQ=11)
                                                                                                                                             ggTGGCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAGGGTCACGACCAAGTCg           >  1:3926837/2‑151 (MQ=11)
                                                                                                                                             aGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCATTGTTGCTCAGGGTCACGATCAGGTCg           <  1:668416/151‑1 (MQ=11)
                                                                                                                                             aGTGGCGTCAACCGCCGAGTTGATGCCCAGACTGATCTGGCCCGCATCCTGGTAGACGTCATCGCCCTGTGCGTCGTGGGTGTACGGCGCGCTGGTTTCGCCAGCCTTGATGGTGACCTGAGCATTGTTGCTCAGGGTCACGATCAGGTCg           >  2:1384817/1‑151 (MQ=11)
                                                                                                                                                ggCGTCTACCGCCGAATTGATGCCCAGGCTGATCTGGCCCGCATCGTTGTAGACGTCATCGCCTTGCGCCGCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGAGCGTTGTTGCTCAGGGTCACGACCAGGTCGTGa        >  1:1191102/1‑151 (MQ=14)
                                                                                                                                                 gCGTCTACCGCCGAATTGATACCAAGGCTGATCTGGCCAGCATCGTTGTAGACGTCATCGCCTTGCGCGTCGTGGGTGTACGGCGCACTGGTTTCACCGGCCTTGATGGTGACCTGGGCGTTGTTGCTCAGGGTCACGACCAAGTCGTggg       >  2:2666649/1‑149 (MQ=2)
                                                                                                                                                      |                                                                                                                                                      
GGTCAGCGTGATGGTGTAAGTGATCTCGCCGCCTTCGGTGACCGAAGGGGTCGCGGTCAGCTTGGCCACCACTTCGTCGGTGGTGTCGGTCACTTGAACCGAAGCCGCACCGCCCAACTCGAGGTTCTCGAACGCAGCGCCAGTGGCGTCCACCGCCGAGTTGATGCCCAGGCTGATTTCGCCCGCATCCTGGTAGACGTCATCGCCTTGCGCATCGTGGGTGTACGGCGCGCTGGTTTCGCCGGCCTTGATGGTGACCTGGGCATTGTTGCTCAGGGTCACGACCAGGTCGTGAGCCAGC  >  CP009974/1895971‑1896271

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: