Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,534,952 C→A 8.8% intergenic (+101/‑82) 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,9520CA 8.8% 332.5 / 13.8 110intergenic (+101/‑82)secA/argJprotein translocase subunit SecA/glutamate N‑acetyltransferase/amino acid acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (51/49);  new base A (5/5);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 113 total aligned reads displayed.

GGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATCGCTTC  >  NC_002947/1534807‑1535101
                                                                                                                                                 |                                                                                                                                                      
ggTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCTAGTGATCTGAACCTTACGTTGTACAGGCACGCCGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTtg                                                                                                                                                    <  1:1853077/150‑1 (MQ=255)
 gTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTTGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTtgt                                                                                                                                                   >  2:237864/1‑150 (MQ=255)
 gTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTtgt                                                                                                                                                   <  1:1888335/150‑1 (MQ=255)
   tCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGTGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTgg                                                                                                                                                 <  2:1398286/150‑1 (MQ=255)
   tCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTgg                                                                                                                                                 <  2:51537/150‑1 (MQ=255)
   tCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTgg                                                                                                                                                 <  1:1616296/150‑1 (MQ=255)
         aagaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCTTCCATGAACGGGCGCGATCCTTGTGGGggcgg                                                                                                                                           >  1:12087/1‑150 (MQ=255)
           gaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGAGGGGGCGGcc                                                                                                                                         >  2:1357515/1‑150 (MQ=255)
             agTTCACGCATTGCCACTGTCAGATCGAGTGATCTGAACCTTACGTTGTCCAAGCCTGGGGTCGCTTTGCGGCCCTCTCGCAACACAAGGCCGCTCCCACAGGTACCGCGATCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCtt                                                                                                                                       <  1:317075/150‑1 (MQ=255)
               ttCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCATCGCGACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGagcg                                                                                                                                            <  2:1704495/142‑1 (MQ=9)
               ttCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCATCGCGACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGagcg                                                                                                                                            >  1:1704495/1‑142 (MQ=9)
                  aaGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTc                                                                                                                                  >  1:925075/1‑150 (MQ=255)
                           cACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATgggg                                                                                                                        >  1:387539/1‑150 (MQ=255)
                                    aTCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGAGGCCTTGTGTCGCGATGGGGTGCGAAGCg                                                                                                               >  2:2363516/1‑150 (MQ=255)
                                     tCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTtgt                                                                                                                                     <  2:2192005/128‑1 (MQ=255)
                                     tCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTtgt                                                                                                                                     >  1:2192005/1‑128 (MQ=255)
                                     tCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCgg                                                                                                              <  1:615935/150‑1 (MQ=255)
                                      cGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGc                                                                                                             <  1:1111804/150‑1 (MQ=255)
                                       gAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCGCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGagcggc                                                                                                                                          <  2:1878239/121‑1 (MQ=9)
                                            aTCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCttt                                                                                                       >  1:853735/1‑150 (MQ=255)
                                             tCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCACCACAAGGCGGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTTTCGCGATGGGGTGCGAATCGGccccaggc                                                                                                      >  2:251050/1‑146 (MQ=255)
                                              cTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCttttt                                                                                                     <  1:1812761/150‑1 (MQ=255)
                                               tGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATc                                                                                                                                                        <  2:2124608/99‑1 (MQ=255)
                                               tGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATc                                                                                                                                                        >  1:2124608/1‑99 (MQ=255)
                                                   ccTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAca                                                                                                <  2:329735/150‑1 (MQ=255)
                                                     ttACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGggcggc                                                                                                                                          >  1:1577146/1‑107 (MQ=255)
                                                     ttACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGggcggc                                                                                                                                          <  2:1577146/107‑1 (MQ=255)
                                                     ttACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAcatc                                                                                              >  1:2383377/1‑150 (MQ=255)
                                                         gTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCg                                                                                          >  2:836641/1‑150 (MQ=255)
                                                         gTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCg                                                                                          >  1:1193995/1‑150 (MQ=255)
                                                         gTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCg                                                                                          >  2:1896274/1‑150 (MQ=255)
                                                              aCAAGCCTGGGGTCGCTTTGCGGCCCTTTCGAAACAAAAAAAGGCTAAAAAAATTAAGGCGGTCGTCCATTTACTGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGttttt                                                                                     >  1:658113/1‑150 (MQ=255)
                                                                aaGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTgg                                                                                                                                                 >  1:1712203/1‑89 (MQ=255)
                                                                   ccTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAg                                                                                >  1:824564/1‑150 (MQ=255)
                                                                     tGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGa                                                                              <  2:853735/150‑1 (MQ=255)
                                                                            gCTTTGCGGCCCTTTCGCAACAAAAGGCCGCTCCCACGGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGctct                                                                       <  2:869420/150‑1 (MQ=255)
                                                                                tGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGGCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATg                                                                   <  1:1348020/150‑1 (MQ=255)
                                                                                   ggCCCTTTCGCAACACAAGGCCACTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCGCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCt                                                                >  2:1103001/1‑150 (MQ=255)
                                                                                        tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTc                                                                                                                                  >  2:834115/1‑80 (MQ=255)
                                                                                        tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTc                                                                                                                                  <  1:834115/80‑1 (MQ=255)
                                                                                        tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCCTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttgg                                                           >  1:1189251/1‑150 (MQ=255)
                                                                                            gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttt                                                                                      >  1:2204498/1‑119 (MQ=255)
                                                                                            gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttt                                                                                      <  2:2204498/119‑1 (MQ=255)
                                                                                            gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACAACATCGGTGTTTCTAGGAGCGCTCTAATGGCTGttggtctt                                                       >  2:1599422/1‑150 (MQ=255)
                                                                                            gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCGTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtctt                                                       >  1:531149/1‑150 (MQ=255)
                                                                                             cAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCAGGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttg                                                      >  2:1389100/1‑150 (MQ=255)
                                                                                              gacacAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttgg                                                     >  1:1877885/2‑150 (MQ=17)
                                                                                              aaCACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttgg                                                     >  1:583750/1‑150 (MQ=255)
                                                                                              aaCACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttgg                                                     <  2:1302466/150‑1 (MQ=255)
                                                                                              aaCACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttgg                                                     <  2:387539/150‑1 (MQ=255)
                                                                                                cacaAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAg                                                                                >  2:1315502/1‑121 (MQ=255)
                                                                                                cacaAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAg                                                                                <  1:1315502/121‑1 (MQ=255)
                                                                                                  caAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCt                                                 <  1:836641/150‑1 (MQ=255)
                                                                                                   aaGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCtt                                                >  2:1668124/1‑150 (MQ=255)
                                                                                                   aaGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCAGCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCtt                                                >  2:1739132/1‑150 (MQ=255)
                                                                                                   aaGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTTGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTaa                                                                     >  1:1732761/1‑129 (MQ=2)
                                                                                                   aaGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGCCTCGTGTTGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTaa                                                                     <  2:1732761/129‑1 (MQ=2)
                                                                                                    aGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGTGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCttt                                               >  2:247881/1‑150 (MQ=255)
                                                                                                    aGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCttt                                               <  1:237864/150‑1 (MQ=255)
                                                                                                    aGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCttt                                               <  1:1896274/150‑1 (MQ=255)
                                                                                                     ggCCGCTCCCAACGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTg                                              >  2:349592/1‑150 (MQ=255)
                                                                                                      gCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGc                                             >  1:981062/1‑150 (MQ=255)
                                                                                                        cGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGccc                                           <  2:583750/150‑1 (MQ=255)
                                                                                                          cTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAc                                         <  1:51370/150‑1 (MQ=255)
                                                                                                             ccACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGcc                                                                                                            <  1:1826855/80‑1 (MQ=255)
                                                                                                             ccACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGcc                                                                                                            >  2:1826855/1‑80 (MQ=255)
                                                                                                             ccACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCt                                      >  2:1187998/1‑150 (MQ=255)
                                                                                                             ccACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCt                                      <  1:699117/150‑1 (MQ=255)
                                                                                                               acaGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAcatca                                                                                             <  2:766470/93‑1 (MQ=255)
                                                                                                               acaGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAcatca                                                                                             >  1:766470/1‑93 (MQ=255)
                                                                                                                caGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCa                                   >  1:2005406/1‑150 (MQ=255)
                                                                                                                  ggTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAcc                                 <  1:1187998/150‑1 (MQ=255)
                                                                                                                   gTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAcc                                        <  1:1156521/142‑1 (MQ=255)
                                                                                                                   gTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAcc                                        >  2:1156521/1‑142 (MQ=255)
                                                                                                                   gTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGCGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTCCCCACCCTGCAccc                                >  2:1347823/1‑150 (MQ=255)
                                                                                                                    tACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTcc                                                  <  1:955640/131‑1 (MQ=255)
                                                                                                                      ccGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGt                             >  1:1973400/1‑150 (MQ=255)
                                                                                                                       cgcgCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGtt                            >  1:1634094/1‑150 (MQ=255)
                                                                                                                       cgcgCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGtt                            >  1:2191881/1‑150 (MQ=255)
                                                                                                                        gcgcTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTc                           <  2:1564431/150‑1 (MQ=255)
                                                                                                                        gcgcTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGCGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTATGCTTGGTCCTTTGCCCACCCTGCAGCAGGTTc                           >  1:2336996/1‑150 (MQ=39)
                                                                                                                          gcTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCggaaccg                         <  1:1389100/150‑6 (MQ=255)
                                                                                                                             cGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGt                      <  1:1599422/150‑1 (MQ=255)
                                                                                                                              gTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtt                     <  2:1634094/150‑1 (MQ=255)
                                                                                                                               ttcATGGACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTTGCGAAAGGGCCGCAAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtt                     <  1:1725401/148‑1 (MQ=18)
                                                                                                                               tCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCCTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGttt                    <  2:1189251/150‑1 (MQ=255)
                                                                                                                                 cATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTg                  <  1:1574807/150‑1 (MQ=255)
                                                                                                                                  aTGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCTTCTTTATCTTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTTCCCACCCTGCACCCGGTTCCGGGTTTTGa                 >  2:719138/1‑150 (MQ=255)
                                                                                                                                  aTGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTTCCCACCCTGCACCCGGTTCCGGGTTTTGa                 >  2:929183/1‑150 (MQ=255)
                                                                                                                                   tGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGaa                <  2:797569/150‑1 (MQ=255)
                                                                                                                                   tGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGaa                <  2:1973400/150‑1 (MQ=255)
                                                                                                                                      aCGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTc             <  2:2383377/150‑1 (MQ=255)
                                                                                                                                       cGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGCAGCGTTACCCTTGGCTAGGACGGCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCg            <  1:1489366/150‑1 (MQ=255)
                                                                                                                                       cGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCg            >  1:155309/1‑150 (MQ=255)
                                                                                                                                          gcgcGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTa         >  2:642587/1‑150 (MQ=255)
                                                                                                                                            gcgATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCgg           >  1:193418/1‑146 (MQ=255)
                                                                                                                                            gcgATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCgg           <  2:193418/146‑1 (MQ=255)
                                                                                                                                             cgATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATCg      <  1:722561/150‑1 (MQ=255)
                                                                                                                                                tCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTcc                                                  <  1:2000830/103‑1 (MQ=255)
                                                                                                                                                 ccTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTATGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATCGCTTc  >  2:424991/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
GGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATCGCTTC  >  NC_002947/1534807‑1535101

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: