Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,534,950 A→G 9.1% intergenic (+99/‑84) secA → / → argJ protein translocase subunit SecA/glutamate N‑acetyltransferase/amino acid acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,534,9500AG 9.1% 247.2 / 9.0 85intergenic (+99/‑84)secA/argJprotein translocase subunit SecA/glutamate N‑acetyltransferase/amino acid acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (39/38);  new base G (3/5);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTC  >  NC_002947/1534806‑1535090
                                                                                                                                                |                                                                                                                                            
gggTTCGGGCAAGAAGTTCAAGCATTGCAACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCtt                                                                                                                                         <  2:360159/149‑1 (MQ=255)
  gTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTtgt                                                                                                                                       <  1:992998/150‑1 (MQ=255)
  gTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTtgt                                                                                                                                       >  2:1618628/1‑150 (MQ=255)
       ggCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGGACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTggggg                                                                                                                                  >  1:966709/1‑150 (MQ=255)
        gCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGggc                                                                                                                                 <  2:1244477/150‑1 (MQ=255)
        gCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGggc                                                                                                                                 <  1:1317473/150‑1 (MQ=255)
          aagaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGggcgg                                                                                                                               >  2:895666/1‑150 (MQ=255)
           agaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGagcggc                                                                                                                              <  2:1044005/150‑1 (MQ=18)
            gaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGcc                                                                                                                             <  1:1223195/150‑1 (MQ=18)
             aagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGCCt                                                                                                                            <  2:1387264/150‑1 (MQ=18)
                ttCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTtgt                                                                                                                         >  1:1073520/1‑150 (MQ=255)
                ttCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGagcg                                                                                                                                <  2:411088/143‑1 (MQ=18)
                ttCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGagcg                                                                                                                                >  1:411088/1‑143 (MQ=17)
                 tCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTtgtg                                                                                                                        >  1:1383255/1‑150 (MQ=255)
                   aaGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTc                                                                                                                      <  1:1146671/150‑1 (MQ=255)
                   aaGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTc                                                                                                                      >  1:1351713/1‑150 (MQ=255)
                    aGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcg                                                                                                                     >  2:890326/1‑150 (MQ=255)
                                 tCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGAt                                                                                                                 >  1:1507289/1‑141 (MQ=255)
                                 tCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGAt                                                                                                                 <  2:1507289/141‑1 (MQ=255)
                                     aTCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGAGTCGCGATGGGGTGCGAAGCg                                                                                                    >  2:675559/1‑150 (MQ=255)
                                      tCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCATAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCgg                                                                                                   >  2:870932/1‑150 (MQ=255)
                                        gAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGcc                                                                                                 <  2:1617159/150‑1 (MQ=255)
                                           tGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCt                                                                                              <  2:966709/150‑1 (MQ=255)
                                              tCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCtttt                                                                                           >  2:1269451/1‑150 (MQ=255)
                                                tGAACCTTACGTTGTACAAGTCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTgg                                                                                                                                     >  1:1498579/1‑106 (MQ=18)
                                                tGAACCTTACGTTGTACAAGTCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTgg                                                                                                                                     <  2:1498579/106‑1 (MQ=18)
                                                tGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTc                                                                                         >  2:1538353/1‑150 (MQ=255)
                                                    ccTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGATGCGGCCCCTTTTTCGAca                                                                                     >  1:1067529/1‑150 (MQ=255)
                                                       tACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAcatca                                                                                  <  1:1618628/150‑1 (MQ=255)
                                                         cGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACatcatc                                                                                >  2:575776/1‑150 (MQ=255)
                                                              tACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttt                                                                           >  1:721968/1‑150 (MQ=255)
                                                                cAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttttt                                                                         >  1:1374811/1‑150 (MQ=255)
                                                                cAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttttt                                                                         <  2:1383255/150‑1 (MQ=255)
                                                                cAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttttt                                                                         >  1:959538/1‑150 (MQ=255)
                                                                   gCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTa                                                                      <  1:895666/150‑1 (MQ=255)
                                                                   gCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTa                                                                      <  2:1351713/150‑1 (MQ=255)
                                                                   gCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTa                                                                      <  2:414781/150‑1 (MQ=255)
                                                                    ccTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAg                                                                     <  2:959538/150‑1 (MQ=255)
                                                                    ccTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAg                                                                     <  1:890326/150‑1 (MQ=255)
                                                                          gTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGttt                                                                            >  2:994068/1‑137 (MQ=255)
                                                                          gTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGttt                                                                            <  1:994068/137‑1 (MQ=255)
                                                                           tCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATTAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCGTTTTAGTCATCATCGTTTTTTCTAGGAGCGct                                                              >  2:257100/1‑150 (MQ=255)
                                                                                ttGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATTAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAAt                                                         <  1:99061/150‑1 (MQ=255)
                                                                                   cGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGc                                                      <  2:1618488/150‑1 (MQ=255)
                                                                                     gCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAg                                                                     <  2:48926/133‑1 (MQ=255)
                                                                                     gCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAg                                                                     >  1:48926/1‑133 (MQ=255)
                                                                                     gCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTg                                                    <  1:1335613/150‑1 (MQ=255)
                                                                                      cccTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGt                                                   <  2:1380158/150‑1 (MQ=255)
                                                                                       ccTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGtt                                                  <  1:1618866/150‑1 (MQ=255)
                                                                                         tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTTGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttgg                                                >  2:262308/1‑150 (MQ=255)
                                                                                         tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttgg                                                >  2:815566/1‑150 (MQ=255)
                                                                                            cGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtct                                             <  1:1269451/150‑1 (MQ=255)
                                                                                             gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCGACGAAGATGGCGTCATCGTTGTTTCTAGGAGCGCTCTACTGGCTGttggtctt                                            >  2:1277082/1‑150 (MQ=255)
                                                                                             gCAACACAAGGCCGCTCCCACAGGGACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtctt                                            >  2:232946/1‑150 (MQ=255)
                                                                                               aaCACAAGGCCGCTTCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttgg                                          <  2:726741/150‑1 (MQ=255)
                                                                                                  acaAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttgg                                          <  1:1152117/147‑1 (MQ=255)
                                                                                                  acaAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttgg                                          >  2:1152117/1‑147 (MQ=255)
                                                                                                  acaAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCACCACTGGTAGGAGCGCTCTAATGGCTGttggtcttggtct                                       >  1:84490/1‑149 (MQ=255)
                                                                                                   caAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCt                                      >  1:952841/1‑150 (MQ=255)
                                                                                                    aaGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGc                                                                                                          <  2:1468439/81‑1 (MQ=255)
                                                                                                    aaGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGc                                                                                                          >  1:1468439/1‑81 (MQ=255)
                                                                                                    aaGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGccc                                                                                                <  2:1470853/91‑1 (MQ=255)
                                                                                                    aaGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGccc                                                                                                >  1:1470853/1‑91 (MQ=255)
                                                                                                    aaGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGa                                                                   <  2:648833/120‑1 (MQ=255)
                                                                                                    aaGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGa                                                                   >  1:648833/1‑120 (MQ=255)
                                                                                                     aGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCttt                                    >  1:961385/1‑150 (MQ=255)
                                                                                                     aGGCCGCTCCCACAGGTACCGCGATCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCttt                                    <  2:379563/150‑1 (MQ=255)
                                                                                                        ccGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGAACATTTGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACAACAGGGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGcc                                 >  1:1532548/1‑150 (MQ=255)
                                                                                                            tCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAcc                             <  1:675559/150‑1 (MQ=255)
                                                                                                             cccACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAccc                            <  2:944104/150‑1 (MQ=255)
                                                                                                                   ggTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTAACAGGCTTGGTCCTTTGCCCACCCTGCAcc                      <  1:1277082/150‑1 (MQ=255)
                                                                                                                   ggTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCCCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAcc                      >  1:670396/1‑150 (MQ=255)
                                                                                                                   ggTACCACGCTCGTCCATGAACGGGCTCCATCCTCGTGACGGCCTCCTCCTGTCGCGATCGCCTGCGACGCGGCCCCTTTTTCCACCTCATCCCTATGTCTAGGAAACCGCTAATGCCTGGTCTCGTTAGTCCTTTGTACAcccagaaca                      >  2:271075/1‑144 (MQ=255)
                                                                                                                       tcgcgCCCGTTCATGGACGAGCGCGGTGCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCgct                  >  2:529064/2‑148 (MQ=38)
                                                                                                                        cgcgCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTg                                                                                                           <  2:1548389/60‑1 (MQ=255)
                                                                                                                        cgcgCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTg                                                                                                           >  1:1548389/1‑60 (MQ=255)
                                                                                                                            cTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGCAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCgg             >  2:563740/1‑150 (MQ=255)
                                                                                                                                tCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcgc                                                                                                                    >  2:613927/1‑43 (MQ=255)
                                                                                                                                tCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcgc                                                                                                                    <  1:613927/43‑1 (MQ=255)
                                                                                                                                   aTGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGttt         >  1:960980/1‑147 (MQ=255)
                                                                                                                                   aTGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGttt         <  2:960980/147‑1 (MQ=255)
                                                                                                                                   aTGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGa      <  2:952841/150‑1 (MQ=255)
                                                                                                                                       aCGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTc  <  1:815566/150‑1 (MQ=255)
                                                                                                                                               gATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAcat                                                                                    >  1:1131564/1‑60 (MQ=255)
                                                                                                                                               gATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAcat                                                                                    <  2:1131564/60‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
TGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTC  >  NC_002947/1534806‑1535090

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: