Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,344,106 G→T 5.5% P13T (CCT→ACT)  PP_1168 ← TRAP dicarboxylate transporter subunit DctQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,344,1060GT 5.5% 277.4 / 2.3 87P13T (CCT→ACT) PP_1168TRAP dicarboxylate transporter subunit DctQ
Reads supporting (aligned to +/- strand):  ref base G (41/41);  new base T (2/3);  total (43/44)
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.78e-01

CAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGG  >  NC_002947/1343963‑1344252
                                                                                                                                               |                                                                                                                                                  
cAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGgcggcg                                                                                                                                              >  1:680410/1‑150 (MQ=255)
cAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGgcggcg                                                                                                                                              >  2:400322/1‑150 (MQ=255)
     cGTAGCGGGCGAAGATGCCCCACGGAATGCTCACTGGCAGGTCGGCGGTCGACAGGCCCGCGATCCAGATGAGACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGgcgc                                                                                                                                         <  1:540270/150‑1 (MQ=255)
      gTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGcgcc                                                                                                                                        >  2:867540/1‑150 (MQ=255)
           gggCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCt                                                                                                                                   <  1:2147469/150‑1 (MQ=255)
             gCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGc                                                                                                                                 >  2:2351010/1‑150 (MQ=255)
              cGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGACCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCTGcc                                                                                                                                <  1:487642/150‑1 (MQ=255)
               gAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCg                                                                                                                               >  1:2371417/1‑150 (MQ=255)
                       cccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGggcg                                                                                                                                                 >  1:2385811/1‑124 (MQ=255)
                       cccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGggcg                                                                                                                                                 <  2:2385811/124‑1 (MQ=255)
                        ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGt                                                                                                                           <  2:165240/145‑1 (MQ=255)
                        ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGt                                                                                                                           >  1:165240/1‑145 (MQ=255)
                        ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAtt                                                                                                                      <  1:867540/150‑1 (MQ=255)
                        ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAtt                                                                                                                      >  2:657376/1‑150 (MQ=255)
                        ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAtt                                                                                                                      >  2:1450604/1‑150 (MQ=255)
                        ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAtt                                                                                                                      >  1:1180302/1‑150 (MQ=255)
                          aCGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTgg                                                                                                                    <  2:2371417/150‑1 (MQ=255)
                            gggATGATCAGCGTCATGGCCAGGATCGACAGGCCCGGGGTCGAGATGCAACTGGGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGc                                                                                                                  <  2:803661/150‑1 (MQ=255)
                            gggATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGc                                                                                                                  >  2:2224001/1‑150 (MQ=255)
                              gatgatCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCaa                                                                                                                >  2:785421/1‑150 (MQ=255)
                                 gatCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAggg                                                                                                             <  2:1440673/150‑1 (MQ=255)
                                   tCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcg                                                                                                           >  2:2142044/1‑150 (MQ=255)
                                   tCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcg                                                                                                           >  2:2013406/1‑150 (MQ=255)
                                      gCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCt                                                                                                        <  1:626909/150‑1 (MQ=255)
                                      gCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCt                                                                                                        >  1:652197/1‑150 (MQ=255)
                                      gCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCt                                                                                                        <  1:1515944/150‑1 (MQ=255)
                                         tCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGtt                                                                                                     >  2:470267/1‑150 (MQ=255)
                                          cATGGTCAGGCTCGACAGGCCCGCGATCCAGATTCAACTGATGTACAGCGCACCATGCTCGCCCCGCAAATGCGCTTTCTTACCCTTCACCACAGCCGCAGGGTGGCCTGCGCCTCGCCAGCGCCGCTGAGTGGGCAAGGGCGGgtgtaa                                                                                                    <  2:2031130/150‑3 (MQ=255)
                                          cATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTa                                                                                                    <  2:1794102/150‑1 (MQ=255)
                                              gCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTgg                                                                                                <  1:1448808/150‑1 (MQ=255)
                                                      cGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTc                                                                                        <  2:1140465/150‑1 (MQ=255)
                                                        aCAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGa                                                                                      >  2:2383258/1‑150 (MQ=255)
                                                         cAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGAt                                                                                     <  1:785421/150‑1 (MQ=255)
                                                          aGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATa                                                                                    >  1:1635709/1‑150 (MQ=255)
                                                              ccGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATa                                                                                    <  2:1154523/146‑1 (MQ=255)
                                                              ccGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATa                                                                                    >  1:1154523/1‑146 (MQ=255)
                                                                gcgATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTtc                                                                              >  2:1377539/1‑150 (MQ=255)
                                                                  gATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTtcat                                                                            >  2:1499228/1‑150 (MQ=255)
                                                                     ccAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAg                                                                         >  1:1446000/1‑150 (MQ=255)
                                                                      cAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGAcc                                                                                             <  2:1649467/129‑2 (MQ=255)
                                                                      cAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACg                                                                                             >  1:1649467/1‑129 (MQ=255)
                                                                      cAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAgg                                                                        <  2:652197/150‑1 (MQ=255)
                                                                       aGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTggg                                                                                                                   <  1:720662/106‑1 (MQ=255)
                                                                       aGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTggg                                                                                                                   >  2:720662/1‑106 (MQ=255)
                                                                       aGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGt                                                                       >  1:1787959/1‑150 (MQ=255)
                                                                        gATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTc                                                                      <  1:1450604/150‑1 (MQ=255)
                                                                                tGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACt                                                              <  1:470267/150‑1 (MQ=255)
                                                                                     tACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTAc                                                                                                   >  1:2210559/1‑108 (MQ=255)
                                                                                     tACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTAc                                                                                                   <  2:2210559/108‑1 (MQ=255)
                                                                                           gtgtCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGt                                                                       >  2:15588/1‑130 (MQ=255)
                                                                                           gtgtCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGt                                                                       <  1:15588/130‑1 (MQ=255)
                                                                                            tgtCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTAtt                                                  >  2:186139/1‑150 (MQ=255)
                                                                                                gTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                 <  2:919861/115‑1 (MQ=255)
                                                                                                gTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                 >  1:919861/1‑115 (MQ=255)
                                                                                                gTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAgg                                                                        >  2:1821325/1‑124 (MQ=255)
                                                                                                gTTCACGCGCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAgg                                                                        <  1:1821325/124‑1 (MQ=255)
                                                                                                 ttCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcg                                             >  1:1209012/1‑150 (MQ=255)
                                                                                                          aGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGaa                                    >  2:1208512/1‑150 (MQ=255)
                                                                                                               aaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGt                               <  1:1377539/150‑1 (MQ=255)
                                                                                                               aaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCAGATCCTGCACGCCGTAGGGATACTTCATTTATGCGTGGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGt                               <  1:657376/150‑1 (MQ=255)
                                                                                                                aGCTATTTCATAGGCATCACCAAAGCGGTAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAATTGGCAAGGGCGGCTGTTACTGTAGGGCGTCGATACGCGTCATCAGGTAGCCGTACTGCGCGTCATATTTTTAGCGTACCCAAGACGtt                              <  2:370694/150‑2 (MQ=255)
                                                                                                                 gCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTgg                             <  2:1446000/150‑1 (MQ=255)
                                                                                                                 gCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTgg                             <  1:400322/150‑1 (MQ=255)
                                                                                                                    aTTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcgcg                                           <  2:1257619/133‑1 (MQ=255)
                                                                                                                    aTTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcgcg                                           >  1:1257619/1‑133 (MQ=255)
                                                                                                                     tttCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAGGCCGTGGcggg                         >  1:36027/1‑148 (MQ=255)
                                                                                                                         aTAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTg                              >  2:350497/1‑141 (MQ=255)
                                                                                                                         aTAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTg                              <  1:350497/141‑1 (MQ=255)
                                                                                                                         aTAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAg                     >  2:1210653/1‑150 (MQ=255)
                                                                                                                            ggCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAAc                  <  1:186139/150‑1 (MQ=255)
                                                                                                                             gCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgtc                         >  1:444919/1‑142 (MQ=255)
                                                                                                                             gCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgtc                         <  2:444919/142‑1 (MQ=255)
                                                                                                                              cATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACgg                <  1:1897031/150‑1 (MQ=255)
                                                                                                                               aTCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGc               >  2:2407054/1‑150 (MQ=255)
                                                                                                                                  aCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCg                                                                                          <  1:2395716/72‑1 (MQ=255)
                                                                                                                                  aCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCg                                                                                          >  2:2395716/1‑72 (MQ=255)
                                                                                                                                  aCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttc            <  1:2013406/150‑1 (MQ=255)
                                                                                                                                    cAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCATCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTg                              >  2:2367619/1‑130 (MQ=255)
                                                                                                                                    cAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCATCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTg                              <  1:2367619/130‑1 (MQ=255)
                                                                                                                                     aaaGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttcttt         >  2:681625/1‑149 (MQ=255)
                                                                                                                                       aGCTGCAGGGCGGCGGTCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTc       <  1:33830/150‑1 (MQ=255)
                                                                                                                                        tcGGCAGTGCGGCGCCCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTTCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCg      >  1:637160/2‑150 (MQ=255)
                                                                                                                                        gCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGc                                  <  2:1052169/122‑1 (MQ=255)
                                                                                                                                        gCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGc                                  >  1:1052169/1‑122 (MQ=255)
                                                                                                                                         cGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTTTCGa     <  1:681625/150‑1 (MQ=255)
                                                                                                                                          ggCAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAAc                  <  2:1733905/136‑1 (MQ=255)
                                                                                                                                          ggCAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAAc                  >  1:1733905/1‑136 (MQ=255)
                                                                                                                                           gCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACg   <  2:1209012/150‑1 (MQ=255)
                                                                                                                                            cAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACgg  <  1:697084/150‑1 (MQ=255)
                                                                                                                                              gggCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTAc                                        <  1:2160025/110‑1 (MQ=255)
                                                                                                                                              gggCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTAc                                        >  2:2160025/1‑110 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGG  >  NC_002947/1343963‑1344252

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: