Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,799,892 G→A 8.3% intergenic (+20/‑167) 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,8920GA 8.3% 200.2 / 5.2 71intergenic (+20/‑167)rnhB/dnaEAribonuclease HII/DNA polymerase III subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (32/33);  new base A (3/3);  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 = 1.00e+00

GGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAGGCTACAATCCCGCCCTTGTCGT  >  NC_002947/1799744‑1800035
                                                                                                                                                    |                                                                                                                                               
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)
                    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                                                                                                                  <  1:1289563/96‑1 (MQ=255)
                                                                                    gcggcTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGa                                                                                                                  >  2:1289563/1‑96 (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                                                                                  <  1:1546932‑M1/120‑41 (MQ=255)
                                                                                            ggAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCccacaggtacggccccgccgctcagaccgatgaatcctgt                                                                                  >  2:1546932‑M1/1‑80 (MQ=255)
                                                                                                gcgcgcGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAagca                                                <  1:210077/150‑1 (MQ=255)
                                                                                                  gcgcgAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAagcagc                                              <  2:1305267/150‑1 (MQ=255)
                                                                                                  gcgcgAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAagcagc                                              <  2:812196/150‑1 (MQ=255)
                                                                                                    gcgAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCccacaggtacggccccgccgctcagaccgatgaatcctgtgagagcgggcgtgcccgcgaa                                                             <  1:1373615‑M1/133‑62 (MQ=255)
                                                                                                    gcgAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCccacaggtacggccccgccgctcagaccgatgaatcctgtgagagcgggcgtgcccgcgaa                                                             >  2:1373615‑M1/1‑72 (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                               >  1:163634/1‑149 (MQ=255)
                                                                                                                  gATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGcccc                               <  2:163634/149‑1 (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                                             <  2:924251/129‑1 (MQ=255)
                                                                                                                        cTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGcc                                             >  1:924251/1‑129 (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)
                                                                                                                                                    |                                                                                                                                               
GGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTGGGCCTACGCCCATTCATCGGCGTTCCTTCGCCCCGGTGCGGGCGGCTTGGGAAGCGCGCGAAGGCGTCACCGATTCCCTGATCTGACCTATTGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTCTCACAGGATTCATCGGTCTGAGCGGCGGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAAAAAAGCCCCTAGTTAGGCTACAATCCCGCCCTTGTCGT  >  NC_002947/1799744‑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: