Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273161,0680AC27.2% 204.8 / 29.2 149intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (50/57);  new base C (5/35);  total (57/92)
Fisher's exact test for biased strand distribution p-value = 1.02e-04
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 151 total aligned reads displayed.

ACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAA  >  CP009273/160918‑161213
                                                                                                                                                      |                                                                                                                                                    
aCCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                      >  1:762729/1‑151 (MQ=255)
aCCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                      >  2:3539597/1‑151 (MQ=255)
    tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                                  >  1:4921597/1‑151 (MQ=255)
     gAGGTGCAAAACGAGATGAATGGGCGTTATCCCACACATGTCTGGCCGGACGACCCGGCTAATACTGCACCGACGCGACGGACGAAATAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTTTGGGAGTGGGCCGGGGTGAGGGCATCAgc                                                                                                                                                 <  2:2155669/151‑1 (MQ=255)
     gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgc                                                                                                                                                 >  2:5124127/1‑151 (MQ=255)
         tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCgcgcgc                                                                                                                                             >  2:5046432/1‑151 (MQ=255)
            aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                          >  1:2023134/1‑151 (MQ=255)
                 gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcac                                                                                                                                     >  2:3988456/1‑151 (MQ=255)
                  agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcc                                                                                                                                    >  1:5032471/1‑146 (MQ=255)
                  agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcc                                                                                                                                    >  1:3199463/1‑146 (MQ=255)
                    attaatgggCGTTACCCCAAAATTGTGTGGCGGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAATGTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTtcaccct                                                                                                                                  <  2:3155097/145‑1 (MQ=25)
                    aTGAATGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAATTTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                  <  2:3155248/151‑1 (MQ=37)
                    aTGAACGGGCGCTAACCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCTCCGCGACGCCCGAAAAAGTATTCGTAAGGTTGGTTTTCCCCATCTCCCTGTGGGCGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGtcaccct                                                                                                                                  <  2:3777185/151‑1 (MQ=255)
                    aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGtggaacca                                                                                                                                  >  2:3301254/1‑144 (MQ=255)
                        aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccca‑cacc                                                                                                                             >  2:3449456/1‑151 (MQ=255)
                        aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtaaccct‑cacc                                                                                                                             >  2:4776249/1‑151 (MQ=255)
                        aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccg‑cacc                                                                                                                             >  1:1189921/1‑151 (MQ=255)
                         aGGGCTTAATCCAAAACATTTCTGGCCGTACTACCCGGAAAAAACTCCACCGACGCTACGGACCAAAAATTATTCTAAAGTTTGTTTTAATCCCACACCCTGTGGGATAGGGCCGGGGTTAGGCCATCAGCGCGCACGTtcaccct‑caccc                                                                                                                            <  2:5256995/151‑1 (MQ=255)
                          gggCGTTATCCCAAACATGTCTGGCCGGCCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTcaccct‑caccct                                                                                                                           <  2:4300358/151‑1 (MQ=255)
                          gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgctccct‑cacccg                                                                                                                           >  2:3225426/1‑150 (MQ=255)
                           ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgcccca‑cacccga                                                                                                                          >  1:3773686/1‑149 (MQ=255)
                            gCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCTTCAGCGCGCACGgtc‑cccc‑caccccgac                                                                                                                        >  2:4188715/1‑147 (MQ=38)
                            gCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCGTCAGCGCGCACGTGCACCCG‑CGCCCaaa                                                                                                                         >  2:3012173/1‑151 (MQ=255)
                            gCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCA‑CCCCCaaa                                                                                                                         >  1:4776025/1‑151 (MQ=255)
                             cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCG‑CACCCAAAc                                                                                                                        >  2:2499161/1‑151 (MQ=255)
                             cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCA‑CACCCGAAc                                                                                                                        >  2:1213711/1‑151 (MQ=255)
                              ggtATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCG‑CACCCaaaac                                                                                                                       >  1:1252616/3‑149 (MQ=255)
                              gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCA‑CAACCAAAcc                                                                                                                       >  2:1847231/1‑151 (MQ=255)
                                tATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGAACCTCCACCCAAAccc                                                                                                                      >  1:300187/1‑151 (MQ=37)
                                  tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCC‑CACCCGCACCcact                                                                                                                   >  2:1348435/1‑151 (MQ=38)
                                   cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGGCCCG‑CCCCCAAAcccacac                                                                                                                  >  1:1010521/1‑147 (MQ=37)
                                     cAAACTTGTCTGGCCGGACGCCCCGGCAAATACGGCCCCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTccc                                                                                                                <  2:1172016/151‑1 (MQ=255)
                                      aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCACCA‑CACCCAAAcccacaacca                                                                                                               >  2:4610824/1‑144 (MQ=38)
                                         cttGTCGGGCGGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTGGGTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCaa                                                                                                            <  2:3542953/149‑1 (MQ=38)
                                          aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCC‑CACCCAAAACCACTaccgccag                                                                                                           >  1:2203375/1‑146 (MQ=37)
                                               tGGCCGGACGCCCCGGCAAATACGGCACCGACGCGACGGACTAAAAATTATTCGTAAGGTGGGTTTTCTCCCTCCCCCGGGGGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGgaga                                                                                                      <  2:634160/151‑1 (MQ=255)
                                               gggCCGGACGCCCCGGCAAATAGTGCACCGACGCGACGGACGAACACGTATTCGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGgaga                                                                                                      <  2:3322409/150‑1 (MQ=38)
                                                ggccggccGACCCGCCAAATACTGCACCGACGCGCCGGACGACAATTTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTTAGGCCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGagag                                                                                                     <  1:785342/151‑1 (MQ=255)
                                                ggCCGGAGGACCCGGCAAATACTGCACCGACGCGACGGACTAAAATGTATTCGTAAGTTTGGTTTTCTCCCTCTCCCTGTGGCAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGagag                                                                                                     <  2:2164232/151‑1 (MQ=255)
                                                 gCCGGACGACCCGGCAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                    <  2:3440637/151‑1 (MQ=255)
                                                  ccGGCCCACCCGGCAAATACTCCACCGACGCGCCGGACCAAAAAGTATTCGTACGGTTGGTTTTCGCCCTCTCCCTGTGGTAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                   <  1:4291872/151‑1 (MQ=255)
                                                  ccGGACGCCCCGCAAAATACGCCACCGCCGCGCGGGACTAAAAATTTTTCTTCGGTTTGTTTTACACCCCCCCCCTGTGGGTGTGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                   <  1:2783612/151‑1 (MQ=255)
                                                  ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                   <  2:4417658/151‑1 (MQ=255)
                                                    ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTATGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                 <  1:3408790/151‑1 (MQ=255)
                                                     gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCACAAGGGAGAGGGGAc                                                                                                >  2:3314993/1‑151 (MQ=255)
                                                       cgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCACCA‑CCCCCATACCCACACCCCCCAGGGCGTGGGGGCCg                                                                                              >  2:4008066/1‑151 (MQ=25)
                                                           ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCA‑CACCCAAAAACAATACCTCAAGGGAGAGGGGTCCGAtag                                                                                          >  2:1929558/1‑149 (MQ=25)
                                                               cAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑ca                                                                                    <  1:3927647/151‑1 (MQ=255)
                                                               cAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCG‑CACCCAAAACCACTCCCCCCAGGGGGAGGGGTCCGagacagcc                                                                                      >  2:2674879/1‑144 (MQ=25)
                                                               cAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCG‑CACCCAAAACCGCTCCCACCAGGGAGAGGGGGCCGGTCGAG‑‑cc                                                                                    >  2:3188303/1‑150 (MQ=37)
                                                                aaaTACGGCACCTACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑cac                                                                                   <  2:4256473/151‑1 (MQ=255)
                                                                 aaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCG‑CAGCCGAAACCACAACCGCAATGGTGGGGGGACCCATCGAG‑‑caca                                                                                  >  1:2193659/1‑151 (MQ=25)
                                                                 aaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGACCCA‑CACCCCAACCCACACCCACAAAGGGGGGGGGTGCGATCGAG‑‑ccca                                                                                  >  1:1716061/1‑151 (MQ=25)
                                                                  aaaCTGCACCGACGCGCGGGACTAAAATTTTTTCGTTCGGTTGTTTTTCCCCCTCCCCCGGTGGGAGAGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACaa                                                                                 <  1:4136774/149‑1 (MQ=255)
                                                                    aCTGCACCGACGCGACGGACGAAAAAGTATGCGTTAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAAt                                                                               <  2:2748654/151‑1 (MQ=255)
                                                                     cTGAACCGACGCTACGTACTAAAATGTTTTCGTATGGTTGTTTTACCCCCTCTCCCTGTGGTTGTGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAAtt                                                                              <  1:1457124/151‑1 (MQ=255)
                                                                      tGCACCTACGCGACGGACGAAAATGTTTTCGTACGGTGGGTTTTCGCCCGCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAAttt                                                                             <  2:454473/151‑1 (MQ=255)
                                                                      tGCACCGACGCGACGGACGAAAAAGTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAAttt                                                                             <  1:799945/151‑1 (MQ=255)
                                                                        cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCT‑CACCCAAACCCACTCCCTCAAGGGAGTGGGGACCGATCGAG‑‑CACAAAATTTg                                                                           >  2:1774733/1‑151 (MQ=255)
                                                                         aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCA‑CACCCAAAACCACACCCACAAGGGGGAGGGGGACGATCGAG‑‑CACAAATTTTgt                                                                          >  1:3413721/1‑150 (MQ=25)
                                                                          ccgccgCCACGGACGCCACAGTATCCGTAAGGATGGTTTTCCCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCT‑CACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTgag                                                                         <  2:4976648/151‑1 (MQ=37)
                                                                          ccGACGCTACGGACGAAAAAGTTTTCGTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGATAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTgag                                                                         <  2:3314660/151‑1 (MQ=255)
                                                                          ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCA‑ACCCCCAACACACAAACAAAACGGGGAGGGGACCGAACGAG‑‑CACAAAATTTtag                                                                         >  1:1071927/1‑151 (MQ=21)
                                                                           cgtcgcgacgGACCAAAATTTTTGCTTATGTTTGTTTTTCACCCTCTCCCTGTGGGATTGGGCGGGGTTGTGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTgaga                                                                        <  1:4594588/148‑1 (MQ=255)
                                                                           cgacgcgacgGACGAAAATTTATTCTTACTGTTGGTTTTCTCCCACTCCCTGGGGGTGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTgaga                                                                        <  1:851945/151‑1 (MQ=255)
                                                                              cgcgaggGCCTAAAATGTTTTTGTTTGTTTGGTTTTCCCCCTCTCCCTGTGGTAGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGat                                                                     <  1:2764599/151‑1 (MQ=255)
                                                                              cgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCA‑CACCCAAAACCACAACCGCACAGGAGAGGGGACCGATAGAGCCCAAAAATTTTgagat                                                                       >  2:1992409/1‑150 (MQ=21)
                                                                              cgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCA‑CACCCAAACCCACTCCCACAAGGGAGATGGGACCGAAAGAG‑‑CACAAAATTTTAGatat                                                                     >  1:3931518/1‑151 (MQ=25)
                                                                              cgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCAACCA‑CACCCAAAACCACTCCCACAAGGGAGAGGGGACCCATCGAT‑‑CACAAAATTTTAGatat                                                                     >  1:3927310/1‑151 (MQ=25)
                                                                               gctagggacgAAAAAGTATTCGTATGGTTGTTTTTCTCCCTCCCCGTGTGGTTGAGGGCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGata                                                                    <  1:4526109/146‑1 (MQ=255)
                                                                                  ccggactAAAATTTTTTCGTATGTTTGTTTTTCCCCCTCTCCCTGGGGTTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATAtct                                                                 <  1:574169/150‑1 (MQ=255)
                                                                                    ggacgACAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATAtcttc                                                               <  2:743817/151‑1 (MQ=255)
                                                                                         accaaGAATTCGTAAGGTTGGTTTTCTCCCCCTCCCTGGGGGAGAGGGCCGGGTGGAGGGCATCCGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAGCGAG‑‑CACAAATTTTGAGAGATAACTTCTTCTg                                                          <  2:5186384/148‑1 (MQ=255)
                                                                                              gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCAACCT‑CACCCAAACCCTATCCCTAAATGGTGTTGGGAACGATAGAT‑‑CACAAATTTTTAGATATATATTCTTCTGTcgtt                                                     >  2:2159732/1‑148 (MQ=21)
                                                                                                 ttgttacggttggttTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                  <  2:3673070/144‑1 (MQ=255)
                                                                                                 ttCTTACGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                  <  2:3673147/151‑1 (MQ=255)
                                                                                                    gTATGGTTGGTTTTCGCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                               <  1:4250084/151‑1 (MQ=255)
                                                                                                      acggttggttCTCTCCCTCTCCCTCGGGGAGTGGGCCGGGTTGCGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                             <  2:2203386/149‑1 (MQ=255)
                                                                                                      aaGGTTGGTTTGCGCCCTCTCCCTGTGGGATTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                             <  1:5076766/151‑1 (MQ=255)
                                                                                                       cggggggttTGCCCCCTCGCCCTGTGGGCGGGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CCCAAATTTTGAGAGCTATCTTCTTCTGTCTTGTAACagaac                                            <  1:5263269/146‑2 (MQ=255)
                                                                                                         gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                          <  2:996862/151‑1 (MQ=255)
                                                                                                              ttttCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaa                                     <  1:3971035/151‑1 (MQ=255)
                                                                                                              ttttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCAACCT‑CACCCAAAACCACTCCCGCAACGGGGGGGGGACCCCTGGCG‑‑CACACATTTTGTGAGTTATCATATTAAGTCTGGTATCAgaagaaaataa                                     >  2:2808938/1‑146 (MQ=12)
                                                                                                                 tctccccctcccTGGGGGAGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                  <  1:2546030/151‑1 (MQ=255)
                                                                                                                     cctctccctGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                              <  2:1819295/151‑1 (MQ=255)
                                                                                                                      ctctccctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGc                             <  2:2285965/151‑1 (MQ=255)
                                                                                                                      ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGc                             <  2:3799738/151‑1 (MQ=255)
                                                                                                                           cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCGGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                        <  2:2157411/151‑1 (MQ=255)
                                                                                                                            ctGTGGGATTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                       <  1:145839/151‑1 (MQ=255)
                                                                                                                            ctGTGGGAGTGGGCCGGGTGGTGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                       <  2:2329973/151‑1 (MQ=255)
                                                                                                                                gggAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                   <  1:3324076/151‑1 (MQ=255)
                                                                                                                                gggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAATTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                   <  2:5157344/151‑1 (MQ=255)
                                                                                                                                 ggAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                  <  1:3332743/151‑1 (MQ=255)
                                                                                                                                 ggAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCT‑CACCCAAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                  <  1:3614613/151‑1 (MQ=255)
                                                                                                                                  gagcGGGCCGGGGTTCGGGCCTCCGCGCGCACGGTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGCTCGAG‑‑CACAAATTGTGAGAGATATCGTCTTCTGTCTTGTAACAGAAGAACAGAACCTCGGGCTTTTGCGCCTGa                 <  2:2564491/151‑1 (MQ=255)
                                                                                                                                  gagaGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                 <  1:4201481/151‑1 (MQ=255)
                                                                                                                                      gggCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat             <  1:1945996/151‑1 (MQ=255)
                                                                                                                                        gCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg           <  1:3059331/151‑1 (MQ=255)
                                                                                                                                             ggTGAGGGCCTCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag      <  2:4483653/151‑1 (MQ=255)
                                                                                                                                             ggTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag      <  1:1572449/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCT‑CACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑‑CACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAA  >  CP009273/160918‑161213

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