Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,738,598 C→T 5.5% A199T (GCT→ACT)  PP_1532 ← putative Phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,738,5980CT 5.5% 224.8 / 2.1 72A199T (GCT→ACT) PP_1532putative Phage integrase
Reads supporting (aligned to +/- strand):  ref base C (29/39);  new base T (2/2);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGCCCCCATGAAAAGGTCGCGCAGTAAGGAG  >  NC_002947/1738450‑1738740
                                                                                                                                                     |                                                                                                                                              
cTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGc                                                                                                                                                <  6:36346/149‑1 (MQ=255)
 tGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGc                                                                                                                                                >  6:18017/1‑148 (MQ=255)
      cTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGcctcct                                                                                                                                           >  7:84288/1‑148 (MQ=255)
      cTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTAATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGcctcct                                                                                                                                           >  4:35069/1‑148 (MQ=255)
       ttAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCa                                                                                                                                         >  6:186133/1‑149 (MQ=255)
          gCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGTTTGAGCCATGGCTGGGCAGCCTCCTCATAg                                                                                                                                      <  1:245694/149‑1 (MQ=255)
           ccTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCCGGTGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAg                                                                                                                                      <  1:126365/149‑1 (MQ=255)
                cccccGTCGTCGATTGGACGACCTGCAGGAATCCGTCGACGTCATGCCGGAAAACGAAGGGAGCAATGGCCT‑CGCGGCGCTGGCCGGTTAGAAGAGCGAGCTCGATGGCGCGCTTGAGGCATCGCTGGGCAGCGTCGACATAGATGGcg                                                                                                                                >  8:343178/1‑148 (MQ=255)
                     gTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCca                                                                                                                           <  6:148947/149‑1 (MQ=255)
                     gTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCca                                                                                                                           <  8:16244/149‑1 (MQ=255)
                        ttCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCa                                                                                                                                         <  4:67671/132‑1 (MQ=255)
                        ttCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCcaca                                                                                                                         >  4:107638/1‑148 (MQ=255)
                          cGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACaa                                                                                                                        <  6:31443/147‑1 (MQ=255)
                           gATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCAt                                                                                                                                        >  3:354391/1‑130 (MQ=255)
                           gATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCAt                                                                                                                                        <  4:354391/130‑1 (MQ=255)
                           gATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAActc                                                                                                                     >  5:142255/1‑149 (MQ=255)
                             ttGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTa                                                                                                                   <  3:107638/149‑1 (MQ=255)
                                gacgacCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGc                                                                                                                 <  5:186133/148‑1 (MQ=255)
                                  cgacCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTc                                                                                                              <  6:365104/149‑1 (MQ=255)
                                    acCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAg                                                                                                            <  8:84288/149‑1 (MQ=255)
                                     ccTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAgg                                                                                                           <  3:347160/149‑1 (MQ=255)
                                       tGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAggcg                                                                                                         >  5:404338/1‑149 (MQ=255)
                                           ggAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAggcgctgg                                                                                                     <  6:404338/149‑1 (MQ=255)
                                            gAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAggcgctggc                                                                                                    >  8:221658/1‑149 (MQ=255)
                                               tCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAggcgctggcgct                                                                                                 <  7:343178/149‑1 (MQ=255)
                                                 cGTCCACCTCATCCCGGAAACGGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGTTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTg                                                                                               <  4:357638/149‑1 (MQ=255)
                                                 cGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGGCCGGTTATAAGAGCGAGCTACATGGCGCTCTTGAGCCATGGCTGGGCAGCCGCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTg                                                                                               <  6:134146/149‑1 (MQ=255)
                                                     cACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCt                                                                                           <  3:35069/149‑1 (MQ=255)
                                                              ccGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTg                                                                                  >  8:314450/1‑149 (MQ=255)
                                                                    taCGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGATTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCtt                                                                            <  6:105085/148‑1 (MQ=255)
                                                                     aCGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCttt                                                                           <  7:221658/149‑1 (MQ=255)
                                                                      cGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGATCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTg                                                                          >  7:261946/1‑149 (MQ=255)
                                                                       gAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTgg                                                                         <  6:142255/149‑1 (MQ=255)
                                                                         aTGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGt                                                                        <  8:130453/148‑1 (MQ=255)
                                                                              gCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTc                                                                  <  7:311674/149‑1 (MQ=255)
                                                                              gCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTc                                                                  >  8:212911/1‑149 (MQ=255)
                                                                                 aTGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGGTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGAc                                                               >  2:132202/1‑149 (MQ=255)
                                                                                 aTGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCtt                                                                   <  6:141182/145‑1 (MQ=255)
                                                                                 aTGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCtt                                                                   >  5:141182/1‑145 (MQ=255)
                                                                                  tGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGTCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACt                                                              >  7:117741/1‑149 (MQ=255)
                                                                                  tGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGTCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACt                                                              >  7:123535/1‑149 (MQ=255)
                                                                                      cT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTa                                                                                                                   >  6:208957/1‑92 (MQ=255)
                                                                                      cT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTa                                                                                                                   <  5:208957/92‑1 (MQ=255)
                                                                                      cT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGt                                                          >  8:365260/1‑149 (MQ=255)
                                                                                             gcgcTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGt                                                    >  3:394921/1‑149 (MQ=255)
                                                                                                   cccGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGAcc                                               <  5:18017/148‑1 (MQ=255)
                                                                                                    ccGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGc                                                                                                                 <  3:318348/81‑1 (MQ=255)
                                                                                                    ccGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGc                                                                                                                 >  4:318348/1‑81 (MQ=255)
                                                                                                     cGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAg                                            <  1:270026/149‑1 (MQ=255)
                                                                                                           taAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGAcc                                               <  7:215564/140‑1 (MQ=255)
                                                                                                           taAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGAcc                                               >  8:215564/1‑140 (MQ=255)
                                                                                                           taAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTg                                      <  3:92660/149‑1 (MQ=255)
                                                                                                            aaGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGt                                                                        >  7:41572/1‑114 (MQ=255)
                                                                                                            aaGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGt                                                                        <  8:41572/114‑1 (MQ=255)
                                                                                                             agagCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTgg                                     <  7:212911/148‑1 (MQ=255)
                                                                                                              gagcgagcTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGc                                    >  4:179161/1‑148 (MQ=255)
                                                                                                                 cgagcTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATg                                >  7:36742/1‑149 (MQ=255)
                                                                                                                  gagcTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATgg                               >  4:63501/1‑149 (MQ=255)
                                                                                                                  gagcTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATgg                               <  4:394921/149‑1 (MQ=255)
                                                                                                                    gcTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGcc                             >  7:280971/1‑149 (MQ=255)
                                                                                                                    gcTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGcc                             <  7:365260/149‑1 (MQ=255)
                                                                                                                        cATGGCGCGCTTGAGCCATGGCTGGGCAGTCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGCCCCCa                         <  8:117741/149‑1 (MQ=255)
                                                                                                                        cATGGCGCGCTTGAGCCATGGCTGGGCAGTCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGCCCCCa                         <  8:123535/149‑1 (MQ=255)
                                                                                                                        cATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGccccc                          <  8:280971/148‑1 (MQ=255)
                                                                                                                         aTGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTggg                                                           <  8:206387/114‑1 (MQ=255)
                                                                                                                         aTGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTggg                                                           >  7:206387/1‑114 (MQ=255)
                                                                                                                           ggCGCGCTTCAGCCAAGGCTGCCCAGCCTCCTCATAGATCGCATTCCACAACTCTAGCCTCAGTCGCTGGCGCTTAATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTAGCTTCTATTCGACCCAGCCCCTGCCCATGCCCCCCATg                       >  4:366328/1‑148 (MQ=255)
                                                                                                                                  ttGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGCCCCCATGAAAAGGTc               >  5:152715/1‑149 (MQ=255)
                                                                                                                                    gAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGCCCCCATGAAAAGGTcgc             >  6:204366/1‑149 (MQ=255)
                                                                                                                                      gCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGCCCCCATGAAAAGGTcgcg            >  3:121411/1‑148 (MQ=255)
                                                                                                                                       ccATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGCCCCCATGAAAAGGTcgc             >  4:225722/1‑146 (MQ=255)
                                                                                                                                       ccATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGCCCCCATGAAAAGGTcgc             <  3:225722/146‑1 (MQ=255)
                                                                                                                                       ccATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGCCCCCATGAAAAGGTCGCGCa          <  5:204366/149‑1 (MQ=255)
                                                                                                                                          tGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGCCCCCATGAAAAGGTCGCGCAGTa       <  8:261946/149‑1 (MQ=255)
                                                                                                                                               gggCAGCATCCTCATAGATAGCTATCAACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCAGTGCGGCTTTTGTAGCTTCTACTGGGTTGTAGTCGACCCAGCCCCTGGCCATAGCCACCATGAAAAGGTCGCGCAGAAAGGAg  >  8:271811/1‑149 (MQ=255)
                                                                                                                                                ggCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGCCCCCATGAAAAGGTCGCGCAGTAAGGa   <  6:152715/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
CTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCT‑CGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTTGAGCCATGGCTGGGCAGCCTCCTCATAGATCGCTTTCCACAACTCTAGCGTCAGGCGCTGGCGCTTGATCTTAACCCGTGCGGCTTTGGTAGCTTCGACTGGGTTGTAGTCGACCCAGCCCCTGGCCATGGCCCCCATGAAAAGGTCGCGCAGTAAGGAG  >  NC_002947/1738450‑1738740

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