Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,534,940 A→G 22.5% intergenic (+89/‑94) 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,9400AG22.5% 293.0 / 77.4 142intergenic (+89/‑94)secA/argJprotein translocase subunit SecA/glutamate N‑acetyltransferase/amino acid acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (49/61);  new base G (13/19);  total (62/80)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

Only 100 of 144 total aligned reads displayed.

GAGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAAC  >  NC_002947/1534792‑1535088
                                                                                                                                                    |                                                                                                                                                     
gAGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGaa                                                                                                                                                      <  2:4997316/150‑1 (MQ=255)
 aGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAAc                                                                                                                                                     <  1:3067375/150‑1 (MQ=255)
 aGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAAc                                                                                                                                                     <  1:874318/150‑1 (MQ=255)
 aGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAAc                                                                                                                                                     <  1:193772/150‑1 (MQ=255)
 aGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAAc                                                                                                                                                     <  1:4225523/150‑1 (MQ=255)
  gCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACg                                                                                                                                                    <  2:2915500/150‑1 (MQ=255)
     gTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACCGGGATCGCGCCCGTTCATGGACGAgc                                                                                                                                                 <  2:1466990/150‑1 (MQ=38)
        cTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGcgcg                                                                                                                                              <  1:653115/150‑1 (MQ=255)
           gtgtGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTATCGCGCCCGTTCATGGACGAGcgcggta                                                                                                                                           <  1:1889216/150‑4 (MQ=38)
            tgtgGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATcc                                                                                                                                          <  2:4137767/150‑1 (MQ=255)
            tgtgGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAAGGATCGCGCCCGTTCATGGACGAGcgcggtac                                                                                                                                          <  1:5037912/150‑5 (MQ=37)
             gtgGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGAtacc                                                                                                                                         <  2:4704183/150‑4 (MQ=255)
              tgGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCtt                                                                                                                                        <  2:2275019/150‑1 (MQ=255)
                ggtCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTCGCAACCCATCGCGACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCtgtg                                                                                                                                     <  1:884867/148‑1 (MQ=14)
                gTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCATCGCGACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCtgtg                                                                                                                                     <  2:706523/150‑1 (MQ=18)
                 ttCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCtgtg                                                                                                                                     <  1:2763224/150‑1 (MQ=25)
                 ttCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCtgtg                                                                                                                                     <  1:1074060/150‑1 (MQ=25)
                  tCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTgg                                                                                                                                    <  1:807290/150‑1 (MQ=37)
                   cGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCACACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTggg                                                                                                                                   <  2:1607391/150‑1 (MQ=255)
                    gggCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTgggg                                                                                                                                  <  1:2661767/150‑1 (MQ=255)
                    gggCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTggga                                                                                                                                  <  2:2330666/150‑2 (MQ=21)
                      gCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGggc                                                                                                                                >  1:2396758/1‑150 (MQ=255)
                      gCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCATTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGggc                                                                                                                                >  2:2076185/1‑150 (MQ=255)
                          gaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGcc                                                                                                                            >  2:352474/1‑150 (MQ=255)
                           aagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCt                                                                                                                           <  2:4580103/150‑1 (MQ=255)
                            agTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGCCtt                                                                                                                          <  1:383537/150‑1 (MQ=18)
                            agTTCAAGCATTGCCACGGTCAGATCGAGTGATCAGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTACGCAACACAAGGCCGCTCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGCCtt                                                                                                                          <  2:2025386/150‑1 (MQ=18)
                              ttCAAGCATGGCCACGGTCAGATCGAGTGATCTGAACCGTACGTTGGACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTtgt                                                                                                                        <  2:1164830/150‑1 (MQ=255)
                                 aaGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTc                                                                                                                     >  1:949369/1‑150 (MQ=255)
                                 aaGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTc                                                                                                                     >  2:2525377/1‑150 (MQ=255)
                                 aaGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTc                                                                                                                     >  1:954959/1‑150 (MQ=255)
                                 aaGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTc                                                                                                                     >  1:3527183/1‑150 (MQ=255)
                                      ttGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTTGCGaa                                                                                                               <  2:4953736/150‑2 (MQ=255)
                                       tGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTTGCGaaa                                                                                                              <  2:4562583/150‑3 (MQ=255)
                                            cGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTg                                                                                                         <  2:1041938/150‑1 (MQ=255)
                                                cAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCGTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGaa                                                                                                     >  2:4014122/1‑150 (MQ=255)
                                                   aTCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCg                                                                                                  <  1:3766038/150‑1 (MQ=255)
                                                    tCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCGTTGTGGGGGCGGGCTTGTGTCGCGATGGGTGGGAAGCGGc                                                                                                >  1:5061260/1‑150 (MQ=255)
                                                     cGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGc                                                                                                <  1:4672147/150‑1 (MQ=255)
                                                       aGTGATCTGAACCTGACGTTGTACCCGCCTGGGGTCGCTGTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGccc                                                                                              <  2:4449596/150‑1 (MQ=255)
                                                         tGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCt                                                                                            >  1:1615991/1‑150 (MQ=255)
                                                             cTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCttttt                                                                                        <  1:3482786/150‑1 (MQ=255)
                                                              tGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTc                                                                                       <  2:199451/150‑1 (MQ=255)
                                                              tGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTACCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTc                                                                                       >  1:4898867/1‑150 (MQ=255)
                                                               gAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCg                                                                                      >  1:1216695/1‑150 (MQ=255)
                                                                  ccTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAca                                                                                   >  1:1640256/1‑150 (MQ=255)
                                                                       cGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACatcatc                                                                              <  1:3827186/150‑1 (MQ=255)
                                                                           gTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGttt                                                                          <  1:1198605/150‑1 (MQ=255)
                                                                            tACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCATCACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttt                                                                         <  1:1569096/150‑1 (MQ=255)
                                                                            tACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttt                                                                         >  1:4506916/1‑150 (MQ=255)
                                                                            tACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttt                                                                         >  1:1182419/1‑150 (MQ=255)
                                                                               aaGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTc                                                                      >  1:1322876/1‑150 (MQ=255)
                                                                                  ccTGGGGTCTCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAg                                                                   <  2:1119898/150‑1 (MQ=255)
                                                                                  ccTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAg                                                                   <  1:3913401/150‑1 (MQ=255)
                                                                                    tGGGGTCGCTTTGCGGCCCTTTCGCAAAACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAgta                                                                 >  1:2188994/1‑148 (MQ=255)
                                                                                           gCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGctct                                                          <  2:4254419/150‑1 (MQ=255)
                                                                                             tttGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTaa                                                        >  1:1225138/1‑150 (MQ=255)
                                                                                                    cccTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGt                                                 <  2:4415337/150‑1 (MQ=255)
                                                                                                     ccTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGtt                                                >  1:26173/1‑150 (MQ=255)
                                                                                                       tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttgg                                              <  2:73776/150‑1 (MQ=255)
                                                                                                           gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCCTCGTTTTGGCTAGGAGCGCTCTAATGGCTGatgggctt                                          >  2:1863417/1‑150 (MQ=255)
                                                                                                            cAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttg                                         >  2:878362/1‑150 (MQ=255)
                                                                                                                acaAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACCTCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTcc                                     >  2:4992362/1‑150 (MQ=255)
                                                                                                                acaAGGCCGCCCCCACAAGGACCGCGCTCGTCCATGAACGGGCGCGa                                                                                                                                             <  2:4306709/47‑1 (MQ=17)
                                                                                                                acaAGGCCGCCCCCACAAGGACCGCGCTCGTCCATGAACGGGCGCGa                                                                                                                                             >  1:4306709/1‑47 (MQ=17)
                                                                                                                  aaGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCtt                                   <  1:1171891/150‑1 (MQ=255)
                                                                                                                    ggCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTg                                 <  2:4875547/150‑1 (MQ=255)
                                                                                                                     gCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGc                                <  2:318050/150‑1 (MQ=255)
                                                                                                                       cgcgCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGccc                              <  2:3367509/150‑1 (MQ=255)
                                                                                                                       cGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGccc                              <  2:3367364/150‑1 (MQ=255)
                                                                                                                        gcccccACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCa                             >  1:3673271/4‑150 (MQ=21)
                                                                                                                        gCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttgg                                        <  2:3082811/139‑1 (MQ=255)
                                                                                                                         cccccACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTTGCGAAAGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCa                             >  1:942085/3‑150 (MQ=17)
                                                                                                                         cccccACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAc                            >  1:672699/3‑150 (MQ=21)
                                                                                                                         cccccACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCGGTGGGAGCGGCCTTGTGTAGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAc                            >  1:3237153/3‑150 (MQ=21)
                                                                                                                          tCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAcc                           <  2:1632116/150‑1 (MQ=255)
                                                                                                                          ccccACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGGGGGAGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAcc                           >  1:4171497/2‑150 (MQ=21)
                                                                                                                            ccACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCt                         >  2:105221/1‑150 (MQ=255)
                                                                                                                             cacaGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTg                        >  1:392090/1‑150 (MQ=255)
                                                                                                                              acaGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGc                       >  2:2476433/1‑150 (MQ=255)
                                                                                                                               caGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCa                      <  1:3553225/150‑1 (MQ=255)
                                                                                                                                aaggatcgcgCCCGTTCATGGACGAGCGCGGTACCTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAc                     >  1:2603389/7‑150 (MQ=37)
                                                                                                                                aaggatcgcgCCCGTTCATGGACGAGCGCGGTACCTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAc                     >  1:50510/7‑150 (MQ=37)
                                                                                                                                aaggatcgcgCCCGTTCATGGACGAGCGCGGTACCTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAc                     >  1:2603585/7‑150 (MQ=37)
                                                                                                                                    atcgcgCCCGTTCATGGACGAGCGCGGTACCTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCgg                 >  2:1361026/3‑150 (MQ=37)
                                                                                                                                    atcgcgCCCGTTCATGGACGAGCGCGGTACCTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCgg                 >  2:1361088/3‑150 (MQ=37)
                                                                                                                                     ccGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGt                >  1:1976167/1‑150 (MQ=255)
                                                                                                                                       gcgcTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTc              >  2:514264/1‑150 (MQ=255)
                                                                                                                                        cgcCCGTTCATGGACGAGCGCGGTACCTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTCGCCCACCCTGCACCCGGTTcc             >  2:3674410/1‑150 (MQ=38)
                                                                                                                                         gcTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCg            <  1:1015145/150‑1 (MQ=255)
                                                                                                                                          cTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCgg           <  2:1887813/150‑1 (MQ=255)
                                                                                                                                            cGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGt         >  1:4770924/1‑150 (MQ=255)
                                                                                                                                             gTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtt        <  2:2396758/150‑1 (MQ=255)
                                                                                                                                              tCCCTGACCGGGCGCGATCCTTGGGGGGGCGGCCGTGGGGCGCGAGGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTCTTTCTCGGCGCGCTCTCATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGttt       <  1:3075634/150‑1 (MQ=255)
                                                                                                                                              tCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGttt       <  2:4898867/150‑1 (MQ=255)
                                                                                                                                              tCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGttt       <  2:916821/150‑1 (MQ=255)
                                                                                                                                               ccATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtttt      <  1:352474/150‑1 (MQ=255)
                                                                                                                                                cATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtttt      >  1:796708/1‑149 (MQ=255)
                                                                                                                                                cATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtttt      <  2:796708/149‑1 (MQ=255)
                                                                                                                                                   gAACGGGCGCGATCCTGGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAAc  >  2:2074457/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAAC  >  NC_002947/1534792‑1535088

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