Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,941,370 T→G 6.3% S187A (TCG→GCG)  PP_2574 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,941,3700TG 6.3% 218.6 / 6.5 81S187A (TCG→GCG) PP_2574LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base G (3/2);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

TCCAGGGCCTGACCGCTGGAACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATGGTCGGT  >  NC_002947/2941221‑2941519
                                                                                                                                                     |                                                                                                                                                     
tCCAGGGCCTGACCGCTGGAACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGtct                                                                                                                                                       <  1:1836106/150‑1 (MQ=255)
 ccAGGGCCTGACCGCTGGAACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGtctc                                                                                                                                                      >  1:2171470/1‑150 (MQ=255)
  cAGGGCCTGACCGCTGGAACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCg                                                                                                                                                     >  2:1590924/1‑150 (MQ=255)
   aGGGCCTGACCGCTGGAACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGc                                                                                                                                                    >  1:906916/1‑150 (MQ=255)
   aGGGCCTGACCGCTGGAACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGc                                                                                                                                                    >  1:81492/1‑150 (MQ=255)
      gCCTGACCGCTGGAACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTgc                                                                                                                                                 <  2:81492/150‑1 (MQ=255)
             cGCTGGAACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCCCCGCTTTGCCAGCCGGGCCCACGTCTCGCTGCGCGAACt                                                                                                                                          <  2:1112715/150‑1 (MQ=255)
              gCTGGAACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCtg                                                                                                                                         >  1:1577700/1‑150 (MQ=255)
                  gAACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCAAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCtcggcc                                                                                                                                     >  2:160379/1‑145 (MQ=255)
                          gACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCg                                                                                                                             >  2:1915863/1‑150 (MQ=255)
                           aCCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCga                                                                                                                            <  1:1590924/150‑1 (MQ=255)
                              tGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCgatga                                                                                                                         >  1:1801394/1‑150 (MQ=255)
                                gCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATc                                                                                                                       <  2:1801394/150‑1 (MQ=255)
                                  aTTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCct                                                                                                                     >  2:921012/1‑150 (MQ=255)
                                     tctcTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCTGGGGGTGATGGTCctgct                                                                                                                  <  2:700408/150‑1 (MQ=255)
                                         tATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGAc                                                                                                              >  2:1460312/1‑150 (MQ=255)
                                         tATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCCGCGCGAc                                                                                                                                            <  1:1548764/120‑1 (MQ=255)
                                         tATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCCGCGCGAc                                                                                                                                            >  2:1548764/1‑120 (MQ=255)
                                           tGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCGCCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGt                                                                                                            >  1:885542/1‑150 (MQ=255)
                                            gACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCCGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTggcagtg                                                                                                           >  2:570519/1‑145 (MQ=255)
                                             cccACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTgg                                                                                                          <  1:1058248/149‑1 (MQ=255)
                                             aCCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCga                                                                                                                            <  2:1442078/132‑1 (MQ=255)
                                             aCCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCga                                                                                                                            >  1:1442078/1‑132 (MQ=255)
                                             aCCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTgg                                                                                                          <  1:803074/150‑1 (MQ=255)
                                               cACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGcc                                                                                                        >  1:1148024/1‑150 (MQ=255)
                                               cACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGcc                                                                                                        >  2:347463/1‑150 (MQ=255)
                                                aCGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGc                                                                                                                               >  2:1885425/1‑126 (MQ=255)
                                                aCGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGc                                                                                                                               <  1:1885425/126‑1 (MQ=255)
                                                  gACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGAc                                                                                                                                            <  2:807455/111‑1 (MQ=255)
                                                  gACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGAc                                                                                                                                            >  1:807455/1‑111 (MQ=255)
                                                     cTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGc                                                                                                  >  1:1482189/1‑150 (MQ=255)
                                                        gATGGCACGATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCCATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCCGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGc                                                                                               <  2:491472/150‑1 (MQ=255)
                                                             cACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTa                                                                                          >  2:1289418/1‑150 (MQ=255)
                                                             cACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTa                                                                                          >  2:51037/1‑150 (MQ=255)
                                                              aCTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTAc                                                                                         >  1:653665/1‑150 (MQ=255)
                                                                tATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACtt                                                                                       <  2:2171470/150‑1 (MQ=255)
                                                                  tCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTg                                                                                     >  1:410046/1‑150 (MQ=255)
                                                                        ccGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCt                                                                               <  1:347463/150‑1 (MQ=255)
                                                                          gAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCTGCACGTACTTTGTCAGCTTg                                                                             <  1:921012/150‑1 (MQ=255)
                                                                            gCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGtt                                                                           >  2:1191972/1‑150 (MQ=255)
                                                                               gCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGCTCCGGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCaa                                                                        <  1:1134981/150‑1 (MQ=255)
                                                                               gCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCACAGGTCGCGCAGCGAGAc                                                                                                                                            <  1:1489301/82‑1 (MQ=38)
                                                                               gCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCACAGGTCGCGCAGCGAGAc                                                                                                                                            >  2:1489301/1‑82 (MQ=38)
                                                                                cTGATGCCGCCGCAGAAGCCCTATATGCTGCTACCCGCGAACCACCGCTTTGCCAGCAAGGCCCAAATCCCGCTGCGAGACCTGTCCACGAAGGATATCCTCCTGCGAGACGTGGCCCCTATCCGTACGTTCTTTGTCAGCTTGTTCAAc                                                                       <  2:402078/150‑1 (MQ=255)
                                                                                   aTGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACgag                                                                    >  1:1488590/1‑150 (MQ=255)
                                                                                      ccgccgCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATg                                                                 <  2:410046/150‑1 (MQ=255)
                                                                                        gccgcAGAAGCCCTATGTGCTGTTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGa                                                                                                                                 >  1:1101184/1‑84 (MQ=255)
                                                                                        gccgcAGAAGCCCTATGTGCTGTTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGa                                                                                                                                 <  2:1101184/84‑1 (MQ=255)
                                                                                         ccgcAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGc                                                              <  2:1644938/150‑1 (MQ=255)
                                                                                         ccgcAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGc                                                              <  2:1148024/150‑1 (MQ=255)
                                                                                            cAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTg                                                           <  1:1191972/150‑1 (MQ=255)
                                                                                               aaGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCa                                                        <  2:1488590/150‑1 (MQ=255)
                                                                                                        gtgCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCTCTTCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATc                                               <  2:906916/150‑1 (MQ=255)
                                                                                                        gtgCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATc                                               <  2:653665/150‑1 (MQ=255)
                                                                                                        ctgctgcTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATc                                               <  1:501321/149‑1 (MQ=255)
                                                                                                          gctgctACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCt                                                            <  2:1492019/135‑1 (MQ=255)
                                                                                                          gctgctACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCt                                                            <  2:1494073/135‑1 (MQ=255)
                                                                                                          gctgctACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCt                                                            >  1:1492019/1‑135 (MQ=255)
                                                                                                          gctgctACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCt                                                            >  1:1494073/1‑135 (MQ=255)
                                                                                                          gctgctACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGt                                             <  1:51037/150‑1 (MQ=255)
                                                                                                          gctgctACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGt                                             >  1:652345/1‑150 (MQ=255)
                                                                                                                  ccGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGtt                                     >  2:1908575/1‑150 (MQ=255)
                                                                                                                   cGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTc                                    >  1:1425392/1‑150 (MQ=255)
                                                                                                                         ccaccGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATg                                                                 >  2:994395/1‑115 (MQ=255)
                                                                                                                         ccaccGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATg                                                                 <  1:994395/115‑1 (MQ=255)
                                                                                                                          caccGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATcc                                                                                                                      >  2:851276/1‑61 (MQ=255)
                                                                                                                          caccGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATcc                                                                                                                      <  1:851276/61‑1 (MQ=255)
                                                                                                                          caccGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCATGCTGGACGTGCGTCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCCCCGtcg                             >  1:885952/1‑150 (MQ=255)
                                                                                                                           accGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATc                                               >  1:2242536/1‑131 (MQ=255)
                                                                                                                           accGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATc                                               <  2:2242536/131‑1 (MQ=255)
                                                                                                                           accGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGtcga                            >  2:1785780/1‑150 (MQ=255)
                                                                                                                            ccGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCGGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGtcgat                           <  2:2334386/150‑1 (MQ=255)
                                                                                                                             cGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGtcgatc                          <  2:1896888/150‑1 (MQ=255)
                                                                                                                                 ttGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCgaga                      <  2:1585926/150‑1 (MQ=255)
                                                                                                                                         cagGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAAc                                                                       <  2:612408/93‑1 (MQ=255)
                                                                                                                                         cagGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAAc                                                                       >  1:612408/1‑93 (MQ=255)
                                                                                                                                            gcacAGGTCGCGCAGCGAGACCTGGGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGc           >  1:2138184/4‑150 (MQ=255)
                                                                                                                                             cccAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGGCGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGc                                                              <  2:388843/98‑1 (MQ=255)
                                                                                                                                             cccAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGGCGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGc                                                              >  1:388843/1‑98 (MQ=255)
                                                                                                                                             cacAGGTCGCGCAGCGAGACCTGGGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGtt                                                                           >  1:1587797/3‑85 (MQ=37)
                                                                                                                                             cacAGGTCGCGCAGCGAGACCTGGGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGtt                                                                           <  2:1587797/83‑1 (MQ=37)
                                                                                                                                               cAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATAATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATg        >  1:785929/1‑150 (MQ=255)
                                                                                                                                                aGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATgg       >  1:409523/1‑150 (MQ=255)
                                                                                                                                                     tcGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATggtcggt  <  2:1577700/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCAGGGCCTGACCGCTGGAACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATGGTCGGT  >  NC_002947/2941221‑2941519

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