Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,896,098 C→A 6.7% R216L (CGG→CTG)  PP_3439 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,896,0980CA 6.7% 227.1 / 5.0 74R216L (CGG→CTG) PP_3439AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base A (3/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

AGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGGCGCGTAT  >  NC_002947/3895950‑3896246
                                                                                                                                                    |                                                                                                                                                     
aGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGccc                                                                                                                                                       >  8:77864/1‑149 (MQ=255)
aGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGccc                                                                                                                                                       >  3:79989/1‑149 (MQ=255)
    gcgcCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGc                                                                                                                                                   <  4:159720/149‑1 (MQ=255)
      gcCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCa                                                                                                                                                <  7:33167/149‑1 (MQ=255)
       cccAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAg                                                                                                                                               >  5:100946/1‑149 (MQ=255)
           ggCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGc                                                                                                                                           <  5:99593/149‑1 (MQ=255)
                   ggggCATTACAGAAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTg                                                                                                                                   >  6:89396/3‑149 (MQ=255)
                    cGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTgg                                                                                                                                  <  7:77864/149‑1 (MQ=255)
                     ggCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGt                                                                                                                                 >  5:256636/1‑149 (MQ=255)
                      gCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTc                                                                                                                                <  5:354023/149‑1 (MQ=255)
                      gCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTc                                                                                                                                <  3:340690/149‑1 (MQ=255)
                       cATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTc                                                                                                                                >  8:338172/1‑148 (MQ=255)
                           tCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGc                                                                                                                            <  6:100946/148‑1 (MQ=255)
                            cTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCa                                                                                                                           <  3:158634/148‑1 (MQ=255)
                             tGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAgg                                                                                                                         >  2:193067/1‑149 (MQ=255)
                                    ccccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCAt                                                                                                                  <  4:362574/149‑1 (MQ=255)
                                     cccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGAGGTAGCCATgc                                                                                                                >  4:291139/1‑149 (MQ=255)
                                      ccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgc                                                                                                                >  7:214574/1‑149 (MQ=255)
                                         gtgtCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgcggc                                                                                                             <  4:79989/149‑1 (MQ=255)
                                             cgccgcCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGt                                                                                                           >  4:28324/1‑147 (MQ=255)
                                              gccgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCa                                                                                                         <  4:330505/148‑1 (MQ=255)
                                               ccgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAgg                                                                                                       >  7:29392/1‑149 (MQ=255)
                                               ccgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAgg                                                                                                       >  8:12438/1‑149 (MQ=255)
                                                     ggTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAg                                                                                                 <  8:170754/149‑1 (MQ=255)
                                                       tGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCg                                                                                               >  7:95082/1‑149 (MQ=255)
                                                        gCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGa                                                                                                   <  1:148280/144‑1 (MQ=255)
                                                        gCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGt                                                                                              >  6:249726/1‑149 (MQ=255)
                                                          tGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTg                                                                                            >  4:298120/1‑149 (MQ=255)
                                                            tttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTgg                                                                                           <  5:89396/148‑1 (MQ=255)
                                                            tttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGc                                                                                          >  5:89554/1‑149 (MQ=255)
                                                            tttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGc                                                                                          >  5:94812/1‑149 (MQ=255)
                                                             ttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGc                                                                                          >  1:18075/1‑148 (MQ=255)
                                                             ttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGc                                                                                          >  1:241489/1‑148 (MQ=255)
                                                             ttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGc                                                                                          >  1:23880/1‑148 (MQ=255)
                                                                 aTGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGaaa                                                                                      <  7:317733/148‑1 (MQ=255)
                                                                  tGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAgg                                                                                    <  1:374319/149‑1 (MQ=255)
                                                                  tGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAgg                                                                                    <  3:28324/149‑1 (MQ=255)
                                                                     ccGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCAc                                                                                 <  8:95082/149‑1 (MQ=255)
                                                                           cAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAAt                                                                             >  3:170995/1‑147 (MQ=255)
                                                                             ggTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGt                                                                        <  3:291139/149‑1 (MQ=255)
                                                                              gtgtGCCTGATCAGCGAACCCGGCCTCCTTTGCAACTTAAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGt                                                                        <  8:104709/149‑1 (MQ=255)
                                                                              gtgtGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGt                                                                        <  5:249726/149‑1 (MQ=255)
                                                                                   ccTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCTCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAc                                                                   <  3:203854/149‑1 (MQ=255)
                                                                                   ccTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTgg                                                                                           >  8:255556/1‑125 (MQ=255)
                                                                                   ccTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTgg                                                                                           <  7:255556/125‑1 (MQ=255)
                                                                                         cAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCgaga                                                              >  2:270742/1‑148 (MQ=255)
                                                                                          aGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGaaa                                                                                      <  4:342058/123‑1 (MQ=255)
                                                                                          aGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTACGACCGAGACa                                                            >  3:2918/1‑149 (MQ=255)
                                                                                                       tcctGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGcccc                                               >  3:369461/1‑149 (MQ=255)
                                                                                                        cctGTGGAACATCAGAAATTGAGCGGCCCTGGCGCAACTGCGACCGGGACAGTTTCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGATGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCg                                              >  2:142844/1‑149 (MQ=255)
                                                                                                        cctGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCg                                              >  4:212367/1‑149 (MQ=255)
                                                                                                          tgtgCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCg                                              >  5:229913/1‑147 (MQ=255)
                                                                                                            tgCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGTGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCa                                            <  1:270742/147‑1 (MQ=255)
                                                                                                            tgCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCa                                            <  8:214574/147‑1 (MQ=255)
                                                                                                                  ttCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCa                                                                                                                   <  2:338425/70‑1 (MQ=255)
                                                                                                                  ttCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCa                                                                                                                   >  1:338425/1‑70 (MQ=255)
                                                                                                                     aGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAg                                                                      <  2:37626/112‑1 (MQ=255)
                                                                                                                     aGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAg                                                                      >  1:37626/1‑112 (MQ=255)
                                                                                                                     aGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGa                                 <  7:338172/149‑1 (MQ=255)
                                                                                                                     aGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGa                                 <  4:369461/149‑1 (MQ=255)
                                                                                                                                 gCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGt                     <  8:29392/149‑1 (MQ=255)
                                                                                                                                 gCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGt                     <  2:23880/149‑1 (MQ=255)
                                                                                                                                 gCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTGGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGt                     <  2:18075/149‑1 (MQ=255)
                                                                                                                                  cccGGGCGCAGTTGCGCCAGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGAAAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGt                     >  6:86928/1‑148 (MQ=255)
                                                                                                                                   ccgggCGCAGTTGCGCCAGGGCCGGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCa                   >  6:326004/4‑149 (MQ=255)
                                                                                                                                   ccgggCGCAGTTGCGCCAGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTa                                                                       >  6:262141/4‑97 (MQ=255)
                                                                                                                                   ccgggCGCAGTTGCGCCAGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTa                                                                       <  5:262141/94‑1 (MQ=255)
                                                                                                                                       gtgcAACTGGCCCGGGCGCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGTGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGc               >  8:19128/3‑149 (MQ=255)
                                                                                                                                       gcgcAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACt                 <  4:2918/147‑1 (MQ=255)
                                                                                                                                       gcgcAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACt                 <  4:170995/147‑1 (MQ=255)
                                                                                                                                       gcgcAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGc               >  2:376446/1‑149 (MQ=255)
                                                                                                                                        cgcAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCg              <  7:12438/149‑1 (MQ=255)
                                                                                                                                          cagttGCGCCAGGGCCGGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGAt            <  5:326004/145‑1 (MQ=255)
                                                                                                                                              tGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGaaa                                                                                      >  2:180785/1‑71 (MQ=255)
                                                                                                                                              tGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGaaa                                                                                      <  1:180785/71‑1 (MQ=255)
                                                                                                                                                cgcCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGGATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGgcgc      <  1:193067/149‑1 (MQ=255)
                                                                                                                                                   ccGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGGCGCGTa   >  4:248503/1‑149 (MQ=255)
                                                                                                                                                    cGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGGCGCGTAt  <  2:241489/149‑1 (MQ=255)
                                                                                                                                                    cGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGGCGCGTAt  <  3:212367/149‑1 (MQ=255)
                                                                                                                                                    cGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGGAGGACCGAGACAGGCCGCAAGCCCCGGAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGGCGCGTAt  <  6:256636/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGGCGCGTAT  >  NC_002947/3895950‑3896246

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