Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,970,603 T→C 5.5% intergenic (‑361/+33) 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,6030TC 5.5% 214.8 / 2.4 73intergenic (‑361/+33)fliC/PP_4379flagellin/3‑oxoacyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base T (31/38);  new base C (2/2);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

ATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAA  >  NC_002947/4970457‑4970745
                                                                                                                                                    |                                                                                                                                              
aTCCCGGCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTg                                                                                                                                               <  1:355668/148‑1 (MQ=255)
 tCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTGATGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTgg                                                                                                                                              >  1:116150/1‑149 (MQ=255)
     gTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGgccgcc                                                                                                                                         <  3:252163/149‑1 (MQ=255)
       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                                                                               <  6:144710/149‑1 (MQ=255)
                                                                cGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGa                                                                               >  1:363828/1‑149 (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                                                                                                            <  8:111847/104‑1 (MQ=255)
                                                                                aCCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAg                                                                                                            >  7:111847/1‑104 (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                             <  6:149636/149‑1 (MQ=255)
                                                                                                                  cTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGa                             <  1:173687/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                          <  2:251266/142‑1 (MQ=255)
                                                                                                                             gCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCa                          >  1:251266/1‑142 (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)
                                                                                                                                                    |                                                                                                                                              
ATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCCCACGACAGGCCCACGCCAAAGCCACTGATCGCCACGCGCTTCCAGCTGCCATCGAGCATGTGCTTCTGCAGCAGCAGCGGAA  >  NC_002947/4970457‑4970745

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: