Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,738,4480AC38.3% 85.9 / 14.3 47V84G (GTG→GGG) trmDtRNA (guanosine(37)‑N1)‑methyltransferase TrmD
Reads supporting (aligned to +/- strand):  ref base A (24/5);  new base C (0/18);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 7.37e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATC  >  NZ_CP009273/2738201‑2738693
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
tGCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCAcc                                                                                                                                                                                                                                                     >  1:539429/1‑250 (MQ=255)
        cTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGcc                                                                                                                                                                                                                                             >  1:302020/1‑250 (MQ=255)
            aTGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCa                                                                                                                                                                                                                                         <  1:434175/250‑1 (MQ=255)
            aTGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCa                                                                                                                                                                                                                                         <  2:9349/250‑1 (MQ=255)
               tCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCAccc                                                                                                                                                                                                                                      >  2:323949/1‑250 (MQ=255)
                     atcccccGGAATAAACCGGGAAACGGAGTCATTAGGCGTTATTGCTGTAAACCCCCCCCCCCGGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCCCCCGCTCATCTATCCCTTCGTAGCGCCCGCACCCCAGAATCAATTTTTGATTCGTTGCCATTTCGCTGACCCCCGCTTGATCAAGCTGGCGTCCGT‑GTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCcgcg                                                                                                                                                                                                                                <  2:362105/246‑1 (MQ=255)
                                      gggAAACGGAGTCAATCGGCGTCATTGCTGGTAACTCACCACCACGGAGAACGTAATCGCCGATTGACCATTCTTCTTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCGGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAAt                                                                                                                                                                                                               <  1:566924/250‑1 (MQ=255)
                                       ggAACCGGGTTCAACCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCGACCGCACCCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCCCCTTTGCGCTTTCACCCGCCGCGGCTTTTGCTGCATGAATg                                                                                                                                                                                                              <  2:40535/250‑1 (MQ=255)
                                               aGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCg                                                                                                                                                                                                      >  1:302090/1‑250 (MQ=255)
                                                     tCAGCTTTATTGTTGGTAACTCACACCCACGGAGACGGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTACGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCACTTGATCAAGCTTGCGTCCCT‑GGGGTGACAAAAAAATCCCCTTGGCGCCTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAggg                                                                                                                                                                                                <  2:279116/250‑1 (MQ=255)
                                                       aGCGTCATTGCTGGTAACCCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTAACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGtt                                                                                                                                                                                              <  1:72726/250‑1 (MQ=255)
                                                        gcgtagttGCGGGTAACTCCCGACCGATGAAAGCTTAAGCGCCGATTGACTATGCTTCGTAAATTTCGGTTTGGACAACGCGCTAATCTATACCTTCGTAGCGACCGCACACAGGAATCAATTTTGGATGCTTTGCCGGTTCGCTGACGCCCGCTTGATCAAGCTGGCTTCCCT‑GTGGGGACAGAAAAAGCCCCTTTGCGCCTTACCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTg                                                                                                                                                                                             <  2:374711/244‑1 (MQ=255)
                                                             aTTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACca                                                                                                                                                                                        >  1:89632/1‑250 (MQ=255)
                                                               tGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcat                                                                                                                                                                                       >  2:387876/1‑249 (MQ=255)
                                                                           ccccccacgtagaACGTAATCCCCATTTGACGATCCTTCCTCAAGTTCGGTTTGGATCACGCGCTCATCTATGCCTTCTTAGCGCCCCCACACCAGAATCAATTTTTGATTCGTTGCCATTTCGCTGACGCCCGCTTGATCACGCTTGCGTCCCG‑GTGTTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATcccc                                                                                                                                                                          <  2:339944/240‑1 (MQ=255)
                                                                                 acgaaaacagTAATCGCCGATCAACCATTCTTCGTAAATTGCGGTTGGGCACCCGCCCTAATCTATCCCTTGGTGGCGCCCGCACACCAGAATCAATTTTTGAGTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GGGGTAACAGAAAAACCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccg                                                                                                                                                                    <  2:481300/241‑1 (MQ=255)
                                                                                      gaACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgccgcc                                                                                                                                                               >  1:302293/1‑250 (MQ=255)
                                                                                        aCGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgccgccg                                                                                                                                                              >  1:122559/1‑249 (MQ=255)
                                                                                               cGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACg                                                                                                                                                      <  1:60136/250‑1 (MQ=255)
                                                                                                      tGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTc                                                                                                                                                >  2:508907/1‑249 (MQ=255)
                                                                                                          cATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacgg                                                                                                                                           >  1:266444/1‑250 (MQ=255)
                                                                                                          cATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGGGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacgg                                                                                                                                           <  1:309500/250‑1 (MQ=255)
                                                                                                                   tCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGcc                                                                                                                                  >  1:554780/1‑250 (MQ=255)
                                                                                                                     aaTTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCgg                                                                                                                                >  1:445487/1‑250 (MQ=255)
                                                                                                                        ttCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGGATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCa                                                                                                                             >  2:400676/1‑250 (MQ=255)
                                                                                                                           ggTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGc                                                                                                                          >  1:20480/1‑250 (MQ=255)
                                                                                                                                    cACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGCGTTGCACCATCATTACCATCCCCGGTCCGCCGCCGTAAGGACGATCGCCCACGGTACGGTGCCGGTCATGCGTGAAGTcg                                                                                                                 >  1:162694/1‑250 (MQ=255)
                                                                                                                                     aCGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTcgc                                                                                                                >  2:74426/1‑250 (MQ=255)
                                                                                                                                           tCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGAc                                                                                                          >  2:383293/1‑250 (MQ=255)
                                                                                                                                                   aCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCACGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGctc                                                                                                  >  2:600279/1‑250 (MQ=255)
                                                                                                                                                     cTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTg                                                                                                >  1:492864/1‑250 (MQ=255)
                                                                                                                                                                ccGCACACCAGAATCATTTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTcagca                                                                                     <  1:233963/250‑1 (MQ=255)
                                                                                                                                                                ccGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGATGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACCCCAGCTCTGGACGCTccg                                                                                       >  1:314657/1‑246 (MQ=255)
                                                                                                                                                                ccGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTcag                                                                                       >  2:228593/1‑248 (MQ=255)
                                                                                                                                                                        aagacacattttttGTTTCGTTGCCTGTTCGCTCACGCCCCTTTGAGCAAG‑TGGGGGCCCTGGTGGTGACAAATAACCCCCCTTTCGCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCAttt                                                                             <  2:515153/241‑1 (MQ=255)
                                                                                                                                                                             tCAATTTTTGATTCGTTGCCGGTTCGCTGACCCCCGCTTGATCAAGCTTGCGTCCTT‑GTGGTAACAGATAAATCCCCTTTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAAc                                                                        <  1:256723/250‑1 (MQ=255)
                                                                                                                                                                                           gtggCCGGTTCCCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GGGGTGACAGATAAATCCCCTTGGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGtt                                                           <  1:384652/246‑1 (MQ=255)
                                                                                                                                                                                             tGCCATTTCGCTGCCCCCCGCTTGATCAA‑CTTGGCTCCCTGGTGGGAACAAATAAATCCCCTTTGCCCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGtt                                                           <  1:358966/247‑1 (MQ=255)
                                                                                                                                                                                                       cTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCgg                                              <  1:318900/250‑1 (MQ=255)
                                                                                                                                                                                                           cGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAAt                                          >  1:123773/1‑250 (MQ=255)
                                                                                                                                                                                                            gCCCGTTTGATCAAGCTTGCGCCCTT‑GGGGTGACAGATAAACCCCCTTGCCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAAtt                                         <  1:286068/250‑1 (MQ=255)
                                                                                                                                                                                                                      caagcctggctccccT‑GTGGTGACAAATAACCCCCCTTGGCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACa                               <  1:153287/238‑1 (MQ=255)
                                                                                                                                                                                                                                 tCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGtt                                                           >  2:553010/1‑211 (MQ=255)
                                                                                                                                                                                                                                  cccT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTa                                                          <  1:553010/211‑1 (MQ=255)
                                                                                                                                                                                                                                    ctggtgaggaaaGAAAAATCCCCTTTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAg                  <  2:275334/239‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAAATCCCCTTTCCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCaa    <  1:484033/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaTCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATc  >  1:347126/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
TGCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATC  >  NZ_CP009273/2738201‑2738693

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