Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 680,671 C→T 6.5% P64L (CCG→CTG)  PP_0583 → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947680,6710CT 6.5% 227.9 / 3.6 77P64L (CCG→CTG) PP_0583AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (35/37);  new base T (3/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CTGCCCAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTGG  >  NC_002947/680526‑680814
                                                                                                                                                 |                                                                                                                                                 
cTGCCCAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCg                                                                                                                                                  >  5:226479/1‑147 (MQ=255)
 tGCCCAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑Ag                                                                                                                                              >  6:355510/1‑149 (MQ=255)
         ccccTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGgt                                                                                                                                      <  1:78030/149‑1 (MQ=255)
          cccTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCCCTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGgtg                                                                                                                                     >  1:71410/1‑149 (MQ=255)
          cccTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACTGC‑GGTggggcaaa                                                                                                                                     <  6:305498/149‑5 (MQ=255)
           ccTCTACGCCCGTGCCGAAAGCCTGGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGgtgt                                                                                                                                    <  7:202703/149‑1 (MQ=255)
            ctctACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGtgtg                                                                                                                                   >  6:116316/1‑149 (MQ=255)
            ctctACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGtgtg                                                                                                                                   >  3:235190/1‑149 (MQ=255)
               tACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCCCTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTgg                                                                                                                                  <  2:71410/147‑1 (MQ=255)
                     cGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCccgc                                                                                                                          >  7:82931/1‑149 (MQ=255)
                      gTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCccgcc                                                                                                                         >  1:214450/1‑149 (MQ=255)
                          cGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATg                                                                                                                     >  6:331834/1‑149 (MQ=255)
                               gCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGAGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAg                                                                                                                >  3:57031/1‑149 (MQ=255)
                                ccTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACTGC‑GGTggggca                                                                                                                                       >  8:186351/1‑123 (MQ=255)
                                ccTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACTGC‑GGTggggca                                                                                                                                       <  7:186351/125‑3 (MQ=255)
                                 cTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGgtgt                                                                                                                                    <  6:150248/127‑1 (MQ=255)
                                 cTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGgtgt                                                                                                                                    >  5:150248/1‑127 (MQ=255)
                                     gcgcCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCTTGaa                                                                                                          <  6:226479/149‑1 (MQ=255)
                                     gcgcCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGaa                                                                                                          <  4:235190/149‑1 (MQ=255)
                                      cgcCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAg                                                                                                         <  3:163310/149‑1 (MQ=255)
                                       gcCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGt                                                                                                        >  4:240839/1‑149 (MQ=255)
                                        ccGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTg                                                                                                       <  4:166996/149‑1 (MQ=255)
                                               gTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCAcc                                                                                                 <  8:82931/148‑1 (MQ=255)
                                                tGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGaa                                                                                                          >  2:241563/1‑138 (MQ=255)
                                                tGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGaa                                                                                                          <  1:241563/138‑1 (MQ=255)
                                                tGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTc                                                                                               >  7:77978/1‑149 (MQ=255)
                                                 ggACAACTCGCCACCAGCACGACTGGGTCCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTTCCACACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTc                                                                                               >  7:143230/1‑148 (MQ=255)
                                                   aCCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGAt                                                                                            <  7:734/149‑1 (MQ=255)
                                                     cACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATg                                                                                           >  1:91302/1‑148 (MQ=255)
                                                       cTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑Caa                                                                                       <  7:299100/149‑1 (MQ=255)
                                                           ccaccaGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCg                                                                                   <  7:318448/149‑1 (MQ=255)
                                                            caccaGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCga                                                                                  <  5:331834/149‑1 (MQ=255)
                                                             accaGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACTGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCgag                                                                                 <  5:195102/149‑1 (MQ=255)
                                                                 gCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATgc                                                                             <  2:214450/149‑1 (MQ=255)
                                                                   aCGACTGGGTGCAATTCTCCTACGCCACCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATgcgc                                                                           <  2:72033/149‑1 (MQ=255)
                                                                       cTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGcc                                                                       >  3:141556/1‑149 (MQ=255)
                                                                         gggTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTAGACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTg                                                                     >  1:69539/1‑149 (MQ=255)
                                                                         gggTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTg                                                                     >  5:104502/1‑149 (MQ=255)
                                                                          ggTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTg                                                                     >  3:1984/1‑148 (MQ=255)
                                                                            tGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTAc                                                                  <  4:141556/149‑1 (MQ=255)
                                                                               aaTTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGAt                                                                                            >  6:83557/1‑121 (MQ=255)
                                                                               aaTTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGAt                                                                                            <  5:83557/121‑1 (MQ=255)
                                                                               aaTTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTg                                                               >  4:132302/1‑149 (MQ=255)
                                                                                    tcCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGc                                                                                                               >  6:342040/1‑97 (MQ=255)
                                                                                    tcCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGc                                                                                                               <  5:342040/97‑1 (MQ=255)
                                                                                     ccTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGcc                                                                       <  6:2658/135‑1 (MQ=255)
                                                                                     ccTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGACGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGcc                                                                       >  5:2658/1‑135 (MQ=255)
                                                                                      cTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGcgccgc                                                         >  4:103712/1‑148 (MQ=255)
                                                                                       tACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAg                                                                                                         <  4:316298/100‑1 (MQ=255)
                                                                                       tACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAg                                                                                                         >  3:316298/1‑100 (MQ=255)
                                                                                          gCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGAtgc                                                    >  6:70462/1‑149 (MQ=255)
                                                                                                    gcgTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATGCCAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATg                                           <  1:259953/149‑1 (MQ=255)
                                                                                                     cgTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATggg                                         >  2:179689/1‑149 (MQ=255)
                                                                                                          tCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCg                                    <  3:240839/149‑1 (MQ=255)
                                                                                                            ggCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCg                                    <  1:179689/147‑1 (MQ=255)
                                                                                                             gCGTGTACACACGCGACGGCAGTTACTTTGCCGCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGaa                                 <  8:143230/149‑1 (MQ=255)
                                                                                                             gCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGaa                                 <  2:91302/149‑1 (MQ=255)
                                                                                                                   acacacGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGc                           >  6:68116/1‑149 (MQ=255)
                                                                                                                    cacacGCGACGGCATTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGTATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGc                           <  5:355510/148‑1 (MQ=255)
                                                                                                                      cacGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCgg                        <  5:116316/149‑1 (MQ=255)
                                                                                                                         gcgACGGCAGTTACTTTGCCCCACCGCGAGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGa                                                       >  1:15973/1‑116 (MQ=255)
                                                                                                                         gcgACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGa                                                       <  2:15973/115‑1 (MQ=255)
                                                                                                                                cAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTg                     <  3:191145/142‑1 (MQ=255)
                                                                                                                                cAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTg                     >  4:191145/1‑142 (MQ=255)
                                                                                                                                 aGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGt             <  6:104502/149‑1 (MQ=255)
                                                                                                                                 aGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGt             <  2:69539/149‑1 (MQ=255)
                                                                                                                                  gTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCgg                        <  8:276959/137‑1 (MQ=255)
                                                                                                                                  gTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCgg                        >  7:276959/1‑137 (MQ=255)
                                                                                                                                     aCTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGaa                                 <  2:173673/125‑1 (MQ=255)
                                                                                                                                     aCTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGaa                                 >  1:173673/1‑125 (MQ=255)
                                                                                                                                      acaccccccACTGC‑GGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGAcgc        >  1:359136/6‑149 (MQ=255)
                                                                                                                                          cccccACTGC‑GGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACTTGCGCCGCGATGCCTGCCCATGGGAGCCGGAACATTGCCGGGTGCTGGAAGTGACGCCGCt    >  7:207190/2‑149 (MQ=255)
                                                                                                                                           ccccACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTg   <  8:77978/149‑1 (MQ=255)
                                                                                                                                            cccACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTgg  >  4:157880/1‑149 (MQ=255)
                                                                                                                                            cccACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTgg  >  3:301012/1‑149 (MQ=255)
                                                                                                                                             ccACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCggg                       >  6:67726/1‑127 (MQ=255)
                                                                                                                                             ccACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCggg                       <  5:67726/127‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTGCCCAGACCCCTCTACGCCCGTGCCGAAAGCCTTGGCGCCGGATCGTGGACCACTCGCCACCAGCACGACTGGGTGCAATTCTCCTACGCCATCAGCGGCGTACTCGGCGTGTACACACGCGACGGCAGTTACTTTGCCCCACCGC‑AGTGGGGGGTGTGGATCCCCGCCGATGCCGAGCATGAAGTGGTCACCTCGATG‑CAAGCCGAGATGCGCAGCCTGTACGTGCGCCGCGATGCCTGCCCATGGGCGCCGGAACAATGCCGGGTGCTGGAAGTGACGCCGCTGG  >  NC_002947/680526‑680814

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