Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,184,221 C→A 5.2% P93P (CCC→CCA PP_5490 → site‑specific DNA recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,184,2210CA 5.2% 252.2 / 2.4 77P93P (CCC→CCAPP_5490site‑specific DNA recombinase
Reads supporting (aligned to +/- strand):  ref base C (40/33);  new base A (2/2);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

AGTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCAAAG  >  NC_002947/2184073‑2184365
                                                                                                                                                     |                                                                                                                                                 
agTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTTGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGccc                                                                                                                                                   >  3:373435/1‑149 (MQ=255)
 gTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCt                                                                                                                                                  <  8:318699/149‑1 (MQ=255)
    gatgaGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGc                                                                                                                                               >  7:320952/1‑149 (MQ=255)
     atgaGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCAGGGCt                                                                                                                                              >  6:365214/1‑149 (MQ=255)
           ctgctTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCAGGGCTCagaca                                                                                                                                        <  5:365214/149‑5 (MQ=255)
             gctTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAAcacca                                                                                                                                      >  6:369252/1‑149 (MQ=255)
              ctTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAAcacca                                                                                                                                      <  6:378744/148‑1 (MQ=255)
                 gATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTAc                                                                                                                                  >  5:387445/1‑149 (MQ=255)
                   tcatcTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCAGGGCTCaga                                                                                                                                          >  4:332450/1‑137 (MQ=255)
                   tcatcTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCAGGGCTCaga                                                                                                                                          <  3:332450/139‑3 (MQ=255)
                       ctcAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGgctgc                                                                                                                            >  6:328922/1‑149 (MQ=255)
                        tcAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGgctgct                                                                                                                           >  1:222584/1‑149 (MQ=255)
                        tcAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGgctgct                                                                                                                           >  3:97839/1‑149 (MQ=255)
                          aggccaggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAac                                                                                                                                          <  3:263373/132‑1 (MQ=255)
                          aggccaggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAac                                                                                                                                          >  4:263373/1‑132 (MQ=255)
                          aggccaggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGgctgct                                                                                                                           <  2:222584/147‑1 (MQ=255)
                          aggccaggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTgg                                                                                                                         >  3:302209/1‑149 (MQ=255)
                             ccaggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTc                                                                                                                       <  1:336154/148‑1 (MQ=255)
                               aggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCACTGGCTCAACACCACTACGGCTGCTGGt                                                                                                                        >  1:105121/1‑146 (MQ=255)
                                ggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAgat                                                                                                                   <  6:387445/149‑1 (MQ=255)
                                 gcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAgat                                                                                                                   >  5:68315/1‑148 (MQ=255)
                                     tgtggttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATaa                                                                                                              >  5:389368/1‑149 (MQ=255)
                                       tggttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATg                                                                                                            <  6:389368/149‑1 (MQ=255)
                                         gttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGa                                                                                                           <  5:328922/148‑1 (MQ=255)
                                               gtgAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGGTGAATCGCTTGGGCGCTGGGTTTCAGTCTCTGTCTGAGCCCTGGCTCaa                                                                                                                                           <  3:124569/110‑1 (MQ=255)
                                                  aGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTg                                                                                                  >  5:152901/1‑148 (MQ=255)
                                                       tcTTGATCGGCTCGCACGGTCGACAATTGATTTTCTGTCTATTGCAGAGCAGCTGTATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTtgcatgca                                                                                            <  8:320952/149‑1 (MQ=255)
                                                       tcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAacac                                                                                                                                        >  4:208525/1‑105 (MQ=255)
                                                       tcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAacac                                                                                                                                        <  3:208525/105‑1 (MQ=255)
                                                       tcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCa                                                                                            >  8:335608/1‑149 (MQ=255)
                                                        cTTGATCGTCGCTCACGGTCGACAATTCATTTGCTGTCTAATGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCAt                                                                                           <  4:43124/149‑1 (MQ=255)
                                                           gATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGc                                                                                        >  6:362043/1‑149 (MQ=255)
                                                                        ggTCGACAATTGATTTGCTGTCTATTGGAGAGCAGCTGAATCGCATGCGCGCGGGGTACCAGTCTCTGTCTGAGCACTGGCTCAACACAACTACCGGAGCAGGTCAGAGGATAAAGAGGAAATTTGAAGGCATCGATAGCTGCGAggg                                                                            >  6:102675/1‑146 (MQ=255)
                                                                               aaTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGa                                                                   <  6:152901/149‑1 (MQ=255)
                                                                               aaTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGa                                                                   >  1:156023/1‑149 (MQ=255)
                                                                               aaTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGa                                                                   >  1:294679/1‑149 (MQ=255)
                                                                                aTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGAt                                                                  <  1:122330/149‑1 (MQ=255)
                                                                                aTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGAt                                                                  >  2:122330/1‑149 (MQ=255)
                                                                                 ttGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATc                                                                 <  4:302209/149‑1 (MQ=255)
                                                                                     aTTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCg                                                                >  3:170668/1‑147 (MQ=255)
                                                                                           tGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCg                                                         >  7:77766/1‑148 (MQ=255)
                                                                                               tATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATc                                                                 <  7:132016/136‑1 (MQ=255)
                                                                                               tATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATc                                                                 >  8:132016/1‑136 (MQ=255)
                                                                                                     agagCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCgg                                              >  7:187376/1‑149 (MQ=255)
                                                                                                      gagCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCggg                                             <  2:294679/149‑1 (MQ=255)
                                                                                                        gcagcTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTc                                           >  2:73411/1‑149 (MQ=255)
                                                                                                           gcTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGTGCGGGTCgag                                        >  7:95857/1‑149 (MQ=255)
                                                                                                           gcTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCgag                                        <  4:373435/149‑1 (MQ=255)
                                                                                                            cTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTa                                                              >  4:349610/1‑126 (MQ=255)
                                                                                                            cTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTa                                                              <  3:349610/126‑1 (MQ=255)
                                                                                                                 tCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAg                                  >  7:32254/1‑149 (MQ=255)
                                                                                                                   gCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAgcg                                >  6:231459/1‑149 (MQ=255)
                                                                                                                   gCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAgc                                 <  6:68315/148‑1 (MQ=255)
                                                                                                                     ttGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAgcg                                >  2:306672/1‑147 (MQ=255)
                                                                                                                            cTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAg                       <  8:77766/149‑1 (MQ=255)
                                                                                                                            cTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAg                       <  2:55078/149‑1 (MQ=255)
                                                                                                                            cTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAg                       <  2:156023/149‑1 (MQ=255)
                                                                                                                             tGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTCATCGTAGACCGGACGAGAGCGCGACGAGAGGCAGCGCAGAGGAAgg                      >  4:245300/1‑149 (MQ=255)
                                                                                                                               ggTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGc                                                >  8:253505/1‑121 (MQ=255)
                                                                                                                               ggTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAggg                     >  2:187324/1‑148 (MQ=255)
                                                                                                                               ggTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAgg                      <  2:105121/147‑1 (MQ=255)
                                                                                                                               ggTTGCAGTCTCTGTCTGAGCCCTGGCTCAACAACACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGATCGCAGACCTGATCGTAGACCGGACGAGAGc                                                <  7:253505/122‑1 (MQ=255)
                                                                                                                                    cAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTa                                                              >  2:308630/1‑102 (MQ=255)
                                                                                                                                    cAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTa                                                              <  1:308630/102‑1 (MQ=255)
                                                                                                                                     aGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATaaaa              <  8:187376/149‑1 (MQ=255)
                                                                                                                                       tctcTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGTGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAAtt            <  5:362043/149‑1 (MQ=255)
                                                                                                                                       tctcTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAAtt            <  1:306672/149‑1 (MQ=255)
                                                                                                                                         tctgtctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAgg                      >  1:191429/1‑137 (MQ=255)
                                                                                                                                         tctgtctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAgg                      <  2:191429/137‑1 (MQ=255)
                                                                                                                                         tctgtctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCg          >  3:258622/1‑149 (MQ=255)
                                                                                                                                          ctgtctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATAATTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGACGGGCATAAAATACgg         >  4:306878/1‑149 (MQ=255)
                                                                                                                                           tgtctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTc           <  8:32254/146‑1 (MQ=255)
                                                                                                                                             tctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTc       <  1:73411/148‑1 (MQ=255)
                                                                                                                                               tgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGGGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCa     >  1:176951/1‑148 (MQ=255)
                                                                                                                                                gAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGt        <  1:91348/144‑1 (MQ=255)
                                                                                                                                                gAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGt        >  2:91348/1‑144 (MQ=255)
                                                                                                                                                gAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCa     >  7:183269/1‑147 (MQ=255)
                                                                                                                                                 aGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCGTCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCAAAg  >  2:234318/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AGTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCAAAG  >  NC_002947/2184073‑2184365

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: