Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099741,897,1240CG67.7% 41.5 / 63.2 65intergenic (‑4902/‑3695)RPPX_08470/RPPX_08480tryptophan synthase subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (14/7);  new base G (18/26);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-04
Rejected as consensus: Frequency below/above cutoff threshold.

CTGACAATGCTTACGGTGACGGCATCGATGCTGTCATTGATCGTGGTCTGAGCCGGCGTCGGGTTCGGGGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTGCTACCGACGTAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACGACGACGGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAGGCCGTCTTTGTTGGTCAGCGTGA  >  CP009974/1896986‑1897271
                                                                                                                                           |                                                                                                                                                   
aagtcgccttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGa                                                                                                                                         <  2:3829023‑M1/82‑1 (MQ=255)
 agtcgccttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGAc                                                                                                                                        <  2:2350643‑M1/83‑1 (MQ=255)
  gtcgccttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACg                                                                                                                                       <  2:1112726‑M1/84‑1 (MQ=255)
   tcgccttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGt                                                                                                                                      >  1:2793184‑M1/67‑151 (MQ=255)
     gccttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTg                                                                                                                                    <  2:2296018‑M1/87‑1 (MQ=255)
     gccttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTg                                                                                                                                    >  2:4454686‑M1/65‑151 (MQ=255)
     gccttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTg                                                                                                                                    <  1:1674905‑M1/87‑1 (MQ=255)
     gccttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTg                                                                                                                                    <  1:2689834‑M1/87‑1 (MQ=255)
     gccttcgttgcccggcgtacctggctcgtcgctcacgctagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGAGATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTg                                                                                                                                    >  2:4261901‑M1/65‑151 (MQ=255)
      ccttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCCCCGTCCGCACCTCCACCGGTAGCGATCGCCATCACGACACGCTCCTTCGCCCCGACgtagacctccc                                                                                                                                   <  2:4685872‑M1/88‑6 (MQ=255)
      ccttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGt                                                                                                                                   <  1:3650339‑M1/88‑1 (MQ=255)
       cttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTc                                                                                                                                  <  2:4694940‑M1/89‑1 (MQ=255)
        ttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTcc                                                                                                                                 >  2:172489‑M1/62‑151 (MQ=255)
        ttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTcc                                                                                                                                 <  2:973141‑M1/90‑1 (MQ=255)
        ttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTcc                                                                                                                                 >  2:381796‑M1/62‑151 (MQ=255)
        ttcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTcc                                                                                                                                 >  1:3681516‑M1/62‑151 (MQ=255)
         tcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCCCGACAGGGTCGTTAGCGCCGACCTAGACGTTGTCCg                                                                                                                                <  2:4697896‑M1/91‑1 (MQ=255)
         tcgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCg                                                                                                                                >  2:518422‑M1/61‑151 (MQ=255)
          cgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCgg                                                                                                                               >  2:4078729‑M1/60‑151 (MQ=255)
          cgttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCgg                                                                                                                               >  1:367128‑M1/60‑151 (MQ=255)
           gttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCggg                                                                                                                              >  1:1765338‑M1/59‑151 (MQ=255)
           gttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCggg                                                                                                                              <  2:289193‑M1/93‑1 (MQ=255)
           gttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCggg                                                                                                                              >  1:2595082‑M1/59‑151 (MQ=255)
           gttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCggg                                                                                                                              <  1:4189532‑M1/93‑1 (MQ=255)
           gttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCggg                                                                                                                              <  1:1241289‑M1/93‑1 (MQ=255)
            ttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCgggg                                                                                                                             >  1:3172576‑M1/58‑151 (MQ=255)
            ttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCgggg                                                                                                                             >  1:1517194‑M1/58‑151 (MQ=255)
            ttgcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCgggg                                                                                                                             >  1:4065488‑M1/58‑151 (MQ=255)
            tggccgggcgtgcctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCgggg                                                                                                                             >  2:2857237‑M1/58‑151 (MQ=255)
              gcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCGGGgca                                                                                                                           <  2:84908‑M1/96‑2 (MQ=255)
              gcccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCGGGgca                                                                                                                           >  2:527628‑M1/56‑150 (MQ=255)
               cccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCGGggcag                                                                                                                          <  1:2603191‑M1/97‑3 (MQ=255)
               cccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCGGggcag                                                                                                                          <  2:1371432‑M1/97‑3 (MQ=255)
               cccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCGGggcag                                                                                                                          <  2:4378012‑M1/97‑3 (MQ=255)
                ccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCGgggcagc                                                                                                                         >  1:4404741‑M1/54‑148 (MQ=255)
                ccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCGgggcagc                                                                                                                         <  2:450771‑M1/98‑4 (MQ=255)
                ccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCGgggcagc                                                                                                                         <  2:1190582‑M1/98‑4 (MQ=255)
                ccggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCGgggcagc                                                                                                                         <  2:4877938‑M1/98‑4 (MQ=255)
                 cggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCggggcagcc                                                                                                                        <  2:1482873‑M1/99‑5 (MQ=255)
                 cggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCggggcagcc                                                                                                                        >  1:5395290‑M1/53‑147 (MQ=255)
                 cggcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCggggcagcc                                                                                                                        <  2:2118409‑M1/99‑5 (MQ=255)
                   gcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTccggggcagccac                                                                                                                      >  2:4333750‑M1/51‑145 (MQ=255)
                   gcgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTccggggcagccac                                                                                                                      <  1:4241518‑M1/101‑7 (MQ=255)
                    cgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCCACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGtccggggcagccacg                                                                                                                     >  2:3241936‑M1/50‑144 (MQ=255)
                    cgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAg                                                                                                                                          <  1:964987‑M1/81‑1 (MQ=255)
                    cgtacctggctcgtcggtcacggtagtgctgaccggggttttgtcgagtGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGtccggggcagccacg                                                                                                                     <  1:1486510‑M1/102‑8 (MQ=255)
                                                                                    ttGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCGGGGCAGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGt                                                     >  2:1226697/1‑151 (MQ=21)
                                                                                           catCCGCGCCTTCGACGGTGGCGATCGACATCACAACCGGGTGGTTGCTACCGACGTAGACGTTGTCCGGGGCAGTTGCAGTCGCCGAGCCCGTGGTGCTGTTCGCTGCCACGACGACGGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTc                                              <  2:1921977/149‑1 (MQ=17)
                                                                                                            gTAGCGATCGACATCACGACAGGGTCGTTAGCGCCGACGTAGACGTTGTCCGGGGCAGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCg                             <  1:1226697/151‑1 (MQ=14)
                                                                                                                              accggcgcGTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTgg           >  1:3164042/8‑151 (MQ=11)
                                                                                                                               cgggggCGTTGGTGCCGACGTAGACGTTGTCCGGGGCGGCTACAGTGGCGGAACCGGTGGTGCTGTTGGCTGCCACGACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAGGCCGTCTTTGTTGGt          <  2:2951724/149‑1 (MQ=11)
                                                                                                                               ccggcgcGTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGt          <  1:5428388/145‑1 (MQ=11)
                                                                                                                               ccggcgcGTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGt          <  1:5428107/145‑1 (MQ=11)
                                                                                                                                 ggcgcGTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGTCa        >  1:2057810/5‑151 (MQ=11)
                                                                                                                                 ggcgcGTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGTCa        <  1:2546570/147‑1 (MQ=11)
                                                                                                                                 ggcgcGTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGTCa        >  2:1504616/5‑151 (MQ=11)
                                                                                                                                  gcgcGTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGTCAg       >  1:2266144/4‑151 (MQ=11)
                                                                                                                                   cgcGTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGTCAGc      >  2:1195128/3‑151 (MQ=11)
                                                                                                                                    gcGTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGTCAGCg     >  1:5112206/2‑151 (MQ=11)
                                                                                                                                     cGTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGTAAGCgg    >  1:3177774/1‑150 (MQ=11)
                                                                                                                                     cGTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGATCCGAGAGGTTGTTGATCGGCAAGCAGTCTTTGTTGGTCAGCGt    >  1:3179622/1‑151 (MQ=11)
                                                                                                                                      gTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGTCAGCGTg   >  2:448176/1‑151 (MQ=11)
                                                                                                                                      gTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGTCAGCGTg   >  1:1307304/1‑151 (MQ=11)
                                                                                                                                      gTTGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGTCAGCGTg   >  1:3439902/1‑151 (MQ=11)
                                                                                                                                       ttGGTACCGACATAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACAACGACAGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAAGCCGTCTTTGTTGGTCAGCGTGa  <  1:3082484/151‑1 (MQ=11)
                                                                                                                                           |                                                                                                                                                   
CTGACAATGCTTACGGTGACGGCATCGATGCTGTCATTGATCGTGGTCTGAGCCGGCGTCGGGTTCGGGGTCAGTTGCTCGAACTTGCCAACGTCCGCACCTTCAACGGTAGCGATCGACATCACGACAGGGTCGTTGCTACCGACGTAGACGTTGTCCGGGGCGGTTGCAGTCGCCGAACCGGTGGTGCTGTTGGCTGCCACGACGACGGTGGTGCCGTCGGTCAGCTTGAAGTACAGCTCCGAGTGGTTGTTGATCGGCAGGCCGTCTTTGTTGGTCAGCGTGA  >  CP009974/1896986‑1897271

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