Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,902,732 T→A 6.3% K102M (AAG→ATG) ‡ PP_5175 ← HlyD family secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,902,7320TA 6.3% 178.9 / 2.7 63K102M (AAG→ATG) ‡PP_5175HlyD family secretion protein
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base A (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

CTTGAGTGCGCTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCT  >  NC_002947/5902584‑5902876
                                                                                                                                                    |                                                                                                                                                 
cTTGAGTGCGCTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCt                                                                                                                                                   >  5:71235/1‑149 (MQ=255)
 ttGAGTGCGCTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCtt                                                                                                                                                  >  4:288656/1‑149 (MQ=255)
 ttGAGTGCGCTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCtt                                                                                                                                                  >  4:85016/1‑149 (MQ=255)
    aGTGCGGTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTCGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTgg                                                                                                                                                <  7:274177/148‑1 (MQ=255)
      tGCGCTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGcc                                                                                                                                              >  1:310331/1‑148 (MQ=255)
       gcgcTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGcc                                                                                                                                              >  1:24328/1‑147 (MQ=255)
             tGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGa                                                                                                                                       >  4:29918/1‑148 (MQ=255)
             tGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGa                                                                                                                                       >  6:73152/1‑148 (MQ=255)
               aTTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACa                                                                                                                                     <  1:34847/148‑1 (MQ=255)
                   ctgcGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGCCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCTAATACCTCGGCCTTGGCCGAGTCGACATTGtt                                                                                                                                <  2:24328/149‑1 (MQ=255)
                    tgcGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTc                                                                                                                               >  5:49678/1‑149 (MQ=255)
                    tgcGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGCGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCAAGGCCGAg                                                                                                                                           <  7:228106/137‑1 (MQ=255)
                             cgGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCt                                                                                                                     <  3:72731/149‑1 (MQ=255)
                              ggCCGCGGCAGTAGGGTACTCGCTCGGGGTGGTGTAGCCCTTGGGCATTTTTTTTTGCGGGCGCACTATGGCCGCTCTGGTCTGGGTGCCCGGCCCCGGCGGGGCTAATACCTGGGCCTGGGCCGACTCGACATTGTTCTGCAGGTCtt                                                                                                                    <  6:49678/149‑1 (MQ=255)
                                ccTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTc                                                                                                                  <  3:85016/149‑1 (MQ=255)
                                 cTCGGCAGCATGGTACTCGCTCTGGCTGGTGTTCCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCg                                                                                                                 <  6:352892/149‑1 (MQ=255)
                                    ggCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCggg                                                                                                               >  7:281818/1‑148 (MQ=255)
                                       agcaTCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCAAGGCCGAg                                                                                                                                           <  2:346617/118‑1 (MQ=255)
                                       agcaTCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCAAGGCCGAg                                                                                                                                           >  1:346617/1‑118 (MQ=255)
                                         caTCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTc                                                                                                             <  6:108463/145‑1 (MQ=255)
                                         caTCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTc                                                                                                             >  5:108463/1‑145 (MQ=255)
                                           tactACTCTCTCTCCCTGCTGTAGCCCTTGGGCAGCACATTTTGCTCGCGCAAGAAGCCCGCACTGGTCTGGGTGACCCGCCCCTGCGCGTCGAATAACTCGGCCTTGGCAGAGTCGACAATCTTATGCAGATCTATCCGGTCAAGccc                                                                                                       >  6:194704/4‑148 (MQ=255)
                                              tACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGc                                                                                                    <  6:111819/149‑1 (MQ=255)
                                                 tCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAg                                                                                                 >  2:109586/1‑149 (MQ=255)
                                                 tCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACACGCGCCTGCGCGGCGAATACCTCGGCCTTGGCAGAGTCGACATTGTTCTTCAGGTCTTTCGGGTAAAGACGGGc                                                                                                    >  7:78130/1‑146 (MQ=255)
                                                   gCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCa                                                                                               >  4:155385/1‑149 (MQ=255)
                                                   gCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCa                                                                                               >  8:1333/1‑149 (MQ=255)
                                                   gCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACATCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCa                                                                                               >  6:2643/1‑149 (MQ=255)
                                                         gctggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGCCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCCGGCGAGCACCTGGt                                                                                         >  5:237414/1‑149 (MQ=255)
                                                           tggTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGAGATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTg                                                                                       >  6:334759/1‑149 (MQ=255)
                                                             gtgtAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGCGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTgg                                                                                      <  1:235259/149‑1 (MQ=255)
                                                                   ccgttGGGCAGCTGTTTTTGGTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGAc                                                                               <  6:71235/146‑1 (MQ=255)
                                                                        gggCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGAt                                                                          >  3:353085/1‑149 (MQ=255)
                                                                         ggCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGAt                                                                          <  3:155385/148‑1 (MQ=255)
                                                                          gcagcaGTTTTTGCTGGCGCACGAAGGCGGCACTGGGCTGGGTGACCCGCGCCTGCGCGGCTAATACCTCGGCATTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGAc                                                                               <  4:252842/142‑1 (MQ=255)
                                                                          gaagcagTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGAGGCGAATACCTCGGCCTTTGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGATATCg                                                                        >  3:252842/3‑149 (MQ=255)
                                                                             gcagTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGccc                                                                     <  5:73152/149‑1 (MQ=255)
                                                                                     tGCTGGCGCACGTAGGCGGCACTGGTCTTGGTGACCGGCGCCTTCGCGGCTAATACCTCTGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGTGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGGGCGATCACGTcaa                                                              <  8:78130/148‑3 (MQ=255)
                                                                                            gcACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTAcgc                                                      >  5:137056/1‑149 (MQ=255)
                                                                                            gcACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTAcgc                                                      >  1:85294/1‑149 (MQ=255)
                                                                                             cACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTAcgc                                                      >  6:167464/1‑148 (MQ=255)
                                                                                             cACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGACGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGc                                                                                                    >  3:111237/1‑102 (MQ=255)
                                                                                             cACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGACTCGACATTGTTCTGCAGGTCTTTCGGGTCAACCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCAAGTCGACGCTAcgc                                                      >  3:337508/1‑148 (MQ=255)
                                                                                                gAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGaaa                                                  <  7:1333/149‑1 (MQ=255)
                                                                                                 aaGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAAt                                                 <  5:167464/149‑1 (MQ=255)
                                                                                                  aGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAAt                                                 >  4:219505/1‑148 (MQ=255)
                                                                                                       gcgtgggTCTGGTGGCCCGGGGCCTGCGGGGCGAAGCCCTGGGCGGTGGCGGCGTCGACATTGTGCTGCAGTGCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCtt                                           <  6:41463/144‑1 (MQ=255)
                                                                                                         aCTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCGTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTgc                                         <  8:281818/149‑1 (MQ=255)
                                                                                                         aCTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTgc                                         >  3:222937/1‑149 (MQ=255)
                                                                                                             gTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTgccgc                                      <  2:310331/148‑1 (MQ=255)
                                                                                                                tGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGAc                                  >  7:224372/1‑149 (MQ=255)
                                                                                                                    cgccccGCGCCTGCGCGGCTTTTACGGGGGCCTTGGCCGTTTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCAGCTGGTTGGGCTTTACTTGATCGCGCACGTCTACGCTACGCGAAATTATCTTGCCGCCGACGCgg                              <  6:237414/146‑1 (MQ=255)
                                                                                                                      aCCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGaa                            >  2:208229/1‑149 (MQ=255)
                                                                                                                       cccGCGCCGGCGCGGCTAATACCTCGGCCTGGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAAt                           <  4:353085/149‑1 (MQ=255)
                                                                                                                        ccGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAAt                           >  7:258374/1‑148 (MQ=255)
                                                                                                                          gcgcCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGAc                        <  1:208229/149‑1 (MQ=255)
                                                                                                                           cgcCTGCGCGGCTAAACCCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACa                       <  8:224372/149‑1 (MQ=255)
                                                                                                                                 cggggccAATCCCTCGGCCTTGGCCTTGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATc                   <  4:337508/144‑1 (MQ=255)
                                                                                                                                  gcggcgAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCgg                 >  5:155568/1‑148 (MQ=255)
                                                                                                                                   cggcgAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTc               <  1:109586/149‑1 (MQ=255)
                                                                                                                                    ggcgAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCt              <  3:288656/149‑1 (MQ=255)
                                                                                                                                      cgAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGc            >  4:260793/1‑149 (MQ=255)
                                                                                                                                       gAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCa           >  3:264896/1‑149 (MQ=255)
                                                                                                                                             cTCGGCCAAGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGCAATGACAAATCGGTCTGCACCcgcg     >  5:216548/1‑149 (MQ=255)
                                                                                                                                                ggCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCTAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCTGAATGACAAATCGGTCTGCACCCgcgc    <  5:2643/147‑1 (MQ=255)
                                                                                                                                                ggCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCt  >  1:68112/1‑149 (MQ=255)
                                                                                                                                                    ttGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCt              >  6:5453/1‑133 (MQ=255)
                                                                                                                                                    ttGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCt              <  5:5453/133‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTTGAGTGCGCTCTGATTGCTGCGCAGCGCGGCCTCGGCAGCATCGTACTCGCTCTGGCTGGTGTAGCCCTTGGGCAGCAGTTTTTGCTGGCGCACGAAGGCCGCACTGGTCTGGGTGACCCGCGCCTGCGCGGCGAATACCTCGGCCTTGGCCGAGTCGACATTGTTCTGCAGGTCTTTCGGGTCAAGCCGGGCGAGCACCTGGTTGGCCTTGACGTGATCGCCCACGTCGACGCTACGCGAAATGATCTTGCCGCCGACGCGGAATGACAAATCGGTCTGCACCCGCGCCT  >  NC_002947/5902584‑5902876

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