Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,577,151 Δ1 bp 6.2% coding (1244/1776 nt) PP_2259 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,577,1510A. 6.2% 206.3 / 6.9 65coding (1244/1776 nt)PP_2259sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base . (2/2);  total (32/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 = 9.30e-01

GGCCGGCGATGCGGTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGCCGAACAGTTC  >  NC_002947/2577005‑2577296
                                                                                                                                                   |                                                                                                                                                 
ggCCGGCGATGCGGTAGTACTGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTAcg                                                                                                                                                 >  6:10876/1‑149 (MQ=255)
 gCCGGCGATGCGGTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTAcgc                                                                                                                                                >  1:357425/1‑149 (MQ=255)
    ggCGATGCGGTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTAcgcg                                                                                                                                               >  6:374273/1‑147 (MQ=255)
    ggCGATGCGGTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCgca                                                                                                                                             <  4:342936/149‑1 (MQ=255)
         tGCGGTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAgg                                                                                                                                        >  4:360405/1‑149 (MQ=255)
           cGGTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCg                                                                                                                                      <  3:111109/149‑1 (MQ=255)
            ggTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCgg                                                                                                                                     <  5:374273/149‑1 (MQ=255)
             gtagtaCAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCggg                                                                                                                                    >  8:49347/1‑149 (MQ=255)
                   cAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCAAGCGGCTGGATGCTGCGCGACTGCACTACGCGCAGCAGGCGGGCCtggg                                                                                                                              >  2:323012/1‑149 (MQ=255)
                       tCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCtgggtggg                                                                                                                          <  2:287106/149‑1 (MQ=255)
                       tCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCtgggtggg                                                                                                                          <  7:49347/149‑1 (MQ=255)
                          tGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAg                                                                                                                       >  5:361169/1‑149 (MQ=255)
                          tGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAg                                                                                                                       >  2:24576/1‑149 (MQ=255)
                                    gcCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCg                                                                                                             >  7:185139/1‑149 (MQ=255)
                                        ggCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCaa                                                                                                         <  5:60271/149‑1 (MQ=255)
                                           gcgcACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAAtct                                                                                                      >  1:285314/1‑149 (MQ=255)
                                           gcgcACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAAtct                                                                                                      >  3:98785/1‑149 (MQ=255)
                                                       ggCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTg                                                                                                                     <  7:339898/122‑1 (MQ=255)
                                                       ggCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTg                                                                                                                     >  8:339898/1‑122 (MQ=255)
                                                        gCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATc                                                                                                  <  3:47087/140‑1 (MQ=255)
                                                        gCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATc                                                                                                  >  4:47087/1‑140 (MQ=255)
                                                        gCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAg                                                                                         <  2:357425/149‑1 (MQ=255)
                                                         ccAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAgg                                                                                        <  3:360405/149‑1 (MQ=255)
                                                         ccAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAgg                                                                                        <  2:105251/149‑1 (MQ=255)
                                                             gTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGc                                                                                    >  6:62708/1‑149 (MQ=255)
                                                                   aTTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCCGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGGGGGCAGTGGCATGTCGCCAATCTCATCCAGACTCAGGATGCCGTggg                                                                              >  8:100159/1‑148 (MQ=255)
                                                                   aTTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGt                                                                              >  5:19351/1‑149 (MQ=255)
                                                                     tgggcGCCGACACGACGCGGATATCCCCCGCCACCGGTTCGCCCCTGCCCAGCGGCTGGATGGTGGGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTc                                                                             <  1:323012/144‑1 (MQ=255)
                                                                      gCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCgg                                                                           >  2:160130/1‑149 (MQ=255)
                                                                          gCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGcc                                                                         <  1:160130/147‑1 (MQ=255)
                                                                            cGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGc                                                                                                                   >  4:44270/1‑103 (MQ=255)
                                                                            cGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGc                                                                                                                   <  3:44270/103‑1 (MQ=255)
                                                                            cGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGc                                                                     <  5:62708/149‑1 (MQ=255)
                                                                            cGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGc                                                                     >  4:254796/1‑149 (MQ=255)
                                                                              acacCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCgg                                                                   <  5:10876/149‑1 (MQ=255)
                                                                                accacGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGAt                                                                 <  5:59572/149‑1 (MQ=255)
                                                                                   acGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAg                                                              >  7:118947/1‑149 (MQ=255)
                                                                                    cgcgGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGt                                                             <  6:361169/149‑1 (MQ=255)
                                                                                     gcgGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCgg                                                                           >  1:182766/1‑134 (MQ=255)
                                                                                        gATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCg                                                         <  8:185139/149‑1 (MQ=255)
                                                                                             ccACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTg                                                    >  4:59644/1‑149 (MQ=255)
                                                                                               accggcaccggTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGcc                                                  >  5:229284/1‑149 (MQ=255)
                                                                                               accggcaccggTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGcc                                                  >  7:283145/1‑149 (MQ=255)
                                                                                                   gcaccggTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTtg                                              >  7:32196/1‑149 (MQ=255)
                                                                                                       cggTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTgtggt                                          <  3:254796/149‑1 (MQ=255)
                                                                                                        ggTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTgtggtg                                         <  6:229284/149‑1 (MQ=255)
                                                                                                         gTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTTGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGtggtgg                                        >  4:338845/1‑149 (MQ=255)
                                                                                                             ggctctgcgctgcgGCTGGATGCGGCGCGCCTGCTGTACGGGCAGCAGGCGGGGCTGGGTGGGCAGTGGCTTGTCGCCATTCTCATCCAGAAAGAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGgcgc                                    <  2:285314/145‑1 (MQ=255)
                                                                                                              ccACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCaa                                                                                                         >  3:337399/1‑79 (MQ=255)
                                                                                                              ccACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCaa                                                                                                         <  4:337399/79‑1 (MQ=255)
                                                                                                                     ccAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAAc                            <  1:69419/149‑1 (MQ=255)
                                                                                                                     ccAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAAc                            >  2:69419/1‑149 (MQ=255)
                                                                                                                         cGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACgcgc                        <  6:19351/149‑1 (MQ=255)
                                                                                                                               gATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGt                  >  8:224543/1‑149 (MQ=255)
                                                                                                                                    tgcGCGTACTGCAGTCGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAg             >  4:375356/1‑149 (MQ=255)
                                                                                                                                    ccgtagaaCTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATa              <  6:146413/142‑1 (MQ=255)
                                                                                                                                     gcgcgtaCTGCAGTCGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAAc                            >  2:212127/7‑133 (MQ=255)
                                                                                                                                     gcgcgtaCTGCAGTCGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGCGTTGCCCTTGTGGTGGGCGCCGGTGAAc                            <  1:212127/127‑1 (MQ=255)
                                                                                                                                      cgcgACTGCAGTACGCGCAGCAGGCGGGCCTGGGGGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGcc           <  4:123017/149‑1 (MQ=255)
                                                                                                                                         gACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGCCGaa        <  3:59644/149‑1 (MQ=255)
                                                                                                                                         gACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGCCGaa        <  8:283145/149‑1 (MQ=255)
                                                                                                                                           aCTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGc            >  1:225635/1‑144 (MQ=255)
                                                                                                                                           aCTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGc            <  2:225635/144‑1 (MQ=255)
                                                                                                                                            cTGCAGTCGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGCCGAACa      <  3:375356/148‑1 (MQ=255)
                                                                                                                                                aGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGCCGAACAGTTc  <  8:118947/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGCCGGCGATGCGGTAGTACAGGTCCTGGCGAAAGCGCCCGGCGCGCACTTCGTCGGCCAGGTTGCGATTGCTCGCCGACACCACGCGGATATCCACCGGCACCGGTTCGCCACTGCCCAGCGGCTGGATGCTGCGCGACTGCAGTACGCGCAGCAGGCGGGCCTGGGTGGGCAGTGGCATGTCGCCAATCTCATCCAGAAACAGGATGCCGTGGTCGGCCTTGCGGATCAGTCCCGGGTTGCCCTTGTGGTGGGCGCCGGTGAACGCGCCTTTGTCATAGCCGAACAGTTC  >  NC_002947/2577005‑2577296

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: