Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,970,611 G→C 6.0% intergenic (‑369/+25) fliC ← / ← PP_4379 flagellin/3‑oxoacyl‑ACP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,970,6110GC 6.0% 214.4 / 3.2 67intergenic (‑369/+25)fliC/PP_4379flagellin/3‑oxoacyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base G (29/34);  new base C (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

CGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTGGAG  >  NC_002947/4970464‑4970753
                                                                                                                                                    |                                                                                                                                              
cGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCgg                                                                                                                                               <  3:308777/149‑1 (MQ=255)
 ggTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCggg                                                                                                                                              <  2:296983/149‑1 (MQ=255)
    ggAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCGTACATTGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGgcgc                                                                                                                                           <  6:322479/149‑1 (MQ=255)
     gacgggatGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCt                                                                                                                                          <  1:305199/142‑1 (MQ=255)
      agggggTGAGGCTGCTAAGAAGGATACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCtt                                                                                                                                         <  4:264940/146‑1 (MQ=255)
         gggTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCCGGCGGCCACCGCCCGGCGCttttt                                                                                                                                      <  1:363824/149‑2 (MQ=38)
           gTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAgg                                                                                                                                    >  1:290821/1‑149 (MQ=255)
           gTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCCGGCGGCCACCGCCCGGCGCttttttt                                                                                                                                    <  6:80209/149‑4 (MQ=38)
               ggCTGTGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCt                                                                                                                                >  4:392318/1‑149 (MQ=255)
                  tGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGAc                                                                                                                             >  6:66743/1‑149 (MQ=255)
                           ggAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAg                                                                                                                    >  5:74987/1‑149 (MQ=255)
                           ggAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAg                                                                                                                    >  4:354501/1‑149 (MQ=255)
                              gccgcgATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTTAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAggg                                                                                                                                   <  3:21127/129‑1 (MQ=255)
                                cgcgATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgcg                                                                                                               <  5:285319/149‑1 (MQ=255)
                                cgcgATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgcg                                                                                                               <  6:74987/149‑1 (MQ=255)
                                cgcgAGGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgcg                                                                                                               <  3:226450/149‑1 (MQ=255)
                                 gcgATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCgg                                                                                                              >  4:348196/1‑149 (MQ=255)
                                   gATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTc                                                                                                                  <  5:398793/143‑1 (MQ=255)
                                   gATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTc                                                                                                                  >  6:398793/1‑143 (MQ=255)
                                       ggTGGCACCGGCTTTGCCGGTGGTCGCGGGTGAACCCGCGCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAgg                                                                                                        <  2:196891/149‑1 (MQ=255)
                                       gcgaacACCGGCAAAGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAg                                                                                                         >  1:206506/6‑148 (MQ=255)
                                        gTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgcg                                                                                                               <  3:181776/141‑1 (MQ=255)
                                        gTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcgcg                                                                                                               >  4:181776/1‑141 (MQ=255)
                                         tGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGAt                                                                                                      <  5:9677/149‑1 (MQ=255)
                                          ggCACCGGCTTTGCCGGTGTTCACGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATg                                                                                                     <  4:271957/149‑1 (MQ=255)
                                               cGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGcc                                                                                                  >  5:144710/1‑147 (MQ=255)
                                                 gCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGc                                                                                              <  8:230863/149‑1 (MQ=255)
                                                       ccGGTGTACGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCAc                                                                                         <  7:17978/148‑1 (MQ=255)
                                                        cGGTGTTCGCGGGTGCACCCGCTCCTACAATGGGCCGCGTTAGGCAGAAGCTTATTTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCggg                                                                                                                                              <  1:71664/95‑1 (MQ=255)
                                                        cGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCggg                                                                                                                                              >  2:71664/1‑94 (MQ=255)
                                                        cGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGa                                                                                       >  1:363828/1‑149 (MQ=255)
                                                        cGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGa                                                                                       <  6:144710/149‑1 (MQ=255)
                                                         ggTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGAc                                                                                      <  7:264865/149‑1 (MQ=255)
                                                            attCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAgg                                                                                   >  5:327752/2‑149 (MQ=255)
                                                               cgcgGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGATGACTGGGATCAGTCGCGGTCCAGGATGGCCGAGCCCCACGCCAGGccc                                                                                >  4:244087/1‑149 (MQ=255)
                                                                gcgGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCa                                                                               >  5:96647/1‑149 (MQ=255)
                                                                    gTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGcc                                                                           <  6:96647/149‑1 (MQ=255)
                                                                        aCCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAg                                                                                                                    >  7:111847/1‑104 (MQ=255)
                                                                        aCCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAg                                                                                                                    <  8:111847/104‑1 (MQ=255)
                                                                              tCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACt                                                                  >  2:173687/1‑148 (MQ=255)
                                                                                cTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGTGCGGTGGCCGACGGGCGCTTTAAGGGAATTAATGGGATCATTCGCGTTACAGGATGGCCGAGCCCAACCACAGGCCAACGCCAAAGCCaataat                                                               >  4:13216/1‑144 (MQ=255)
                                                                                cTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGAt                                                               >  5:149636/1‑149 (MQ=255)
                                                                                  aCAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGAAGAAAAAAAAGAGCAGGGCGGTGGCCGCCGGGCGCTTATAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCg                                                             >  8:247603/1‑149 (MQ=255)
                                                                                   cAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGc                                                            <  3:152002/149‑1 (MQ=255)
                                                                                   cAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGc                                                            <  2:363828/149‑1 (MQ=255)
                                                                                    aaTGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCCATGATCGcc                                                           >  3:252515/1‑149 (MQ=255)
                                                                                      tGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCAc                                                         >  5:111903/1‑149 (MQ=255)
                                                                                        ggCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCAcgc                                                       <  2:290821/149‑1 (MQ=255)
                                                                                         gCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCAcgcg                                                      >  2:408042/1‑149 (MQ=255)
                                                                                            gCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCCCGCGCtt                                                   >  7:277614/1‑149 (MQ=255)
                                                                                             cGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCt                                                    <  6:327752/147‑1 (MQ=255)
                                                                                                  ggCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGc                                              <  3:392318/148‑1 (MQ=255)
                                                                                                   gCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTg                                            <  8:54989/149‑1 (MQ=255)
                                                                                                    cAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGACGAGCCCCACGACAGGCCCACGCCAAAGCCAGTGATCGCCACGCGCTTCCAGCTGc                                           >  8:188202/1‑149 (MQ=255)
                                                                                                     aGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGc                                           >  3:127622/1‑148 (MQ=255)
                                                                                                       aaGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCAt                                        >  6:52260/1‑149 (MQ=255)
                                                                                                        aGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATc                                       <  8:154286/149‑1 (MQ=255)
                                                                                                          cTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGa                                     <  1:173687/149‑1 (MQ=255)
                                                                                                          cTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGa                                     <  6:149636/149‑1 (MQ=255)
                                                                                                          cTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGa                                     <  6:111903/149‑1 (MQ=255)
                                                                                                               ttttCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCAt                                 >  3:41715/1‑149 (MQ=255)
                                                                                                                     gcattaaaaaaaGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCtt                           <  3:13216/144‑1 (MQ=255)
                                                                                                                     gCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCa                                  >  1:251266/1‑142 (MQ=255)
                                                                                                                     gCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCa                                  <  2:251266/142‑1 (MQ=255)
                                                                                                                       aGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCATCTGCCATCGAGCATGTGCTTCt                         >  8:329271/1‑149 (MQ=255)
                                                                                                                       aGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCt                         <  2:206506/149‑1 (MQ=255)
                                                                                                                           ctaaaaGCGCCCGGCGGCCACCGCCCGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTgcag                     >  1:31586/3‑149 (MQ=39)
                                                                                                                               aaGCGCCCGGCGGCCACCGCCCGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTgcagcag                  >  5:87723/1‑148 (MQ=255)
                                                                                                                                  cgcCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCagcagcagc              <  7:370013/149‑1 (MQ=255)
                                                                                                                                      gggCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGAGCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGaa          <  2:116150/149‑1 (MQ=255)
                                                                                                                                           gccaccgcccGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTg     <  6:87723/145‑1 (MQ=255)
                                                                                                                                              gccgccGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTggag  <  4:127622/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAACACTGGAG  >  NC_002947/4970464‑4970753

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