Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,552,399:1 +G 5.1% coding (56/978 nt) PP_3137 → group 2 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,552,3991.G 5.1% 325.6 / 10.7 79coding (56/978 nt)PP_3137group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base . (35/40);  new base G (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCACGCGC  >  NC_002947/3552252‑3552546
                                                                                                                                                    |                                                                                                                                                    
gTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAg                                                                                                                                                     >  3:115929/1‑149 (MQ=255)
   gATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAgct                                                                                                                                                   <  5:353783/148‑2 (MQ=255)
   gATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAgcg                                                                                                                                                   <  7:149439/148‑1 (MQ=255)
       tgCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCTCGAGGCGCTgg                                                                                                                                              <  5:129280/148‑2 (MQ=255)
       ggCCTGACAATTCACTGGCTGGGGGGGTGCCGGCCATCATCAAGGCCCGGGATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGTTATTCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGcc                                                                                                                                             <  8:165716/148‑1 (MQ=255)
         ccTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCtt                                                                                                                                           >  5:303629/1‑149 (MQ=255)
          cTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTg                                                                                                                                          <  1:271785/149‑1 (MQ=255)
                aTTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGc                                                                                                                                    <  3:310972/149‑1 (MQ=255)
                  tCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTg                                                                                                                                  >  2:294735/1‑149 (MQ=255)
                    aCTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTgg                                                                                                                                 <  2:116373/148‑1 (MQ=255)
                    aCTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGgc                                                                                                                                >  2:111992/1‑149 (MQ=255)
                     cTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTCAGGCCTATGGCAACGGCGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGgcg                                                                                                                               >  3:316723/1‑149 (MQ=255)
                        gCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCACCGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCAc                                                                                                                            >  5:281218/1‑149 (MQ=255)
                        gCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCAc                                                                                                                            >  5:207645/1‑149 (MQ=255)
                              ggggTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGcc                                                                                                                      <  8:43591/149‑1 (MQ=255)
                               gggTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCg                                                                                                                     <  1:353201/149‑1 (MQ=255)
                                    ccGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACaaa                                                                                                                <  2:120342/149‑1 (MQ=255)
                                    ccGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACa                                                                                                                  >  4:270752/1‑147 (MQ=255)
                                      ggCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGt                                                                                                              >  5:143405/1‑149 (MQ=255)
                                        ccATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCa                                                                                                            >  1:374790/1‑149 (MQ=255)
                                         catcatCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCATTGTTTCGCTGGCGCTCGGGGCCGACAAAGTCAt                                                                                                           >  6:87901/1‑149 (MQ=255)
                                             atcaAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTAc                                                                                                       >  4:232002/1‑149 (MQ=255)
                                              tcaAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGCGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACg                                                                                                      >  3:331448/1‑149 (MQ=255)
                                              tcaAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACg                                                                                                      >  5:267247/1‑149 (MQ=255)
                                                  ggCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGa                                                                                                  <  6:303629/149‑1 (MQ=255)
                                                     ccGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCttca                                                                                                                                         >  4:364111/1‑105 (MQ=255)
                                                     ccGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCttca                                                                                                                                         <  3:364111/107‑3 (MQ=255)
                                                     ccGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTc                                                                                               >  7:76936/1‑149 (MQ=255)
                                                     ccGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTc                                                                                               <  4:239841/149‑1 (MQ=255)
                                                      cGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGGTGGCGCACGGGGCCGACAAAGTCATGTACGTCGCCTCg                                                                                              >  6:199613/1‑149 (MQ=255)
                                                         cATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGCGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCgtgt                                                                                           <  4:316723/149‑1 (MQ=255)
                                                          aTCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGtt                                                                                          <  2:287899/149‑1 (MQ=255)
                                                              ggCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGTGTTCCAc                                                                                      <  6:207645/149‑1 (MQ=255)
                                                              ggCCGCCGGGGGATGTGTGGAGTGTGGCCTTTGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCTCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCAc                                                                                      <  8:334455/149‑1 (MQ=255)
                                                                         cttCGGTCGAGTGAGGCCTATGGCAACTGTTATGGTGGGGGTGGTAATCGGCATCAATGAAGGCCAGCGCCTCGAGCGGTGCGTTGTTTCTCTGGCTCACGGGGCCGACTAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTctc                                                                           <  4:331448/147‑3 (MQ=255)
                                                                          aTCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCt                                                                          <  4:115929/149‑1 (MQ=255)
                                                                          aTCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCt                                                                          >  7:27167/1‑149 (MQ=255)
                                                                             ggTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCACACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCa                                                                       <  1:294735/149‑1 (MQ=255)
                                                                              gTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAg                                                                      <  3:270752/149‑1 (MQ=255)
                                                                                     gAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGgccg                                                               >  5:1263/1‑149 (MQ=255)
                                                                                      aGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCg                                                                            <  5:245504/135‑1 (MQ=255)
                                                                                      aGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCg                                                                            >  6:245504/1‑135 (MQ=255)
                                                                                        gtcTATGGCAACGGTGATCCGGGTGGTGGTAATCGGCGGCAATGAAGGCCAGCCCCTCGAGCGCAGCCTTGTTTCGCTGGCGCCAGGGGCCGACAAAGTCATGTCCGTCGACTCGGGTTCCACCGAAGGTTCGGTGCAGTAGgccgcc                                                             >  6:110858/3‑148 (MQ=255)
                                                                                          cTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCt                                                                          >  4:56337/1‑133 (MQ=255)
                                                                                          cTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCt                                                                          <  3:56337/133‑1 (MQ=255)
                                                                                           tATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAAcc                                                         >  4:251943/1‑149 (MQ=255)
                                                                                                 aaCGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCg                                                   >  2:263783/1‑149 (MQ=255)
                                                                                                  aCGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGt                                                  <  3:361845/149‑1 (MQ=255)
                                                                                                   cGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGgccgcc                                                             <  1:164492/137‑1 (MQ=255)
                                                                                                   cGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGgccgcc                                                             >  2:164492/1‑137 (MQ=255)
                                                                                                      tGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCa                                              <  2:374790/149‑1 (MQ=255)
                                                                                                        aTCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCtt                                                       >  5:183020/1‑138 (MQ=255)
                                                                                                        aTCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCtt                                                       <  6:183020/138‑1 (MQ=255)
                                                                                                        aTCCGGGTGGTGGTAATCGGCCGCAATGAAGGCCATCGCCTCGAGCGCAGCCTTGTTTCGCTGGCGCACGCGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTAGGCCGCCAACCTTGGCGTGGCAg                                             >  2:315921/1‑148 (MQ=255)
                                                                                                           ggggTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACt                                         <  6:1263/149‑1 (MQ=255)
                                                                                                            gggTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTg                                        >  7:93070/1‑149 (MQ=255)
                                                                                                             ggtggtggtAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCCACCTTGGCGTGGCAGTACTgg                                       >  5:263540/1‑149 (MQ=255)
                                                                                                                     tAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGTCCGACAAAGTCAGGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGAc                               <  3:232002/149‑1 (MQ=255)
                                                                                                                      actCGGCCCCAATGAAGGCCAGCGCCGCGAGCGCTGCCTTGTTTCGCTGCCGCCCGGGGCCGACAAAGTCATGTCCGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGGGGCAGTACTGGCGCTGGACa                              <  6:263540/147‑1 (MQ=255)
                                                                                                                       ttCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACATAGTCATGTACGTCGACTGGGGTTCCACCGTTGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACAt                             <  5:110858/148‑1 (MQ=255)
                                                                                                                       aTCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACa                              <  1:263783/148‑1 (MQ=255)
                                                                                                                       aTCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACAt                             <  8:76936/149‑1 (MQ=255)
                                                                                                                       aTCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACAt                             <  1:111992/149‑1 (MQ=255)
                                                                                                                          ggCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                           <  3:251943/148‑1 (MQ=255)
                                                                                                                          ggCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                           <  8:27167/148‑1 (MQ=255)
                                                                                                                          ggCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCCCACGGGGCCGACAAACTCATATACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                           <  6:281218/148‑1 (MQ=255)
                                                                                                                             cGCAATGAAGGCCAGCGCGTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCg                                                   <  8:122112/121‑1 (MQ=255)
                                                                                                                             cGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCg                                                   >  7:122112/1‑121 (MQ=255)
                                                                                                                             cGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACa                              >  6:120076/1‑142 (MQ=255)
                                                                                                                             cGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACa                              <  5:120076/142‑1 (MQ=255)
                                                                                                                                   gAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGAc                               >  1:376810/1‑135 (MQ=255)
                                                                                                                                   gAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGAc                               <  2:376810/135‑1 (MQ=255)
                                                                                                                                   gAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGtt                 >  7:281176/1‑149 (MQ=255)
                                                                                                                                      ggCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAc              <  6:267247/149‑1 (MQ=255)
                                                                                                                                      aggcAGCGCTCGAGGCGCTGGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCa               >  5:218248/4‑149 (MQ=255)
                                                                                                                                      aggcAGCGCTCGAGGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                           <  5:81186/133‑1 (MQ=255)
                                                                                                                                      aggcAGCGCTCGAGGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                           >  6:81186/4‑136 (MQ=255)
                                                                                                                                        ccAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAccg            >  3:71188/1‑149 (MQ=255)
                                                                                                                                          aGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTAGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAccgcc          >  2:245382/1‑149 (MQ=255)
                                                                                                                                                 cGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGAc                                                                                                 >  3:130497/1‑55 (MQ=255)
                                                                                                                                                 cGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGAc                                                                                                 <  4:130497/55‑1 (MQ=255)
                                                                                                                                                  gAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCACgcgc  <  8:281176/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCACGCGC  >  NC_002947/3552252‑3552546

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