Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,224,401 G→T 6.8% A49E (GCA→GAA)  PP_1962 ← site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,224,4010GT 6.8% 227.3 / 4.9 74A49E (GCA→GAA) PP_1962site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (34/34);  new base T (2/3);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

ACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCCGGCGGTCTCCTGCAT  >  NC_002947/2224258‑2224546
                                                                                                                                               |                                                                                                                                                 
aCCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGTGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGCATGTAACATTCGCtt                                                                                                                                              <  7:359565/149‑1 (MQ=255)
   gATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGCATGTAACATACGCTttga                                                                                                                                           <  6:259910/149‑2 (MQ=255)
       gTCGTGAGTTTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCTCCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTttgtgtt                                                                                                                                        <  1:208499/148‑1 (MQ=255)
       gTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTttgtgttg                                                                                                                                       <  1:347789/149‑1 (MQ=255)
         cGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTttgtgttgt                                                                                                                                      <  8:114773/148‑1 (MQ=255)
          gTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTtgtgttgtgt                                                                                                                                    >  6:366248/1‑149 (MQ=255)
          gTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTtgtgttgtgt                                                                                                                                    >  6:140281/1‑149 (MQ=255)
               aTTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTttgtgtt                                                                                                                                        <  5:341888/140‑1 (MQ=255)
               aTTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTttgtgtt                                                                                                                                        >  6:341888/1‑140 (MQ=255)
               aTTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGATTTGTGTTGTGTCCGc                                                                                                                                >  5:49770/1‑148 (MQ=255)
                           atgatACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATg                                                                                                                   >  5:353679/1‑149 (MQ=255)
                               tACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGt                                                                                                               >  3:27352/1‑149 (MQ=255)
                                  gAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAAt                                                                                                                    <  7:10524/141‑1 (MQ=255)
                                  gAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAAt                                                                                                                    >  8:10524/1‑141 (MQ=255)
                                  gAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGtctc                                                                                                            <  3:63409/149‑1 (MQ=255)
                                      gAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGaa                                                                                                        >  1:276986/1‑149 (MQ=255)
                                            aTAAGATTGAAGGGAGTAAATGCAGTTAGCAAATGTTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCg                                                                                                   <  5:151270/148‑1 (MQ=255)
                                                gATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTttgtgtt                                                                                                                                        <  2:372973/107‑1 (MQ=255)
                                                gATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTttgtgtt                                                                                                                                        >  1:372973/1‑107 (MQ=255)
                                                  ttGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTaa                                                                                            >  8:276504/1‑149 (MQ=255)
                                                    gAAGGGAGTAAATGCAGTGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAAt                                                                                           <  5:140281/148‑1 (MQ=255)
                                                    gAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAAt                                                                                           <  7:276504/148‑1 (MQ=255)
                                                               aTGCAGGGAGCAAATGGTAAACAATTGCCATTTTCAGGTCGTTCCCCCCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGaaa                                                                               >  1:239536/1‑149 (MQ=255)
                                                                 gCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTc                                                                                                     <  2:330288/125‑1 (MQ=255)
                                                                 gCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTc                                                                                                     >  1:330288/1‑125 (MQ=255)
                                                                         cAAATGGCAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTg                                                                     <  7:140308/149‑1 (MQ=255)
                                                                          aaaTGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTgg                                                                    <  6:353679/149‑1 (MQ=255)
                                                                          aaaTGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTgg                                                                    >  5:164965/1‑149 (MQ=255)
                                                                          aaaTGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTgg                                                                    <  6:164965/149‑1 (MQ=255)
                                                                           aaTGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATg                                                                                                                   <  1:221761/101‑1 (MQ=255)
                                                                           aaTGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATg                                                                                                                   >  2:221761/1‑101 (MQ=255)
                                                                                  aaCAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTggg                                                                   <  3:38231/142‑1 (MQ=255)
                                                                                  aaCAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTggg                                                                   >  4:38231/1‑142 (MQ=255)
                                                                                  aaCAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGt                                                              >  7:122586/1‑147 (MQ=255)
                                                                                     aaTTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTaag                                                           >  6:70646/1‑147 (MQ=255)
                                                                                       ttGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGc                                                       <  2:276986/149‑1 (MQ=255)
                                                                                          cccTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCtt                                                     >  1:217409/1‑148 (MQ=255)
                                                                                                  aGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTa                                             <  5:366248/148‑1 (MQ=255)
                                                                                                  aGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTa                                             <  8:122586/148‑1 (MQ=255)
                                                                                                   ggTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAgg                                           <  5:70646/149‑1 (MQ=255)
                                                                                                    gTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGc                                          >  5:83060/1‑149 (MQ=255)
                                                                                                    gTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGc                                          <  2:217409/149‑1 (MQ=255)
                                                                                                         tCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGctc                                     >  6:91399/1‑149 (MQ=255)
                                                                                                           ccACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCtt                                   >  6:65288/1‑149 (MQ=255)
                                                                                                             aCCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTc                                  >  7:49126/1‑148 (MQ=255)
                                                                                                             aCCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTc                                  >  6:269602/1‑148 (MQ=255)
                                                                                                             aCCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCg                                 <  5:65288/149‑1 (MQ=255)
                                                                                                             aCCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTACGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCg                                 >  6:228002/1‑149 (MQ=255)
                                                                                                               cTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAg                               <  2:239536/149‑1 (MQ=255)
                                                                                                                  gTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGAc                            <  5:269602/149‑1 (MQ=255)
                                                                                                                    aTCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACt                           >  6:391974/1‑148 (MQ=255)
                                                                                                                     tctctcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTc                          <  8:49126/148‑1 (MQ=255)
                                                                                                                       tctcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGa                                >  8:141987/1‑140 (MQ=255)
                                                                                                                       tctcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGa                                <  4:21964/140‑1 (MQ=255)
                                                                                                                       tctcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGa                                <  7:141987/140‑1 (MQ=255)
                                                                                                                       tctcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGa                                >  3:21964/1‑140 (MQ=255)
                                                                                                                       tctcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGa                             <  1:180231/143‑1 (MQ=255)
                                                                                                                       tctcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGa                             >  2:180231/1‑143 (MQ=255)
                                                                                                                         tcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTcc                                                                          >  7:42643/1‑96 (MQ=255)
                                                                                                                         tcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTcc                                                                          <  8:42643/96‑1 (MQ=255)
                                                                                                                         tcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGtt                     >  1:383885/1‑149 (MQ=255)
                                                                                                                          aacacAAAGCGCATGTAACATTCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTc                                  >  7:187985/5‑135 (MQ=255)
                                                                                                                          aacacAAAGCGCATGTAACATTCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTc                                  <  8:187985/131‑1 (MQ=255)
                                                                                                                               aaaGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCcgg               >  2:91264/1‑149 (MQ=255)
                                                                                                                               aaaGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCcgg               <  5:391974/149‑1 (MQ=255)
                                                                                                                               aaaGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCcg                <  2:383885/148‑1 (MQ=255)
                                                                                                                                aaGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCcggc              <  5:91399/149‑1 (MQ=255)
                                                                                                                                aaGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGAGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCcggc              >  3:343360/1‑149 (MQ=255)
                                                                                                                                     aaTGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCCGGCGGtct         >  7:14595/1‑149 (MQ=255)
                                                                                                                                      aTGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCCGGCGGtctc        <  4:27352/149‑1 (MQ=255)
                                                                                                                                       tGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCCGGCGGTCTcc       >  8:101898/1‑149 (MQ=255)
                                                                                                                                        gtaaCATTCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTa                                             >  2:344215/4‑110 (MQ=255)
                                                                                                                                        gtaaCATTCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTa                                             <  1:344215/107‑1 (MQ=255)
                                                                                                                                            cATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGTGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCCGGCGGTCTCCTGCAt  <  4:343360/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
ACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCCGGCGGTCTCCTGCAT  >  NC_002947/2224258‑2224546

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: ‑