Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 627,782 G→T 5.9% P29P (CCC→CCA PP_0541 ← acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947627,7820GT 5.9% 214.6 / 2.9 67P29P (CCC→CCAPP_0541acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base T (2/2);  total (34/33)
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.58e-01

CCGCTTCAACGACAGCGACGAGTAGCTCGGGTACAGCTGGCAAAAGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTGGGAAACAGCGGAATGCCAGCAGCGTTCAAGGGCATGCT  >  NC_002947/627635‑627929
                                                                                                                                                     |                                                                                                                                                   
ccGCTTCAACGACAGCGACGAGTAGCTCGGGTACAGCTGGCAAAAGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAgg                                                                                                                                                    <  2:232068/149‑1 (MQ=255)
          gACAGCGACGAGTAGCTCGGGTACAGCTGGCAAAAGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGccc                                                                                                                                          <  2:14749/149‑1 (MQ=255)
          gACAGCGACGAGTAGCTCGGGTACAGCTGGCAAAAGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGccc                                                                                                                                          <  3:354361/149‑1 (MQ=255)
           aCAGCGACGAGTAGCTCGGGTACAGCTGGCAAAAGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGGGCCAGGTAACGCACCGATTCGTCCCGCTTCAGGCGCGTTTCCAACAACCTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCg                                                                                                                                         >  8:296350/1‑149 (MQ=255)
                    aGTAGCTCGGGTACAGCTGGCAAAAGCCCAGCAGATTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCa                                                                                                                                >  6:351926/1‑149 (MQ=255)
                               tACAGCTGGCAAAAGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGAc                                                                                                                     <  5:126615/149‑1 (MQ=255)
                                  aGCTGGCAAAAGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGaa                                                                                                                                    <  6:136176/131‑1 (MQ=255)
                                   gCTGGCAAAAGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGaa                                                                                                                  >  8:66121/1‑148 (MQ=255)
                                       gCAAAAGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACg                                                                                                             >  5:53795/1‑149 (MQ=255)
                                         aaaaGCCCATGAGTTTGCGCTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGACCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACggg                                                                                                           <  6:208341/149‑1 (MQ=255)
                                         aaaaGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACggg                                                                                                           <  7:66121/149‑1 (MQ=255)
                                         aaaaGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACggg                                                                                                           >  7:61226/1‑149 (MQ=255)
                                         aaaaGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCAATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACggg                                                                                                           >  8:29625/1‑149 (MQ=255)
                                          aaaGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTAGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACggg                                                                                                           <  1:92096/149‑1 (MQ=255)
                                                cagcagTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAg                                                                                                      >  6:356656/1‑147 (MQ=255)
                                                         gCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAg                                                                                             <  3:328591/148‑1 (MQ=255)
                                                                tcatcaTCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAg                                                                                    >  3:133032/1‑149 (MQ=255)
                                                                 catcatCCGGAAGTGCCAGGTCAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTCGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAgg                                                                                   <  7:296350/149‑1 (MQ=255)
                                                                 catcatCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGAAGGTCGAGGTGTTCCAgg                                                                                   >  2:222309/1‑149 (MQ=255)
                                                                  atcatcCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCATTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAggg                                                                                  >  3:279364/1‑149 (MQ=255)
                                                                    catcCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGc                                                                                                     >  8:64708/1‑128 (MQ=255)
                                                                     atcCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGt                                                                                 >  8:293107/1‑147 (MQ=255)
                                                                     atcCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGt                                                                                 <  7:293107/147‑1 (MQ=255)
                                                                     atcCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTg                                                                               >  4:130323/1‑149 (MQ=255)
                                                                     atcCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTg                                                                               >  3:234130/1‑149 (MQ=255)
                                                                     atcCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTg                                                                               >  4:79444/1‑149 (MQ=255)
                                                                          gAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAAGGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTtg                                                                          <  5:351926/149‑1 (MQ=255)
                                                                           aaGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTtga                                                                         <  6:53795/149‑1 (MQ=255)
                                                                                       aaTCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATgc                                                                   <  3:4026/143‑1 (MQ=255)
                                                                                       aaTCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATgc                                                                   >  4:4026/1‑143 (MQ=255)
                                                                                       aaTCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCAtgt                                                             <  8:61226/149‑1 (MQ=255)
                                                                                       aaTCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCAtgt                                                             <  4:133032/149‑1 (MQ=255)
                                                                                          cACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCg                                                          >  8:320881/1‑149 (MQ=255)
                                                                                           aCCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTATGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCgg                                                         >  3:151064/1‑149 (MQ=255)
                                                                                           aCCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTATGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCgg                                                         >  3:145495/1‑149 (MQ=255)
                                                                                             cGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGa                                                                                                                      <  4:35389/86‑1 (MQ=255)
                                                                                             cGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGa                                                                                                                      >  3:35389/1‑86 (MQ=255)
                                                                                                   gTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCAtgt                                                             <  6:318603/137‑1 (MQ=255)
                                                                                                   gTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCaga                                                 <  7:29625/149‑1 (MQ=255)
                                                                                                   gTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCa                                                   >  3:185540/1‑147 (MQ=255)
                                                                                                   gTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCa                                                   >  2:254381/1‑147 (MQ=255)
                                                                                                   gTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTCTTCCAGGGTTGCTTTGATGATGCGCAtgt                                                             >  5:318603/1‑137 (MQ=255)
                                                                                                    tCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCagag                                                <  1:254381/149‑1 (MQ=255)
                                                                                                      ccGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTtga                                                                         >  8:258265/1‑122 (MQ=255)
                                                                                                      ccGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTtga                                                                         <  7:258265/122‑1 (MQ=255)
                                                                                                      ccGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTc                                              >  6:96170/1‑149 (MQ=255)
                                                                                                      ccGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTAGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGt                                               >  7:303993/1‑148 (MQ=255)
                                                                                                       cGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTcc                                             <  3:130323/149‑1 (MQ=255)
                                                                                                        gCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTAACGGTATTTGACGAACAACGGGGGGAGCAGCTCGAGGTGTTCCAGGGTTGCTTTCATGATGCGCAAGTGCGGCCCTCAGACTCCa                                            >  3:324837/1‑149 (MQ=255)
                                                                                                         cTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTATGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCAt                                           <  4:145495/149‑1 (MQ=255)
                                                                                                         cTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTATGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCAt                                           <  4:151064/149‑1 (MQ=255)
                                                                                                            cAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCAttt                                         <  3:79444/148‑1 (MQ=255)
                                                                                                                        cAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTGGGAAACAGCgg                             <  5:356656/148‑1 (MQ=255)
                                                                                                                          aaaaaaCTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATgcgc                                                                 >  5:111200/1‑111 (MQ=255)
                                                                                                                          aaaaaaCTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATgcgc                                                                 <  6:111200/111‑1 (MQ=255)
                                                                                                                          aaaaaaCTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGc                                                           >  8:3254/1‑117 (MQ=255)
                                                                                                                          aaaaaaCTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGc                                                           <  7:3254/117‑1 (MQ=255)
                                                                                                                          aaaaaaCTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCaga                                                 >  4:270145/1‑127 (MQ=255)
                                                                                                                          aaaaaaCTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTgg                                       >  5:266332/1‑137 (MQ=255)
                                                                                                                          aaaaaaCTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTgg                                       <  6:266332/137‑1 (MQ=255)
                                                                                                                          aaaaaaCTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTGGGAAACAGCGGaa                           >  5:115339/1‑149 (MQ=255)
                                                                                                                                   cgcgAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTGGGAAACAGCGGAATGCcagca                   <  5:96170/148‑1 (MQ=255)
                                                                                                                                    gcgAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTGGGAAACAGCGGAATGCCagcagc                 >  6:170662/1‑149 (MQ=255)
                                                                                                                                     cgAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTGGGAAACAGCGGAATGCCAGCAGCg                >  7:262362/1‑149 (MQ=255)
                                                                                                                                       aGCTGTCGTGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTGGGAAACAGCGGAATGCCAGCAGCGtt              <  4:279364/149‑1 (MQ=255)
                                                                                                                                       aGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTGGGAAACAGCGGAATGCCAGCAGCGtt              <  5:340259/149‑1 (MQ=255)
                                                                                                                                             tCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTGGGAAACAGCGGAATGCCAGCAGCGTTCaa           <  7:320881/147‑1 (MQ=255)
                                                                                                                                             tCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTGGGAAACAGCGGAATGCCAGCAGCGTTCAAg          <  4:138030/148‑1 (MQ=255)
                                                                                                                                                     gggCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTGGGAAACAGCGGAATGCCAGCAGCGTTCAAGGGCAtgct  <  4:185540/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGCTTCAACGACAGCGACGAGTAGCTCGGGTACAGCTGGCAAAAGCCCAGCAGTTTGCCGTCGTCATCATCCGGAAGTGCCAGGTAAATCACCGATTCGTCCCGCTTCAGGCGCTTTTCCAAAAAACTGCGCGAGCTGTCGGGGTAGGGCAGTTGCCCGTAGAACTCACGGTATTTGACGAACAACGGGGTGAGCAGGTCGAGGTGTTCCAGGGTTGCTTTGATGATGCGCATGTGCGGCCCTCAGAGTCCATTTGGGAAACAGCGGAATGCCAGCAGCGTTCAAGGGCATGCT  >  NC_002947/627635‑627929

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: