Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,078,665 G→A 5.0% T192T (ACC→ACT)  sdaC ← serine:H+ symport permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,078,6650GA 5.0% 302.9 / 4.2 99T192T (ACC→ACT) sdaCserine:H+ symport permease
Reads supporting (aligned to +/- strand):  ref base G (42/52);  new base A (2/3);  total (44/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GCGCCGTAACGGCGCTTCTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCGAGGATCAGCACGAATGACAGGATTGCC  >  NC_002947/4078516‑4078812
                                                                                                                                                     |                                                                                                                                                   
gcgcCGTAACGGCGCTTCTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGCGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCg                                                                                                                                                     >  2:1898820/1‑150 (MQ=255)
  gcCGTAACGGCGCTTCTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGt                                                                                                                                                   >  1:1040665/1‑150 (MQ=255)
       aaCGGCGCTTCTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGt                                                                                                                                              >  2:1952689/1‑150 (MQ=255)
         cGGCGCTTCTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGc                                                                                                                                            <  1:993894/150‑1 (MQ=255)
          ggCGCTTCTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCg                                                                                                                                           >  1:1805371/1‑150 (MQ=255)
          ggCGCTTCTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCg                                                                                                                                           >  1:214281/1‑150 (MQ=255)
          ggCGCTTCTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCg                                                                                                                                           <  2:381605/150‑1 (MQ=255)
             gcTTCTGGTCGACGGCACAGGCCGAGATGATCGGCGAATCGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGCGTGCAGGAAGACCGAGGCAGGCACCACCTCTGTGGCCCTGTCCAGAATGCCACCGGCCCAGTGCCGCa                                                                                                                                        >  2:130747/1‑150 (MQ=255)
              cTTCTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCAc                                                                                                                                       <  1:2368981/150‑1 (MQ=255)
               ttCTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCAcc                                                                                                                                      <  2:1805371/150‑1 (MQ=255)
                 cTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAg                                                                                                                                    <  2:709799/150‑1 (MQ=255)
                 cTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAg                                                                                                                                    <  1:1417331/150‑1 (MQ=255)
                     tCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTAc                                                                                                                                >  1:344749/1‑150 (MQ=255)
                     tCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTAc                                                                                                                                <  2:2390259/150‑1 (MQ=255)
                     tCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTAc                                                                                                                                >  2:2350928/1‑150 (MQ=255)
                       gACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACa                                                                                                                               <  1:1611956/149‑1 (MQ=255)
                       gACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACa                                                                                                                               >  2:1611956/1‑149 (MQ=255)
                       gACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAg                                                                                                                              <  1:1334104/150‑1 (MQ=255)
                       gACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAg                                                                                                                              >  2:1694751/1‑150 (MQ=255)
                         cGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGc                                                                                                                            <  2:801341/150‑1 (MQ=255)
                          ggCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAgg                                                                                                                             <  2:1247866/148‑1 (MQ=255)
                          ggCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAgg                                                                                                                             >  1:1247866/1‑148 (MQ=255)
                            cAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGACGCCGGCGGCACCTCTGCGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCa                                                                                                                         <  2:146355/150‑1 (MQ=255)
                             aaaGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAg                                                                                                                        >  1:1677304/1‑150 (MQ=255)
                                 gCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAAc                                                                                                                    <  2:1361455/150‑1 (MQ=255)
                                 gCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAAc                                                                                                                    >  2:2275617/1‑150 (MQ=255)
                                   cGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGc                                                                                                                  >  1:1845191/1‑150 (MQ=255)
                                   cGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGc                                                                                                                  <  2:344749/150‑1 (MQ=255)
                                     agaTGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCa                                                                                                                <  2:1686343/150‑1 (MQ=255)
                                           tCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGcga                                                                                                          <  2:1845191/150‑1 (MQ=255)
                                           tCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGcga                                                                                                          >  1:1666711/1‑150 (MQ=255)
                                             ggCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGcgacg                                                                                                        >  1:2264797/1‑150 (MQ=255)
                                              gCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCTCAGCAACGTCAGAGcgacga                                                                                                       <  2:1040665/150‑1 (MQ=255)
                                                gAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATg                                                                                                     <  2:2264797/150‑1 (MQ=255)
                                                   tGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAAc                                                                                                  <  1:2275617/150‑1 (MQ=255)
                                                   tGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAAc                                                                                                  <  1:1309918/150‑1 (MQ=255)
                                                   tGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAAc                                                                                                  <  2:1761147/150‑1 (MQ=255)
                                                     gTTGAACGAGAACACCATCACCGGAATGGCCAGCTACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACgg                                                                                                <  1:1169967/150‑1 (MQ=255)
                                                        gAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTa                                                                                             <  2:299151/150‑1 (MQ=255)
                                                           cGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGCGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGAc                                                                                          <  1:1898820/150‑1 (MQ=255)
                                                             agaACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACca                                                                                        <  2:1590997/150‑1 (MQ=255)
                                                               aaCACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACcagc                                                                                      >  1:601188/1‑150 (MQ=255)
                                                               aaCACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACcagc                                                                                      >  1:1772639/1‑150 (MQ=255)
                                                                 cacCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACcagcag                                                                                    >  1:1899910/1‑150 (MQ=255)
                                                                  acCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTAc                                                                                                                                >  1:527577/1‑105 (MQ=255)
                                                                  acCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTAc                                                                                                                                <  2:527577/105‑1 (MQ=255)
                                                                    cATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGcgacga                                                                                                       <  1:77199/128‑1 (MQ=255)
                                                                    cATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGcgacga                                                                                                       >  2:77199/1‑128 (MQ=255)
                                                                    cATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCt                                                                                 >  2:332922/1‑150 (MQ=255)
                                                                     aTCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTc                                                                                >  1:502946/1‑150 (MQ=255)
                                                                     aTCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTc                                                                                >  1:1456009/1‑150 (MQ=255)
                                                                       cACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCAt                                                                              <  2:502946/150‑1 (MQ=255)
                                                                           ggAATGGCCAGCCATAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACATCTGTGGCGCTGTCGAGAAGGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGAcc                                                                          <  2:1087579/150‑1 (MQ=255)
                                                                            gAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGc                                                                                                                                            >  1:1175488/1‑83 (MQ=255)
                                                                            gAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGc                                                                                                                                            <  1:1241638/83‑1 (MQ=255)
                                                                            gAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGc                                                                                                                                            >  2:1241638/1‑83 (MQ=255)
                                                                            gAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGc                                                                                                                                            <  2:1175488/83‑1 (MQ=255)
                                                                             gatctccagccaCAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCtt                                                                        <  1:1720037/145‑1 (MQ=255)
                                                                             aaTGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCtt                                                                        <  2:930880/150‑1 (MQ=255)
                                                                               tGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGa                                                                      <  1:332922/150‑1 (MQ=255)
                                                                                 gccagccaCAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCaga                                                                                                              <  2:184055/108‑1 (MQ=255)
                                                                                 gccagccaCAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCaga                                                                                                              >  1:184055/1‑108 (MQ=255)
                                                                                  ccagccaCAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGAcc                                                                          >  2:1103840/1‑143 (MQ=255)
                                                                                  ccagccaCAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGAcc                                                                          <  1:1103840/143‑1 (MQ=255)
                                                                                   cagccaCAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGt                                                                  >  2:1937838/1‑150 (MQ=255)
                                                                                    agccaCAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTg                                                                 <  2:382678/150‑1 (MQ=255)
                                                                                    agccaCAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTg                                                                 <  1:1924641/150‑1 (MQ=255)
                                                                                      ccaCAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTAGCCCATGACCGCATTCCGCAGTTGCGCCAGGTACAGTCCCAGCCCCGCCAGAGCGACGATGACCGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGc                                                               <  1:2350928/150‑1 (MQ=255)
                                                                                       cacaGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACATTGACCGTGGcc                                                              <  2:1899910/150‑1 (MQ=255)
                                                                                        acaGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCt                                                             >  1:366714/1‑150 (MQ=255)
                                                                                                  cAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCa                                                   <  2:601188/150‑1 (MQ=255)
                                                                                                    ggAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCggcca                                                 <  2:1772639/150‑5 (MQ=255)
                                                                                                          ccgaggccgGCGGCACCTCTGTGGCGCTGTCGAGAATTCCACCGGTCCAGTGCGGCACCAGGTACAGACCGATCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGa                                           <  1:2222342/150‑1 (MQ=255)
                                                                                                           cgaggccgGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGAt                                          >  2:1702729/1‑150 (MQ=255)
                                                                                                           cgaggccgGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGAt                                          >  2:545093/1‑150 (MQ=255)
                                                                                                            gaggccgGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGAt                                                                                                      <  1:1342290/89‑1 (MQ=255)
                                                                                                            gaggccgGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGAt                                                                                                      >  2:1342290/1‑89 (MQ=255)
                                                                                                               gccgGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGcc                                      <  2:1666711/150‑1 (MQ=255)
                                                                                                                 cggcggCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCcag                                    <  2:1677304/150‑1 (MQ=255)
                                                                                                                   gcggcACCTCTGTGGCGCTGTCGAGAATGTCACCAGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGc                                                                                  <  1:2314649/102‑1 (MQ=255)
                                                                                                                   gcggcACCTCTGTGGCGCTGTCGAGAATGTCACCAGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGc                                                                                  >  2:2314649/1‑102 (MQ=255)
                                                                                                                   gcggcACCTCTGTGGCGCTGTCGAGAATGTCACCAGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGc                                                                                  >  2:2315920/1‑102 (MQ=255)
                                                                                                                   gcggcACCTCTGTGGCGCTGTCGAGAATGTCACCAGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGc                                                                                  <  1:2315920/102‑1 (MQ=255)
                                                                                                                        aCCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCg                             >  1:2237915/1‑150 (MQ=255)
                                                                                                                        aCCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCg                             >  2:757946/1‑150 (MQ=255)
                                                                                                                        aCCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCg                             >  1:2390517/1‑150 (MQ=255)
                                                                                                                          ctctGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACcagc                                                                                      >  1:1943785/1‑91 (MQ=255)
                                                                                                                          ctctGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACcagc                                                                                      <  2:1943785/91‑1 (MQ=255)
                                                                                                                          ctctGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAgc                                                 >  2:2183885/1‑128 (MQ=255)
                                                                                                                          ctctGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAgc                                                 <  1:2183885/128‑1 (MQ=255)
                                                                                                                          ctctGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCGAg                           <  1:1952689/150‑1 (MQ=255)
                                                                                                                          ctctGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCGAg                           >  1:1060215/1‑150 (MQ=255)
                                                                                                                                  cgcTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGc                                                               <  2:1748289/106‑1 (MQ=255)
                                                                                                                                  cgcTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGc                                                               >  1:1748289/1‑106 (MQ=255)
                                                                                                                                         gagaATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCGAGGATCAGCACGAATGa            <  2:2237915/150‑1 (MQ=255)
                                                                                                                                         gagaATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCGAGGATCAGCACGAATGa            >  2:2172393/1‑150 (MQ=255)
                                                                                                                                                ccACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCGAGGATCAGCACGAATGACAGGAtt     >  1:2376904/1‑150 (MQ=255)
                                                                                                                                                 cACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCGAGGATCAGCACGaa               <  1:464277/139‑1 (MQ=255)
                                                                                                                                                 cACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCGAGGATCAGCACGaa               >  2:464277/1‑139 (MQ=255)
                                                                                                                                                   ccGGTCCAGTGAGGCACCAGGTACAGGCCCAGAAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCGAGGATCAGCACGAATGACAGGATTGcc  <  1:1694751/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCGCCGTAACGGCGCTTCTGGTCGACGGCAAAGGCCGAGATGATCGGCGAATGGTTGAACGAGAACACCATCACCGGAATGGCCAGCCACAGGGTGTGCAGGAAGGCCGAGGCCGGCGGCACCTCTGTGGCGCTGTCGAGAATGCCACCGGTCCAGTGCGGCACCAGGTACAGGCCCAGCAACGCCAGAGCGACGATGAACGGGTAGACCAGCAGGCTCATGACCTTGACCGTGGCCTGCTCCCCGCAGCGCACGATGGCCAGCAGACCGAGGATCAGCACGAATGACAGGATTGCC  >  NC_002947/4078516‑4078812

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: ‑