Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 977,876 T→A 5.0% intergenic (+39/‑141) cysE → / → iscR serine acetyltransferase/DNA‑binding transcriptional regulator IscR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947977,8760TA 5.0% 2365.8 / 66.3 684intergenic (+39/‑141)cysE/iscRserine acetyltransferase/DNA‑binding transcriptional regulator IscR
Reads supporting (aligned to +/- strand):  ref base T (290/358);  new base A (16/19);  total (306/378)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 717 total aligned reads displayed.

GCGCTGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGCGACTG  >  NC_002947/977727‑978025
                                                                                                                                                       |                                                                                                                                                     
gcgcTGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCgtgt                                                                                                                                                       <  2:6379216/150‑1 (MQ=255)
 cgcTGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGtgtg                                                                                                                                                      <  2:2369374/150‑1 (MQ=255)
  gcTGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGc                                                                                                                                                     <  2:7955960/150‑1 (MQ=255)
        cAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGc                                                                                                                                               >  2:66829/1‑150 (MQ=255)
        cAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGc                                                                                                                                               >  1:5931427/1‑150 (MQ=255)
          aGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGAGCACACGccc                                                                                                                                             <  1:1637532/150‑1 (MQ=255)
          aGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGAGCACACGccc                                                                                                                                             <  2:3714916/150‑1 (MQ=255)
           gATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGAGCACACGcccc                                                                                                                                            <  1:5681360/150‑1 (MQ=255)
           gATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCCGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGAGCACACGcccc                                                                                                                                            <  1:7007643/150‑1 (MQ=255)
               ggCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGt                                                                                                                                        <  2:4148191/150‑1 (MQ=255)
                gCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGAGGGGCGTGAGCACACGccccgcatg                                                                                                                                       <  1:522334/150‑6 (MQ=255)
                 cAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCt                                                                                                                                      <  1:1892739/150‑1 (MQ=255)
                        tATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACaa                                                                                                                               >  1:6529703/1‑150 (MQ=255)
                        tATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACaa                                                                                                                               <  2:2506340/150‑1 (MQ=255)
                          ttGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATc                                                                                                                             <  1:2266308/150‑1 (MQ=255)
                          ttGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATc                                                                                                                             <  1:329826/150‑1 (MQ=255)
                          ttGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATc                                                                                                                             >  2:5677685/1‑150 (MQ=255)
                           tgtgCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATcg                                                                                                                            >  1:6053787/1‑150 (MQ=255)
                             tgCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATcgcg                                                                                                                          >  1:2012791/1‑150 (MQ=255)
                               ccAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGcgcc                                                                                                                        >  1:7279331/1‑150 (MQ=255)
                                 aaGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGcgcccg                                                                                                                      >  2:4934046/1‑150 (MQ=255)
                                  aGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGcgcccgc                                                                                                                     >  2:4819445/1‑150 (MQ=255)
                                     aaCTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTc                                                                                                                  <  2:7628564/150‑1 (MQ=255)
                                      aCTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTcc                                                                                                                 <  2:2017684/150‑1 (MQ=255)
                                        tGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCt                                                                                                               >  1:4199832/1‑150 (MQ=255)
                                           cGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGAt                                                                                                            <  1:1167262/150‑1 (MQ=255)
                                           cGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGAt                                                                                                            <  2:450815/150‑1 (MQ=255)
                                            ggCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTg                                                                                                              >  2:8414332/1‑147 (MQ=255)
                                            ggCCTTGCCGGAGGACGACTTCAGCGAGGCTGCACAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATg                                                                                                           <  1:8593988/150‑1 (MQ=255)
                                                  gCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATcc                                                                                                     <  2:3527831/150‑1 (MQ=255)
                                                   ccGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATccc                                                                                                    <  1:3915079/150‑1 (MQ=255)
                                                   ccGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATccc                                                                                                    <  2:3254041/150‑1 (MQ=255)
                                                     ggaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGa                                                                                                  <  1:4740818/150‑1 (MQ=255)
                                                     ggaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGa                                                                                                  <  2:6072822/150‑1 (MQ=255)
                                                     ggaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGa                                                                                                  >  2:4348918/1‑150 (MQ=255)
                                                      gaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGAc                                                                                                 >  2:737445/1‑150 (MQ=255)
                                                          acgacTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTaaa                                                                                             <  2:6887299/150‑1 (MQ=255)
                                                           cgacTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTTGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAg                                                                                            <  2:2271703/150‑1 (MQ=255)
                                                             acTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCa                                                                                          <  1:5576529/150‑1 (MQ=255)
                                                               ttCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGt                                                                                        >  2:2556929/1‑150 (MQ=255)
                                                               ttCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGt                                                                                        <  1:8057321/150‑1 (MQ=255)
                                                                   gCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGt                                                                                    <  2:5116751/150‑1 (MQ=255)
                                                                   gCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGAC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGt                                                                                    >  1:100801/1‑150 (MQ=255)
                                                                       ggCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCAGACGAGGGGCGTGAGCACACGccccgcat                                                                                                                                        >  2:1077498/1‑92 (MQ=25)
                                                                       ggCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTaa                                                                                <  1:4535824/150‑1 (MQ=255)
                                                                        gCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAAt                                                                               >  2:3060816/1‑150 (MQ=255)
                                                                        gCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAAt                                                                               >  2:541313/1‑150 (MQ=255)
                                                                          tGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACGATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAg                                                                             >  2:5992315/1‑150 (MQ=255)
                                                                            cccAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGtt                                                                           >  1:4465313/1‑150 (MQ=255)
                                                                              cAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGa                                                                         >  1:3805940/1‑150 (MQ=255)
                                                                              cAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGa                                                                         <  1:5036034/150‑1 (MQ=255)
                                                                               aaGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGAc                                                                        <  1:755422/150‑1 (MQ=255)
                                                                                    gCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACACTTGCGCCCGCCTCCCTGATGCAATCCCGACTCACGCAGTAGGTCTAAGAGGTGACTTaaa                                                                   >  1:7138852/1‑150 (MQ=255)
                                                                                    gCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTaaa                                                                   <  2:8205045/150‑1 (MQ=255)
                                                                                    gCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTaaa                                                                   <  1:2485448/150‑1 (MQ=255)
                                                                                    gCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTaaa                                                                   >  1:6754637/1‑150 (MQ=255)
                                                                                    gCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTaaa                                                                   <  1:8267824/150‑1 (MQ=255)
                                                                                     cccAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAAt                                                                  <  2:271242/150‑1 (MQ=255)
                                                                                      ccAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATg                                                                 >  2:4616988/1‑150 (MQ=255)
                                                                                      ccAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATg                                                                 >  1:3614995/1‑150 (MQ=255)
                                                                                               cacTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAAt                                                        <  2:1199112/150‑1 (MQ=255)
                                                                                                acTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATc                                                       <  2:1384950/150‑1 (MQ=255)
                                                                                                acTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATc                                                       >  1:995211/1‑150 (MQ=255)
                                                                                                acTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATc                                                       <  1:3487050/150‑1 (MQ=255)
                                                                                                   cAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCa                                                    <  2:7787648/150‑1 (MQ=255)
                                                                                                     gCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATa                                                  <  2:693587/150‑1 (MQ=255)
                                                                                                     gCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATa                                                  <  2:7464943/150‑1 (MQ=255)
                                                                                                        gCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGcccc                                                                                                                                            >  2:1583972/1‑57 (MQ=255)
                                                                                                        gCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGcccc                                                                                                                                            <  1:1583972/57‑1 (MQ=255)
                                                                                                           ttGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGca                                            <  2:8058568/150‑1 (MQ=255)
                                                                                                           ttGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGca                                            <  1:6418282/150‑1 (MQ=255)
                                                                                                           ttGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGca                                            <  2:4716787/150‑1 (MQ=255)
                                                                                                                cgcgCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACAt                                       >  1:4072724/1‑150 (MQ=255)
                                                                                                                 gcgcCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATc                                      <  2:2133846/150‑1 (MQ=255)
                                                                                                                  cgccgcACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATcc                                     >  1:2934989/1‑150 (MQ=255)
                                                                                                                       cACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTc                                                              >  2:8724269/1‑120 (MQ=255)
                                                                                                                       cACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTc                                                              <  1:8724269/120‑1 (MQ=255)
                                                                                                                       cACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGaaa                                <  2:28784/150‑1 (MQ=255)
                                                                                                                             aGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCg                          <  1:1103821/150‑1 (MQ=255)
                                                                                                                             aGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCg                          <  1:2191185/150‑1 (MQ=255)
                                                                                                                              gTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGt                         >  1:2650440/1‑150 (MQ=255)
                                                                                                                               tCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTATGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTg                        <  1:6256124/150‑1 (MQ=255)
                                                                                                                                 aTGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGt                      >  1:7809712/1‑150 (MQ=255)
                                                                                                                                   ggtggC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTAc                    <  1:3162224/150‑1 (MQ=255)
                                                                                                                                    gtggC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTAcc                   <  2:6310415/150‑1 (MQ=255)
                                                                                                                                       gC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCaa                >  2:1683697/1‑150 (MQ=255)
                                                                                                                                        c‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAAt               >  1:7948416/1‑150 (MQ=255)
                                                                                                                                            tGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAg             <  2:7063251/150‑1 (MQ=255)
                                                                                                                                             gaggggCGTGAGCACACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGc            >  2:8312453/3‑150 (MQ=255)
                                                                                                                                              cGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGcc           <  2:6563090/150‑1 (MQ=255)
                                                                                                                                               ggggCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCa          >  2:4026706/1‑150 (MQ=255)
                                                                                                                                               ggggCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCa          >  2:2661449/1‑150 (MQ=255)
                                                                                                                                                gggCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCAt         <  1:5420592/150‑1 (MQ=255)
                                                                                                                                                gggCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCAt         <  1:7399524/150‑1 (MQ=255)
                                                                                                                                                gggCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCAt         <  1:1769117/150‑1 (MQ=255)
                                                                                                                                                gggCGTGAGCACACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCAt         >  1:2750229/1‑150 (MQ=255)
                                                                                                                                                  gCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGc       <  2:7279331/150‑1 (MQ=255)
                                                                                                                                                  gCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGc       <  1:835951/150‑1 (MQ=255)
                                                                                                                                                    gtgtGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGCga     <  2:5861067/150‑1 (MQ=255)
                                                                                                                                                     tgtgCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGCgac    >  1:4804966/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GCGCTGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGC‑‑ATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGCGACTG  >  NC_002947/977727‑978025

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