Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,970,517 T→A 6.7% intergenic (‑275/+119) 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,5170TA 6.7% 207.4 / 2.4 69intergenic (‑275/+119)fliC/PP_4379flagellin/3‑oxoacyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base A (3/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

GCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCC  >  NC_002947/4970375‑4970663
                                                                                                                                               |                                                                                                                                                  
gCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGgt                                                                                                                                              <  5:94221/149‑1 (MQ=255)
gCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGgt                                                                                                                                              >  2:179262/1‑149 (MQ=255)
  aTGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCAAAGCCGgtgc                                                                                                                                            <  1:369391/149‑2 (MQ=255)
    ggggCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTc                                                                                                                                          >  5:362905/1‑149 (MQ=255)
    ggggCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCAAAGCCGgtg                                                                                                                                             <  1:217035/146‑1 (MQ=255)
    ggggCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCAAAGCCGgtg                                                                                                                                             >  2:217035/1‑146 (MQ=255)
      ggCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTcgc                                                                                                                                        >  7:247171/1‑149 (MQ=255)
       gCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTcgcg                                                                                                                                       >  8:132096/1‑149 (MQ=255)
        ccAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                                      <  6:362905/149‑1 (MQ=255)
          aGCCCTTATAATTCTTGGTGTTAAGTGCTGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGTTGGAGCGGGTGTGGCTGTGAAGATGTTGACGCGTTGGGTGGTTCCGGCTTTGCCGGTGTTCGCGGGt                                                                                                                                    <  7:243331/149‑1 (MQ=255)
              cTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGTATGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGaa                                                                                                                                 >  1:106366/1‑149 (MQ=255)
              cTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAAc                                                                                                                                >  1:258181/1‑149 (MQ=255)
               tGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAAcc                                                                                                                               >  4:231710/1‑149 (MQ=255)
                gagaAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAAccc                                                                                                                              >  6:28651/1‑149 (MQ=255)
                  gaAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGc                                                                                                                            <  8:359593/149‑1 (MQ=255)
                    aaTCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTc                                                                                                                          <  2:258181/149‑1 (MQ=255)
                    aaTCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTc                                                                                                                          <  5:28651/149‑1 (MQ=255)
                     aTCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTcc                                                                                                                         <  8:313721/149‑1 (MQ=255)
                     aTCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTcc                                                                                                                         >  7:313721/1‑149 (MQ=255)
                        aaGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTa                                                                                                                       <  6:351207/148‑1 (MQ=255)
                             gTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATgg                                                                                                                 >  5:374045/1‑149 (MQ=255)
                             gTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATgg                                                                                                                 >  7:188319/1‑149 (MQ=255)
                              ttAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATggg                                                                                                                <  3:269007/149‑1 (MQ=255)
                                  gggCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGc                                                                                                            <  5:215504/149‑1 (MQ=255)
                                  gggCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGc                                                                                                            >  5:177387/1‑149 (MQ=255)
                                        tttATTGGCGTCATAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAgg                                                                                                      <  7:132096/149‑1 (MQ=255)
                                           aTTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAgg                                                                                                      >  6:90496/1‑146 (MQ=255)
                                                 gTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCtt                                                                                             <  2:161839/149‑1 (MQ=255)
                                                 gTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCt                                                                                              >  4:147916/1‑148 (MQ=255)
                                                      aaaCTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTAtttt                                                                                        <  7:29569/149‑1 (MQ=255)
                                                       aaCTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTc                                                                                       <  3:141876/149‑1 (MQ=255)
                                                       aaCTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTc                                                                                       <  2:103962/149‑1 (MQ=255)
                                                        aCTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCg                                                                                      <  6:177387/149‑1 (MQ=255)
                                                                aaCTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGaaa                                                                               >  2:21550/1‑148 (MQ=255)
                                                                 aCTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGCGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGaaaa                                                                              >  2:122362/1‑148 (MQ=255)
                                                                     ggCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGaaaaaaaa                                                                          <  3:231710/148‑1 (MQ=255)
                                                                     ggCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAg                                                                         >  2:366627/1‑149 (MQ=255)
                                                                              gCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGATGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCCCTCCTACAATGGGCTCCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAACCGCCGGGCg                                                                >  3:224292/1‑149 (MQ=38)
                                                                                cTCATCCCGTCGGTGGAGCGGGTGAGGCTGGGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCgg                                                               >  1:34398/1‑149 (MQ=255)
                                                                                  cATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTg                                                             >  8:318288/1‑149 (MQ=255)
                                                                                   aTCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGt                                                              <  8:131044/147‑1 (MQ=255)
                                                                                     cccGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGgcc                                                          >  5:88442/1‑149 (MQ=255)
                                                                                      ccGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGgcc                                                          <  1:21550/148‑1 (MQ=255)
                                                                                        gTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGgccgcc                                                       >  6:40574/1‑149 (MQ=255)
                                                                                             tGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGgcg                                                  >  7:170280/1‑149 (MQ=255)
                                                                                               gAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCt                                                >  6:221892/1‑149 (MQ=255)
                                                                                               gAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCt                                                <  2:106366/149‑1 (MQ=255)
                                                                                               gAGCGGGTGAGGCTGCGAAGAAGGAGACGCGACGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCt                                                >  4:200509/1‑149 (MQ=255)
                                                                                                aGCGGGGGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCtt                                               <  8:247171/149‑1 (MQ=255)
                                                                                                  cGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCttt                                              <  2:34398/148‑1 (MQ=255)
                                                                                                      tGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAggg                                         <  1:179262/149‑1 (MQ=255)
                                                                                                               gaagaaGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTgg                                <  6:88442/149‑1 (MQ=255)
                                                                                                               gaagaaGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTgg                                >  7:335471/1‑149 (MQ=255)
                                                                                                               caaccagGAGACGCGAGGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGGTTGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCCGGCGGCCACCGCCCGGCGCttttt                                            <  5:128705/132‑2 (MQ=25)
                                                                                                                aagaagGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTggg                               <  2:258530/149‑1 (MQ=255)
                                                                                                                  gaagGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGAt                             <  7:318288/149‑1 (MQ=255)
                                                                                                                      gagaCGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCGTATTTTCGGCAGAAAAAAAAGAGCCGGGCCGTGGCCGCCGGGCGCATTTAGGGGCTGACTGGGATCAGt                         >  8:166608/1‑149 (MQ=255)
                                                                                                                       agaCGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGTCTGACTGGGATCAGTc                        >  4:76281/1‑149 (MQ=255)
                                                                                                                        gaCGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTcg                       >  2:276585/1‑149 (MQ=255)
                                                                                                                               tGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACa                <  6:374045/149‑1 (MQ=255)
                                                                                                                                cgcgaacACCGGCAAAGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAg               >  1:197442/7‑149 (MQ=39)
                                                                                                                                 ggTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAg               <  5:40574/148‑1 (MQ=255)
                                                                                                                                 gcgaacACCGGCAAAGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAg               >  8:215621/6‑148 (MQ=255)
                                                                                                                                   tGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGAt            >  4:104536/1‑149 (MQ=255)
                                                                                                                                   tGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGAt            >  7:91875/1‑149 (MQ=255)
                                                                                                                                       aCCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGcc        <  3:200509/149‑1 (MQ=255)
                                                                                                                                         cGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGa      <  5:90496/149‑1 (MQ=255)
                                                                                                                                           gCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCGGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGc    <  7:276400/149‑1 (MQ=255)
                                                                                                                                            cTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGcc   <  1:366627/149‑1 (MQ=255)
                                                                                                                                             tttGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGccc  <  8:170280/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTGCGAAGAAGGAGACGCGATGGGTGGCACCGGCTTTGCCGGTGTTCGCGGGTGAACCCGCTCCTACAATGGGCTGCGTTAGGCAGAAGCTTATTTTCGGCAGAAAAAAAAGCGCCGGGCGGTGGCCGCCGGGCGCTTTTAGGGGCTGACTGGGATCAGTCGCGGTACAGGATGGCCGAGCCC  >  NC_002947/4970375‑4970663

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: