Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,866,969 G→A 9.0% intergenic (‑48/‑159) PP_3415 ← / → gnuK LacI family transcriptional regulator/D‑gluconate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,866,9690GA 9.0% 243.2 / 10.7 87intergenic (‑48/‑159)PP_3415/gnuKLacI family transcriptional regulator/D‑gluconate kinase
Reads supporting (aligned to +/- strand):  ref base G (34/45);  new base A (3/5);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

AAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTCATCGGCGAGACAGCGCTATCTCAC  >  NC_002947/3866827‑3867112
                                                                                                                                               |                                                                                                                                               
aaGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCCGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCt                                                                                                                                          <  1:710662/150‑1 (MQ=255)
   gcgcgcGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGtc                                                                                                                                       <  2:1049740/150‑1 (MQ=255)
     gcgcgAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGtctc                                                                                                                                     >  1:1555097/1‑150 (MQ=255)
     gcgcgAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGACCGGGGTTGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGtctc                                                                                                                                     >  1:62122/1‑150 (MQ=255)
       gcgAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGCACGGCCACTAAGATCGCGCTGTCTCAg                                                                                                                                   <  2:1842871/150‑1 (MQ=255)
       gcgAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCa                                                                                                                                    >  1:939422/1‑149 (MQ=255)
       gcgAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCa                                                                                                                                    <  2:939422/149‑1 (MQ=255)
       gcgAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAg                                                                                                                                   <  2:599711/150‑1 (MQ=255)
       gcgAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAg                                                                                                                                   <  2:368029/150‑1 (MQ=255)
          aaGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGAt                                                                                                                                <  1:1990676/150‑1 (MQ=255)
            gCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGc                                                                                                                              <  1:1237978/150‑1 (MQ=255)
             cGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCa                                                                                                                             >  1:1591485/1‑150 (MQ=255)
             cGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCa                                                                                                                             <  2:1555097/150‑1 (MQ=255)
              ggTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAg                                                                                                                            >  2:708272/1‑150 (MQ=255)
               gTAATCGGGGAAGCCCCGGAAAGCATGGCAAGGTTCCACAGTGTAGGCCGCGTGGTAGTTCGGGAGCTGTGGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGATAGCGCTGTCTCAGGATGCAgg                                                                                                                           <  2:1728140/150‑1 (MQ=255)
               gTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAgcg                                                                                                                                            <  1:1649074/133‑1 (MQ=255)
               gTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAgcg                                                                                                                                            >  2:1649074/1‑133 (MQ=255)
                    cGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATg                                                                                                                      <  2:3502/150‑1 (MQ=255)
                     ggggAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGt                                                                                                                     >  1:213701/1‑150 (MQ=255)
                     ggggAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTACGCCCCTGACAAGATAGAGTTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGt                                                                                                                     >  2:695164/1‑150 (MQ=255)
                         aaaCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATc                                                                                                                 <  2:1056889/150‑1 (MQ=255)
                             cccGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTg                                                                                                             >  1:536969/1‑150 (MQ=255)
                             cccGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTAAGGATGCAGGCAATGTAATCTTTg                                                                                                             >  2:2287367/1‑150 (MQ=255)
                              ccGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTgg                                                                                                            <  1:787370/150‑1 (MQ=255)
                              ccGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTgg                                                                                                            <  2:536969/150‑1 (MQ=255)
                                ggAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTa                                                                                                          <  1:168869/150‑1 (MQ=255)
                                   aaGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGACAAGTAATCTTTGGTagca                                                                                                       <  1:695164/150‑1 (MQ=255)
                                        tGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCttt                                                                                                              <  1:259317/138‑1 (MQ=255)
                                        tGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCttt                                                                                                              >  2:259317/1‑138 (MQ=255)
                                          gCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGATAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGccacca                                                                                                >  2:2104291/1‑150 (MQ=255)
                                              cTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCAcgcg                                                                                            <  2:1591485/150‑1 (MQ=255)
                                              cTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCAcgcg                                                                                            >  1:1098177/1‑150 (MQ=255)
                                               ttCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTa                                                                                                          >  2:2006911/1‑135 (MQ=255)
                                               ttCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTa                                                                                                          <  1:2006911/135‑1 (MQ=255)
                                               ttCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGcc                                                                                                    >  1:989687/1‑141 (MQ=255)
                                               ttCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGcc                                                                                                    <  2:989687/141‑1 (MQ=255)
                                               ttCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACgcgc                                                                                           <  2:372121/150‑1 (MQ=255)
                                                   gCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGAAATGATGCGACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGc                                                                                       >  1:2100955/1‑150 (MQ=255)
                                                     cAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGcca                                                                                     <  2:1098177/150‑1 (MQ=255)
                                                           gggACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGACCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGc                                                                               <  1:1457618/150‑1 (MQ=255)
                                                            ggACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCg                                                                              >  2:740775/1‑150 (MQ=255)
                                                                gacCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGATAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCa                                                                          >  2:1615094/1‑150 (MQ=255)
                                                                   cTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTagca                                                                                                       <  2:2058507/118‑1 (MQ=255)
                                                                   cTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTagca                                                                                                       >  1:2058507/1‑118 (MQ=255)
                                                                   cTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATAc                                                                       <  1:708272/150‑1 (MQ=255)
                                                                    tGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACt                                                                      <  2:1628506/150‑1 (MQ=255)
                                                                        gTTCGGGGGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGt                                                                  >  2:1089726/1‑150 (MQ=255)
                                                                         ttCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTg                                                                                                             >  2:76520/1‑106 (MQ=255)
                                                                         ttCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTg                                                                                                             <  1:76520/106‑1 (MQ=255)
                                                                         ttCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTa                                                                 <  1:122238/150‑1 (MQ=255)
                                                                          tCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTaa                                                                <  2:610168/150‑1 (MQ=255)
                                                                          tCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTaa                                                                <  2:2100955/150‑1 (MQ=255)
                                                                             ggAGCCGGTGCGGGACATGAATGTTGTAATTTTAGTACTTGCAAGGATTGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTAGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGa                                                             >  2:2082511/1‑150 (MQ=255)
                                                                              gggCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGAc                                                            <  1:1089726/148‑1 (MQ=255)
                                                                              gAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGAc                                                            >  2:789081/1‑150 (MQ=255)
                                                                              gAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTAAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGAc                                                            <  1:2287367/150‑1 (MQ=255)
                                                                                gCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGATAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACaa                                                          <  2:1956899/150‑1 (MQ=255)
                                                                                       cGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAAt                                                   <  1:740775/150‑1 (MQ=255)
                                                                                            cATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAacac                                              <  2:213701/150‑1 (MQ=255)
                                                                                              tGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCg                                            <  1:1941606/150‑1 (MQ=255)
                                                                                               gAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAAt                                                                                                                       >  2:2132108/1‑74 (MQ=255)
                                                                                               gAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAAt                                                                                                                       <  1:2132108/74‑1 (MQ=255)
                                                                                                 atttttgtAATTTTACTACTTGCACGGAGAGGGAACGGCGACTACTGTCGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACACGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGTCCAATAACACCggg                                          <  2:916561/146‑1 (MQ=255)
                                                                                                   tgttgtAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTgg                                                                                                            <  1:646196/82‑1 (MQ=255)
                                                                                                   tgttgtAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTgg                                                                                                            >  2:646196/1‑82 (MQ=255)
                                                                                                           ttttACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGACTg                                >  2:1425394/1‑150 (MQ=255)
                                                                                                             ttACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGtt                              >  1:388666/1‑150 (MQ=255)
                                                                                                              tactacTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCCTTCCCCTTGTGGAGGTGGCTGTTc                             <  2:1419232/150‑1 (MQ=255)
                                                                                                              tactacTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTc                             <  2:1416415/150‑1 (MQ=255)
                                                                                                               actactTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCg                            >  1:1317104/1‑150 (MQ=255)
                                                                                                                 tactTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTc                          >  2:781069/1‑150 (MQ=255)
                                                                                                                    ttGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGt                               >  2:1876394/1‑142 (MQ=255)
                                                                                                                    ttGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGt                               <  1:1876394/142‑1 (MQ=255)
                                                                                                                    ttGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTCATc                       >  2:97862/1‑150 (MQ=255)
                                                                                                                    ttGCAAGGAGAGGGAACGGCCACTAAGATAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTCATc                       >  1:1848346/1‑150 (MQ=255)
                                                                                                                        aaGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTCATCGGCg                   >  1:696419/1‑150 (MQ=255)
                                                                                                                                ggAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTCATCGGCGAGACAgcg           >  2:2132256/1‑150 (MQ=255)
                                                                                                                                  aaCGGCCACTAAGATAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTCATCGGCGAGACAGCGCt         <  1:2104291/150‑1 (MQ=255)
                                                                                                                                     ggCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTCATCGGCGAGACAGCGCTAtc      <  2:1317104/150‑1 (MQ=255)
                                                                                                                                      gCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCACCACCAGGCGTc                                                                            >  1:2265825/1‑79 (MQ=255)
                                                                                                                                      gCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCACCACCAGGCGTc                                                                            <  2:2265825/79‑1 (MQ=255)
                                                                                                                                       ccACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTCATCGGCGAGACAGCGCTAtctc    >  1:1344724/1‑150 (MQ=255)
                                                                                                                                        cACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGa                                        >  2:2276935/1‑113 (MQ=255)
                                                                                                                                        cACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGa                                        <  1:2276935/113‑1 (MQ=255)
                                                                                                                                         aCTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTCATCGGCGAGACAGCGCTATCTCAc  <  1:486045/150‑1 (MQ=255)
                                                                                                                                         aCTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTCATCGGCGAGACAGCGCTATCTCAc  <  1:2132256/150‑1 (MQ=255)
                                                                                                                                         aCTAAGATAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTCATCGGCGAGACAGCGCTATCTCAc  <  2:1848346/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
AAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAATCTTTGGTAGCAGCCACCACGCGCTCGCCACCAGGCGTCCATACTGCGTAAGGACAAGACCAATAACACCGGGGAGGTGGCTGTTCGTCATCGGCGAGACAGCGCTATCTCAC  >  NC_002947/3866827‑3867112

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: