Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,951,040 C→A 8.7% R658L (CGT→CTT)  ppk ← polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,951,0400CA 8.7% 335.6 / 12.2 110R658L (CGT→CTT) ppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (48/52);  new base A (6/4);  total (54/56)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

GGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAATGCCC  >  NC_002947/5950891‑5951189
                                                                                                                                                     |                                                                                                                                                     
ggTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGacac                                                                                                                                                       >  1:1032505/1‑150 (MQ=255)
 gTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACg                                                                                                                                                      >  2:1926682/1‑150 (MQ=255)
  ttCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGc                                                                                                                                                     <  1:1926682/150‑1 (MQ=255)
  ttCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGc                                                                                                                                                     >  1:1443171/1‑150 (MQ=255)
          gCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAg                                                                                                                                             <  2:526417/150‑1 (MQ=255)
           ccGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAgg                                                                                                                                            <  1:1867852/150‑1 (MQ=255)
           ccGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAgg                                                                                                                                            <  1:1734911/150‑1 (MQ=255)
            cGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGTTCTCGACAAGCGTGTCGAgac                                                                                                                                           <  1:2158303/150‑3 (MQ=255)
             ggTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGtt                                                                                                                                          >  1:1746228/1‑150 (MQ=255)
                cGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTgcg                                                                                                                                       <  1:16764/150‑1 (MQ=255)
                        tgcGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtccatcc                                                                                                                               >  2:693426/1‑150 (MQ=255)
                            cACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTc                                                                                                                           <  1:326487/150‑1 (MQ=255)
                                    gcCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGAga                                                                                                                                            <  2:56432/125‑2 (MQ=255)
                                    gcCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGAga                                                                                                                                            >  1:56432/1‑124 (MQ=255)
                                     cccGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGa                                                                                                                  >  2:256817/1‑150 (MQ=255)
                                      ccGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGAc                                                                                                                 >  1:1444426/1‑150 (MQ=255)
                                        gTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTgcgc                                                                                                                                      >  2:913725/1‑127 (MQ=255)
                                        gTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTgcgc                                                                                                                                      <  1:913725/127‑1 (MQ=255)
                                        gTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAg                                                                                                               >  1:1432142/1‑150 (MQ=255)
                                          cTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGt                                                                                                             <  1:256817/150‑1 (MQ=255)
                                            ggCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGCCTCGACAAGCGTGTCGAgacc                                                                                                                                          >  2:1381106/1‑116 (MQ=255)
                                            ggCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGCCTCGACAAGCGTGTCGAgacc                                                                                                                                          <  1:1381106/119‑4 (MQ=255)
                                               tGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATc                                                                                                        >  1:2175761/1‑150 (MQ=255)
                                               tGCAAGGTCCAGGCGTGGGTATTGTAGGCCAGGTAACCTTCCAGCTCCTTCTTCACCCGCAGCGACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCTGGTCGAGGTTGCGCTCCACCCAGTCGCCGCCGGACAGGGAGATc                                                                                                        >  1:564456/1‑150 (MQ=255)
                                                     gTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGctct                                                                                                  >  1:1353290/1‑150 (MQ=255)
                                                      tCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCtt                                                                                                 >  1:2014754/1‑150 (MQ=255)
                                                       ccAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTc                                                                                                <  2:240350/150‑1 (MQ=255)
                                                       ccAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTc                                                                                                <  2:1104353/150‑1 (MQ=255)
                                                        cAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGAt                                                                                                         <  2:1581243/140‑1 (MQ=255)
                                                        cAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGAt                                                                                                         >  1:1581243/1‑140 (MQ=255)
                                                          ggCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTcgcc                                                                                             <  2:2239713/150‑1 (MQ=255)
                                                             gTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGa                                                                                                                                              >  1:803126/1‑98 (MQ=255)
                                                             gTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGtt                                                                                                                                          <  2:1767054/102‑1 (MQ=255)
                                                             gTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCgccgcc                                                                                          >  2:1925582/1‑150 (MQ=255)
                                                                 gTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTg                                                                                      <  1:2202154/150‑1 (MQ=255)
                                                                    ttGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAgg                                                                                   <  2:1032505/150‑1 (MQ=255)
                                                                      gccGGTCGTGTAGCTTTCCAGCTCCTTCTTCACACGCAGCAACTGTTTCCTGACTTCCGTCGTGAGGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGTCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGac                                                                                 <  2:564456/148‑2 (MQ=255)
                                                                      gTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGaa                                                                                 <  1:1636427/150‑1 (MQ=255)
                                                                         ggTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTa                                                                              >  2:444472/1‑150 (MQ=255)
                                                                         ggTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTa                                                                              >  1:864440/1‑150 (MQ=255)
                                                                          gTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCa                                                                                         >  1:20474/1‑138 (MQ=255)
                                                                              ggTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTGGTCGAGGTTGCGCTCAATCCAGTCGGCGCTAGTCATTTAGTTCTGCTCTTCGCCGCCATTGAGGAAGTAGAaca                                                                         >  1:1345063/1‑150 (MQ=255)
                                                                               gTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGCCAGTTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAacac                                                                        >  1:2231378/1‑150 (MQ=255)
                                                                                   ccTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACAccc                                                                      >  2:836908/1‑148 (MQ=255)
                                                                                   ccTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACAccc                                                                      <  1:836908/148‑1 (MQ=255)
                                                                                   ccTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCcgc                                                                    <  2:730799/150‑1 (MQ=255)
                                                                                   ccTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCcgc                                                                    <  2:2175761/150‑1 (MQ=255)
                                                                                    cTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGc                                                                                                                                                     <  1:1077665/68‑1 (MQ=255)
                                                                                    cTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGc                                                                                                                                                     >  2:1077665/1‑68 (MQ=255)
                                                                                     ttCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCgt                                                                  <  2:824707/150‑1 (MQ=255)
                                                                                         aGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGc                                                              <  2:864440/150‑1 (MQ=255)
                                                                                          gCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAg                                                                                                                                             <  1:1432895/70‑1 (MQ=255)
                                                                                          gCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAg                                                                                                                                             >  2:1432895/1‑70 (MQ=255)
                                                                                          gCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGCGTGcc                                                             >  2:366123/1‑149 (MQ=255)
                                                                                             ccTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTg                                                                                                                                         <  1:1080215/71‑1 (MQ=255)
                                                                                             ccTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTg                                                                                                                                         >  2:1080215/1‑71 (MQ=255)
                                                                                              cttcttCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAg                                                         >  2:1102777/1‑150 (MQ=255)
                                                                                               ttcttcACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAgg                                                        >  2:2184278/1‑150 (MQ=255)
                                                                                                  ttcACCCGCAGCAACAGTTTCTTGCCTTACACCGTGAATCTGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGTCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAg                                                     <  2:1378740/150‑1 (MQ=255)
                                                                                                   tcACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGc                                                    >  1:1405756/1‑150 (MQ=255)
                                                                                                    cACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCg                                                   <  2:714478/150‑1 (MQ=255)
                                                                                                     aCCCGCAGCAACAGTTTCTGGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCgg                                                  >  1:566116/1‑150 (MQ=255)
                                                                                                      cccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGc                                                 >  2:398198/1‑150 (MQ=255)
                                                                                                       ccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGAACAGGTAGATCTGCTCTTCGCCGCCATTGAg                                                                                    >  1:202915/1‑114 (MQ=255)
                                                                                                       ccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGAACAGGTAGATCTGCTCTTCGCCGCCATTGAg                                                                                    <  2:202915/114‑1 (MQ=255)
                                                                                                       ccGCAGCAACAGTTTCTTGCCTTCCACCAGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCATCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGcc                                                >  2:1256713/1‑150 (MQ=255)
                                                                                                          cagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAAt                                             <  1:1096527/150‑1 (MQ=255)
                                                                                                          cagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAAt                                             >  1:2286909/1‑150 (MQ=255)
                                                                                                            gcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGa                                           <  1:1102777/150‑1 (MQ=255)
                                                                                                              aaCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGaa                                                                           >  1:379531/1‑116 (MQ=255)
                                                                                                              aaCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGaa                                                                           <  2:379531/116‑1 (MQ=255)
                                                                                                              aaCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATc                                         >  2:696374/1‑150 (MQ=255)
                                                                                                              aaCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATc                                         <  1:444472/150‑1 (MQ=255)
                                                                                                               aCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGa                                                                            >  1:819211/1‑114 (MQ=255)
                                                                                                               aCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGa                                                                            <  2:819211/114‑1 (MQ=255)
                                                                                                               aCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCg                                        >  1:2318140/1‑150 (MQ=255)
                                                                                                                 aGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACTCGCTTGTCGAGGTTGCGCTCCTTCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTTAGGAATTAGAACACCCGCGTGGGCTCTAGGAAGCGGCCAATGATCGAg                                      <  2:492913/150‑1 (MQ=255)
                                                                                                                  gTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGa                                                          <  1:1988040/129‑1 (MQ=255)
                                                                                                                  gTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGa                                                          >  2:1988040/1‑129 (MQ=255)
                                                                                                                       ttGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCGCGTGTTAGCGGCCAATGATCGGGCCCAcc                                >  2:1136962/1‑150 (MQ=255)
                                                                                                                        tGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCt                               <  2:2286909/150‑1 (MQ=255)
                                                                                                                         gCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAgg                                                                                                                                            <  1:1674148/40‑1 (MQ=255)
                                                                                                                         gCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAgg                                                                                                                                            >  2:1674148/1‑40 (MQ=255)
                                                                                                                                aCCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtgttg                       >  2:1315236/1‑150 (MQ=255)
                                                                                                                                 ccGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtgttgt                      <  1:1925582/150‑1 (MQ=255)
                                                                                                                                  cGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCAcc                                >  2:920470/1‑139 (MQ=255)
                                                                                                                                  cGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCAcc                                <  1:920470/139‑1 (MQ=255)
                                                                                                                                  cGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTtgtg                     <  2:1432142/150‑1 (MQ=255)
                                                                                                                                   gggAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGc                    <  1:1453136/150‑1 (MQ=255)
                                                                                                                                     gAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGa                  >  2:929393/1‑150 (MQ=255)
                                                                                                                                          aGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGaa                                                                                 <  1:2136939/82‑1 (MQ=255)
                                                                                                                                          aGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGaa                                                                                 >  2:2136939/1‑82 (MQ=255)
                                                                                                                                             tctcGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCg          <  2:1443171/150‑1 (MQ=255)
                                                                                                                                             cctcGACAAGCGTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGc                                                              <  2:1376286/97‑1 (MQ=255)
                                                                                                                                             cctcGACAAGCGTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAAcacacg          >  2:426560/2‑145 (MQ=255)
                                                                                                                                             cctcGACAAGCGTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGCAGATCTGCTCTTCGCCGCCAGTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCg          >  2:1326029/2‑150 (MQ=255)
                                                                                                                                               tcGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGa        >  2:2325006/1‑150 (MQ=255)
                                                                                                                                                  acacGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAAccc             <  2:1631629/142‑1 (MQ=255)
                                                                                                                                                  acacGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAAccc             >  1:1631629/1‑142 (MQ=255)
                                                                                                                                                     cGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAATGccc  <  1:1256713/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAATGCCC  >  NC_002947/5950891‑5951189

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: