Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,896,127 T→G 9.4% Q221P (CAG→CCG)  PP_3439 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,896,1270TG 9.4% 170.2 / 10.2 64Q221P (CAG→CCG) PP_3439AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base G (2/4);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

TGCACGGCTCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGG  >  NC_002947/3895980‑3896269
                                                                                                                                                   |                                                                                                                                               
tGCACGGCTCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTgg                                                                                                                                               <  1:1219374/150‑1 (MQ=255)
 gCACGGCTCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTggt                                                                                                                                              >  1:269251/1‑149 (MQ=255)
 gCACGGCTCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGgc                                                                                                                                              <  2:942334/150‑1 (MQ=255)
     ggCTCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCGGGGCAACTTCAGCAATTGAGCGGCCCTGGCGCaa                                                                                                                                          <  1:1403726/150‑1 (MQ=255)
         ccTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTgc                                                                                                                                      <  2:350153/150‑1 (MQ=255)
         ccTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCagttgc                                                                                                                                      >  2:202078/1‑145 (MQ=255)
               ggAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGCCAgg                                                                                                                                <  2:1291390/150‑1 (MQ=255)
               ggAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGCCAgg                                                                                                                                <  2:1291464/150‑1 (MQ=255)
               ggAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGCCAgg                                                                                                                                <  1:202078/150‑1 (MQ=255)
                gAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGAAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCggg                                                                                                                               >  2:1047821/1‑150 (MQ=255)
                    gcCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAg                                                                                                                          >  1:454448/1‑150 (MQ=255)
                           cTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGgcg                                                                                                                                             >  2:483152/1‑125 (MQ=255)
                           cTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGgcg                                                                                                                                             <  1:483152/125‑1 (MQ=255)
                             aGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGTGCAACTGGCCCGGGCGCAGTTGcgccag                                                                                                                 <  2:269251/150‑6 (MQ=255)
                                      aTTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAg                                                                                                        >  1:1348753/1‑150 (MQ=255)
                                        ttCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCa                                                                                                      <  1:241777/150‑1 (MQ=255)
                                          cTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAgg                                                                                                    <  2:1061302/150‑1 (MQ=255)
                                           tGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCGTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCAGGCCAGTTGCACGCGCTGGTCTAGCAGGt                                                                                                   >  2:122923/1‑150 (MQ=255)
                                           tGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGt                                                                                                   >  1:546309/1‑150 (MQ=255)
                                               aTGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGcc                                                                                               >  2:1058417/1‑150 (MQ=255)
                                                                 cAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGaa                                                                             <  2:546309/150‑1 (MQ=255)
                                                                          tttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCa                                                                                                                           >  1:1163684/1‑95 (MQ=255)
                                                                          tttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCa                                                                                                                           <  2:1163684/95‑1 (MQ=255)
                                                                          tttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCg                                                                    <  2:812776/150‑1 (MQ=255)
                                                                          tttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCg                                                                    >  1:476086/1‑150 (MQ=255)
                                                                              taTGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAg                                                                >  2:1159161/2‑150 (MQ=255)
                                                                                 gCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCa                                                             <  1:1132010/150‑1 (MQ=255)
                                                                                 gCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCa                                                             >  1:1255161/1‑150 (MQ=255)
                                                                                   ccGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACg                                                           >  2:1067072/1‑150 (MQ=255)
                                                                                          aGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAgg                                                    <  2:454448/150‑1 (MQ=255)
                                                                                           ggTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGt                                                   <  1:1058417/150‑1 (MQ=255)
                                                                                           ggTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGt                                                   <  1:1159161/150‑1 (MQ=255)
                                                                                            gtgtGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTa                                                  >  2:1347443/1‑150 (MQ=255)
                                                                                            gtgtGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTAGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTa                                                  >  2:919563/1‑150 (MQ=255)
                                                                                             tgtgCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAgg                                                    >  2:518786/1‑147 (MQ=255)
                                                                                             tgtgCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAgg                                                    <  1:518786/147‑1 (MQ=255)
                                                                                                     aTCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCgaga                                         <  2:1072285/150‑1 (MQ=255)
                                                                                                      tCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATc                                                       <  1:72424/135‑1 (MQ=255)
                                                                                                      tCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATc                                                       >  2:72424/1‑135 (MQ=255)
                                                                                                         gCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGtt                                                                                                                        <  1:130388/67‑1 (MQ=255)
                                                                                                         gCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGtt                                                                                                                        >  2:130388/1‑67 (MQ=255)
                                                                                                         gCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAgg                                     >  2:257820/1‑150 (MQ=255)
                                                                                                         gCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGAAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAgg                                     <  2:341237/150‑1 (MQ=255)
                                                                                                             aCCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGACCCTGGCGCTCAAACCAGGGGCCCGGTTGCACACCCGGGACGAGCATATAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGc                                 >  1:1243554/1‑150 (MQ=255)
                                                                                                                cGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGc                                                                                                                      >  1:110731/1‑62 (MQ=255)
                                                                                                                cGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGc                                                                                                                      <  2:110731/62‑1 (MQ=255)
                                                                                                                     tcctGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCg                         >  2:209694/1‑150 (MQ=255)
                                                                                                                     tcctGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCg                         >  2:745866/1‑150 (MQ=255)
                                                                                                                          tgCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCa                       <  1:1255972/147‑1 (MQ=255)
                                                                                                                          tgCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCa                       >  2:1255972/1‑147 (MQ=255)
                                                                                                                          tgCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGAt                    >  2:740343/1‑150 (MQ=255)
                                                                                                                              aCTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTcc                >  2:25830/1‑150 (MQ=255)
                                                                                                                                  cAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGa            <  2:476086/150‑1 (MQ=255)
                                                                                                                                  cAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGa            <  1:257820/150‑1 (MQ=255)
                                                                                                                                    gCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAgg          >  1:790416/1‑150 (MQ=255)
                                                                                                                                     cAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAgg                                                    <  2:75115/107‑1 (MQ=255)
                                                                                                                                     cAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAgg                                                    >  1:75115/1‑107 (MQ=255)
                                                                                                                                      aaTTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCg                                                                    <  2:515893/90‑1 (MQ=255)
                                                                                                                                      aaTTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCg                                                                    >  1:515893/1‑90 (MQ=255)
                                                                                                                                           aGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCg   >  1:302438/1‑150 (MQ=255)
                                                                                                                                           aGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCg   >  1:441940/1‑150 (MQ=255)
                                                                                                                                            gCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCgg  <  1:209694/150‑1 (MQ=255)
                                                                                                                                               gCCCGGGCGCAGTTGCGCCAGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGcc                                   <  2:54654/114‑1 (MQ=255)
                                                                                                                                               gCCCGGGCGCAGTTGCGCCAGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGcc                                   >  1:54654/1‑114 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TGCACGGCTCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGG  >  NC_002947/3895980‑3896269

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