Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,892,771 G→T 5.7% intergenic (+81/+39) PP_5164 → / ← plpB penicillin amidase family protein/NLPA lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,892,7710GT 5.7% 521.2 / 13.2 158intergenic (+81/+39)PP_5164/plpBpenicillin amidase family protein/NLPA lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (80/69);  new base T (2/7);  total (82/76)
Fisher's exact test for biased strand distribution p-value = 8.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 162 total aligned reads displayed.

CAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTTGTTGTC  >  NC_002947/5892622‑5892920
                                                                                                                                                     |                                                                                                                                                     
cAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAg                                                                                                                                                       <  2:1018222/150‑1 (MQ=255)
  gTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGgt                                                                                                                                                     <  2:2501127/150‑1 (MQ=255)
  gTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGgt                                                                                                                                                     <  2:2044145/150‑1 (MQ=255)
     tCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTc                                                                                                                                                  >  1:4264135/1‑150 (MQ=255)
      cAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCg                                                                                                                                                 >  1:5268572/1‑150 (MQ=255)
        gCGGGCGCTGATGGGGGTGATTCCGGCTCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCTGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGc                                                                                                                                               <  1:4611657/150‑1 (MQ=255)
          gggCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGt                                                                                                                                             >  2:3056138/1‑150 (MQ=255)
          gggCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGt                                                                                                                                             >  2:3056193/1‑150 (MQ=255)
           ggCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGtt                                                                                                                                            <  2:927375/150‑1 (MQ=255)
           ggCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGtt                                                                                                                                            <  1:3390906/150‑1 (MQ=255)
            gcgcAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTg                                                                                                                                           <  1:1276124/150‑1 (MQ=255)
             cgcAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTgg                                                                                                                                          >  1:2199653/1‑150 (MQ=255)
              gcAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGc                                                                                                                                         >  2:3361856/1‑150 (MQ=255)
                      gggTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCa                                                                                                                                 >  2:4429593/1‑150 (MQ=255)
                           aTTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGt                                                                                                                            >  1:4441576/1‑150 (MQ=255)
                           aTTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGt                                                                                                                            <  2:462806/150‑1 (MQ=255)
                              ccGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGc                                                                                                                         <  2:1534353/150‑1 (MQ=255)
                               cGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTggt                                                                                                                        >  1:4551317/1‑148 (MQ=255)
                               cGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCt                                                                                                                        >  1:2338685/1‑150 (MQ=255)
                               cGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGGTGGCACAGATCAGGCGTTGCt                                                                                                                        >  1:1015940/1‑150 (MQ=255)
                                 gcgcAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAg                                                                                                                      <  1:4510648/150‑1 (MQ=255)
                                 gcgcAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAg                                                                                                                      <  2:1529053/150‑1 (MQ=255)
                                 gcgcAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAg                                                                                                                      >  1:4978764/1‑150 (MQ=255)
                                   gcAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCg                                                                                                                    <  2:4978764/150‑1 (MQ=255)
                                     aaaGTACCTTACGGCTGGTACCCCCGAAGTATGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGaa                                                                                                                  <  1:5046390/150‑1 (MQ=255)
                                     aaaGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGaa                                                                                                                  <  1:3356580/150‑1 (MQ=255)
                                     aaaGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGaa                                                                                                                  >  2:3942117/1‑150 (MQ=255)
                                      aaGTACCTTACGGCTGGTACCCCCGAGGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGaaa                                                                                                                 >  1:5301317/1‑150 (MQ=255)
                                            cTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGaa                                                                                                           >  1:3704330/1‑150 (MQ=255)
                                             ttACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAAc                                                                                                          >  2:2259195/1‑150 (MQ=255)
                                               aCGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGc                                                                                                        >  1:3765970/1‑150 (MQ=255)
                                                cGGCTGGTACCCCCGACGTACGCGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGcc                                                                                                       <  1:2164191/150‑1 (MQ=255)
                                                  gCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCgg                                                                                                     >  2:1429641/1‑150 (MQ=255)
                                                  gCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCgg                                                                                                     <  2:239181/150‑1 (MQ=255)
                                                    tGGTACCCCCGATGTATGTGGCTGATTTGCGGCCGCTTCGCTGCAAATGGCTTCTACGGTGCTGTGGGGGCAGGGATTAATAAAAACGCCGACACAAGGTTTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCa                                                                                                   <  1:1532799/150‑1 (MQ=255)
                                                    tGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCa                                                                                                   <  1:3964110/150‑1 (MQ=255)
                                                      gTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGt                                                                                                 >  2:1136571/1‑150 (MQ=255)
                                                       tACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTa                                                                                                >  1:4495330/1‑150 (MQ=255)
                                                           cccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGc                                                                                            >  1:1440485/1‑150 (MQ=255)
                                                           cccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGc                                                                                            >  2:5271499/1‑150 (MQ=255)
                                                            ccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCa                                                                                           <  2:2380726/150‑1 (MQ=255)
                                                                 gTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTg                                                                                      <  1:3146660/150‑1 (MQ=255)
                                                                  tAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTgg                                                                                     <  2:3462801/150‑1 (MQ=255)
                                                                        gctgctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTActtc                                                                               >  1:2707883/1‑150 (MQ=255)
                                                                        gctgctTTGCGGCCGCTTCGCAGAAAAAGGCTTCTAGGGTGCTGTGCGGGCAGGGCTGAATAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTActtc                                                                               <  1:1488701/150‑1 (MQ=255)
                                                                           gctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTTAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTg                                                                            <  1:1628922/150‑1 (MQ=255)
                                                                                   gCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGc                                                                                                                     >  2:2627857/1‑101 (MQ=255)
                                                                                    ccGATTCGCAGCAAAAGGATTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCaatgaat                                                                   <  1:1519633/150‑1 (MQ=255)
                                                                                       cTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGcc                                                                <  1:3497189/150‑1 (MQ=255)
                                                                                       cTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGGTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGcc                                                                >  2:1621720/1‑150 (MQ=255)
                                                                                        ttCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTAGTTCTTGGCAATGAATGCCt                                                               >  2:4820060/1‑150 (MQ=255)
                                                                                        ttCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCt                                                               >  1:2543767/1‑150 (MQ=255)
                                                                                        ttCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGAAGACACAAGGTGTAGTCGTTGGAACATATCAGGCGTTGCTAGCGAAAACAGAACGCCGGCAGTACGGCACCTTGGTACATATTGGAAATGAATGACt                                                               >  2:36085/1‑150 (MQ=255)
                                                                                        ttCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAACAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCt                                                               <  1:1047669/150‑1 (MQ=255)
                                                                                         tCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCtt                                                              <  2:1710218/150‑1 (MQ=255)
                                                                                          cGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTg                                                             >  2:630567/1‑150 (MQ=255)
                                                                                           gcagcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGa                                                            <  1:3226427/150‑1 (MQ=255)
                                                                                           gcagcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGa                                                            >  1:1619624/1‑150 (MQ=255)
                                                                                             agcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACt                                                          <  2:2361064/150‑1 (MQ=255)
                                                                                             agcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACt                                                          <  1:4494124/150‑1 (MQ=255)
                                                                                              gcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGcc                                                                <  1:543426/143‑1 (MQ=255)
                                                                                              gcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGcc                                                                >  2:543426/1‑143 (MQ=255)
                                                                                                aaaaGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCgg                                                                                                     >  1:1147917/1‑104 (MQ=255)
                                                                                                aaaaGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCgg                                                                                                     <  2:1147917/104‑1 (MQ=255)
                                                                                                      cTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAATCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTg                                                 >  2:788313/1‑150 (MQ=255)
                                                                                                      cTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTg                                                 >  1:2625816/1‑150 (MQ=255)
                                                                                                       ttCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAAt                                                                                                                >  1:3895296/1‑86 (MQ=255)
                                                                                                       ttCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAAt                                                                                                                <  2:3895296/86‑1 (MQ=255)
                                                                                                       ttCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTgg                                                >  2:4886454/1‑150 (MQ=255)
                                                                                                          tACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATAACAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCATGGTACTTCTTGGCACTGAATGCCTTGACTTCCGGGCTGGTCa                                             >  1:3674330/1‑150 (MQ=255)
                                                                                                            cGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAgg                                           <  1:5154890/150‑1 (MQ=255)
                                                                                                                 cTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCtt                                      <  2:28686/150‑1 (MQ=255)
                                                                                                                  ggtgCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTg                                     <  2:4579721/149‑1 (MQ=255)
                                                                                                                   gtgCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTgg                                    >  1:2628605/1‑150 (MQ=255)
                                                                                                                       gggCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAg                                >  1:2232235/1‑150 (MQ=255)
                                                                                                                         gCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCtt                                                              >  2:4091058/1‑118 (MQ=255)
                                                                                                                         gCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCtt                                                              <  1:4091058/118‑1 (MQ=255)
                                                                                                                             ggCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTg                          <  2:1389994/150‑1 (MQ=255)
                                                                                                                             ggCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTg                          <  1:1202190/150‑1 (MQ=255)
                                                                                                                                tGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGAt                       >  1:1993938/1‑150 (MQ=255)
                                                                                                                                tGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGAt                       >  1:1714101/1‑150 (MQ=255)
                                                                                                                                  aagaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATgg                     >  2:5428365/1‑150 (MQ=255)
                                                                                                                                   agaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCGTGGAATTCCGGGATGTTCAGGGCCTGGGCCTGTTTGTGTATGGc                    >  2:4575911/1‑150 (MQ=255)
                                                                                                                                   agaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAg                                            >  1:5130798/1‑126 (MQ=255)
                                                                                                                                   agaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAg                                            <  2:5130798/126‑1 (MQ=255)
                                                                                                                                       aaaCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCa                <  2:2199653/150‑1 (MQ=255)
                                                                                                                                           gCCGACACCTTGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCa                                             <  1:469388/117‑1 (MQ=255)
                                                                                                                                           gCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCt                                                               <  2:3786288/99‑1 (MQ=255)
                                                                                                                                           gCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCt                                                               >  1:3786288/1‑99 (MQ=255)
                                                                                                                                           gCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGt            >  2:5273912/1‑150 (MQ=255)
                                                                                                                                            ccGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTc           >  1:2140730/1‑150 (MQ=255)
                                                                                                                                              gACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCt         <  2:3600988/150‑1 (MQ=255)
                                                                                                                                              gACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCt         <  2:3396290/150‑1 (MQ=255)
                                                                                                                                                 acaAGGTGTCGGCGTTGGCACATACCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCttgt      >  1:5390367/1‑150 (MQ=255)
                                                                                                                                                   aaGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAg                                            >  2:2244832/1‑110 (MQ=255)
                                                                                                                                                   aaGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGc                                   >  1:823129/1‑119 (MQ=255)
                                                                                                                                                   aaGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGc                                   <  2:823129/119‑1 (MQ=255)
                                                                                                                                                    aGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTtgttgt   >  2:2423242/1‑150 (MQ=255)
                                                                                                                                                     ggTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTTGTTGTc  >  2:1382306/1‑150 (MQ=255)
                                                                                                                                                     ggTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTTGTTGTc  >  2:2905888/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTTGTTGTC  >  NC_002947/5892622‑5892920

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