Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,799,889 C→T 5.6% intergenic (+17/‑170) rnhB → / → dnaEA ribonuclease HII/DNA polymerase III subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,799,8890CT 5.6% 211.2 / 2.3 71intergenic (+17/‑170)rnhB/dnaEAribonuclease HII/DNA polymerase III subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (33/33);  new base T (2/2);  total (35/36)
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.19e-01

AAGGGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAGGCTACAATCCCGCCCTTGTCGT  >  NC_002947/1799741‑1800035
                                                                                                                                                    |                                                                                                                                                  
aaGGGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTcg                                                                                                                                                   >  2:1121138/1‑150 (MQ=255)
   ggCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTGGGCCTACGTCCATTCATCGGCGTTCCTTCGCCCCGGTGCTGGCGGCTTGGGAAGCGTGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGG                                                                                                                                                <  2:88314‑M1/150‑1 (MQ=255)
   ggCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGG                                                                                                                                                >  1:1511641‑M1/1‑150 (MQ=255)
     cTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCC                                                                                                                                              <  2:610948‑M1/150‑1 (MQ=255)
              cccGGTGCATCTCGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTAACAGGGCGCTGTTGTGGGGGCCGCTTCGCGGCCCTTCGCGGG                                                                                                                                     >  2:151382‑M1/1‑150 (MQ=255)
               ccGGTGCATCTCGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGC                                                                                                                                    >  2:619881‑M1/1‑150 (MQ=255)
                      aTCTCGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCG                                                                                                                             >  1:1220687‑M1/1‑150 (MQ=255)
                       tctcGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTcgc                                                                                                                                                  <  1:1165755/128‑1 (MQ=255)
                       tctcGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTcgc                                                                                                                                                  >  2:1165755/1‑128 (MQ=255)
                       tctcGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGC                                                                                                                            >  2:1861934‑M1/1‑150 (MQ=255)
                            aGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTca                                                                                                                       <  1:1121138/150‑1 (MQ=255)
                                 ttGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGa                                                                                                                  <  1:1870879/150‑1 (MQ=255)
                                  tGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCACAGTCACGCCCGCTCTCACAGGAt                                                                                                                 <  1:997778/150‑1 (MQ=255)
                                      aCGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCAt                                                                                                             >  1:1369852/1‑150 (MQ=255)
                                        gTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCg                                                                                                           >  1:1559687/1‑150 (MQ=255)
                                        gTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGCGTGTGGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCACGCTCTCACAGGATTCATCg                                                                                                           <  1:962514/150‑1 (MQ=39)
                                               gCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCAAGCCCGCTCTCACAGGATTCATCGGTCTGAg                                                                                                    <  1:991208/150‑1 (MQ=255)
                                                 cTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGTCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAgcg                                                                                                  <  2:1806386/150‑1 (MQ=255)
                                                     gCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGcggtgg                                                                                              <  1:1510181/150‑1 (MQ=255)
                                                      cccATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCggg                                                                                             <  1:95605/150‑1 (MQ=255)
                                                        cATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGc                                                                                           <  1:866751/150‑1 (MQ=255)
                                                             aTCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTAc                                                                                      >  1:643220/1‑150 (MQ=255)
                                                                ggCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCtg                                                                                   >  1:812196/1‑150 (MQ=255)
                                                                 gCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCtgt                                                                                  <  2:1369852/150‑1 (MQ=255)
                                                                  cGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCtgtg                                                                                 <  1:619881/150‑1 (MQ=255)
                                                                   gTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACGGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTgg                                                                                <  2:800305/150‑1 (MQ=255)
                                                                   gTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTgg                                                                                <  1:553802/150‑1 (MQ=255)
                                                                      ccTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTCTTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGgag                                                                             >  2:270901‑M2/1‑147 (MQ=255)
                                                                        ttCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCccacaggtacggccccgccgctcagaccgatgaatcctgtgagagcg                                                                           <  1:4923‑M1/150‑48 (MQ=255)
                                                                        ttCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGgagcg                                                                           <  1:1305267‑M2/150‑6 (MQ=255)
                                                                         tCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGgagcgg                                                                          >  1:1407170‑M2/1‑144 (MQ=255)
                                                                             cccGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTgg                                                                                >  2:210057/1‑140 (MQ=255)
                                                                             cccGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTgg                                                                                <  1:210057/140‑1 (MQ=255)
                                                                              ccGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGgagcgggcgtg                                                                     >  2:1559687‑M2/1‑139 (MQ=255)
                                                                                       gcggcTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGa                                                                                                                  >  2:1289563/1‑96 (MQ=255)
                                                                                       gcggcTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGa                                                                                                                  <  1:1289563/96‑1 (MQ=255)
                                                                                         ggcTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGgagcgggcgtgcccgcgaaggg                                                          >  1:756534‑M2/1‑128 (MQ=255)
                                                                                         ggcTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGgagcgggcgtgcccgcgaaggg                                                          >  2:1220687‑M2/1‑128 (MQ=255)
                                                                                           cTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCt                                                        <  2:643220/150‑1 (MQ=255)
                                                                                               ggAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCccacaggtacggccccgccgctcagaccgatgaatcctgt                                                                                  >  2:1546932‑M1/1‑80 (MQ=255)
                                                                                               ggAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCccacaggtacggccccgccgctcagaccgatgaatcctgt                                                                                  <  1:1546932‑M1/120‑41 (MQ=255)
                                                                                                   gcgcgcGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAagca                                                <  1:210077/150‑1 (MQ=255)
                                                                                                     gcgcgAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAagcagc                                              <  2:812196/150‑1 (MQ=255)
                                                                                                     gcgcgAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAagcagc                                              <  2:1305267/150‑1 (MQ=255)
                                                                                                       gcgAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCccacaggtacggccccgccgctcagaccgatgaatcctgtgagagcgggcgtgcccgcgaa                                                             >  2:1373615‑M1/1‑72 (MQ=255)
                                                                                                       gcgAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCccacaggtacggccccgccgctcagaccgatgaatcctgtgagagcgggcgtgcccgcgaa                                                             <  1:1373615‑M1/133‑62 (MQ=255)
                                                                                                               gTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCaaaaaaa                                    <  2:1511641/150‑1 (MQ=255)
                                                                                                                tCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAg                                   >  1:476653/1‑150 (MQ=255)
                                                                                                                tCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTAAAAAGGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAg                                   >  1:1366295/1‑150 (MQ=39)
                                                                                                                 cACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCAGCCCTTCGCGGGCACGCCCGCTCccacaggtacggccccgccgctcagaccgatgaatcctgtgagagcgggc                                                                        >  2:609530‑M1/1‑62 (MQ=255)
                                                                                                                 cACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCAGACCTTCGCGGGCACGCCCGCTCccacaggtacggccccgccgctcagaccgatgaatcctgtgagagcgggc                                                                        <  1:609530‑M1/112‑51 (MQ=255)
                                                                                                                    cGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGAt                                                                                                                 >  2:1291556/1‑68 (MQ=255)
                                                                                                                    cGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGAt                                                                                                                 <  1:1291556/68‑1 (MQ=255)
                                                                                                                     gATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGcccc                               <  2:163634/149‑1 (MQ=255)
                                                                                                                     gATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGcccc                               >  1:163634/1‑149 (MQ=255)
                                                                                                                      aTTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCGGCAAAGCAGCCCCAAAAAAAGCCCCTa                             >  1:995/1‑150 (MQ=255)
                                                                                                                        tCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGGACCTGTGGGAGCGGGCGGGCCCGCGAAGGGCTGCAAAGCAGCCCCCAAAAAAGCCCCTAgg                           >  1:1051740/1‑149 (MQ=255)
                                                                                                                          ccTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTa                         >  2:1822107/1‑150 (MQ=255)
                                                                                                                           cTGATCTGACCTATTGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCccacaggtacggccccgccgctcagaccgatgaatcctgtgagagc                                                                            >  2:328925‑M1/1‑52 (MQ=255)
                                                                                                                           cTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGcc                                             >  1:924251/1‑129 (MQ=255)
                                                                                                                           cTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGcc                                             <  2:924251/129‑1 (MQ=255)
                                                                                                                           cTGATCTGACCTATTGGGCCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCccacaggtacggccccgccgctcagaccgatgaatcctgtgagagc                                                                            <  1:328925‑M1/98‑47 (MQ=255)
                                                                                                                            tggtCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCCGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAgg                       >  2:1425062/4‑150 (MQ=255)
                                                                                                                                  gACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGCTTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAGGCTACaa                 >  2:1661114/1‑150 (MQ=255)
                                                                                                                                      tATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAGGCTACAATccc             >  2:908728/1‑150 (MQ=255)
                                                                                                                                        ttGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGcccc                               >  1:679867/1‑130 (MQ=255)
                                                                                                                                        ttGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGcccc                               <  2:679867/130‑1 (MQ=255)
                                                                                                                                        gaagggCCGCGAAGCGGCCCCTCGCGGGCACGCCCGCTCTCCCAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAagc                                                 <  2:1808683/109‑1 (MQ=37)
                                                                                                                                          ggggCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAGGCTACAATCCCGccc         <  1:1861934/150‑1 (MQ=255)
                                                                                                                                              ccGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAGGCTACAATCCCGCCCTTGt     >  2:669618/1‑150 (MQ=255)
                                                                                                                                               cGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAGGCTACAATCCCGCCCTTGTc    <  2:995/150‑1 (MQ=255)
                                                                                                                                                gCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAGGCTACAATCCCGCCCTTGTCg   >  1:1242595/1‑150 (MQ=255)
                                                                                                                                                gCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAGGCTACAATCCCGCCCTTGTCg   >  2:572565/1‑150 (MQ=255)
                                                                                                                                                 cTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAGGCTACAATCCCGCCCTTGTCGt  <  1:572565/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAGGGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAGGCTACAATCCCGCCCTTGTCGT  >  NC_002947/1799741‑1800035

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