Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,577,657 A→G 6.0% intergenic (‑65/+74) gntK ← / ← gntR D‑gluconate kinase, thermostable/DNA‑binding transcriptional repressor GntR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,577,6570AG 6.0% 303.5 / 6.7 100intergenic (‑65/+74)gntK/gntRD‑gluconate kinase, thermostable/DNA‑binding transcriptional repressor GntR
Reads supporting (aligned to +/- strand):  ref base A (48/46);  new base G (2/4);  total (50/50)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACG  >  NC_000913/3577511‑3577805
                                                                                                                                                  |                                                                                                                                                    
cTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGTTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGcc                                                                                                                                                   <  1:505080/150‑1 (MQ=255)
  tCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCgg                                                                                                                                                 <  2:806909/150‑1 (MQ=255)
  tCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCgg                                                                                                                                                 >  1:1151045/1‑150 (MQ=255)
       ggCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAat                                                                                                                                            >  2:150377/1‑150 (MQ=255)
        gCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatg                                                                                                                                           <  2:593412/150‑2 (MQ=255)
        gCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAata                                                                                                                                           <  1:75891/150‑1 (MQ=255)
         cGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGCGGTGATCATTGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAatat                                                                                                                                          <  1:809777/150‑1 (MQ=255)
          gACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATc                                                                                                                                         <  2:398365/150‑1 (MQ=255)
            ccGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatgttac                                                                                                                                       <  1:1616278/150‑6 (MQ=255)
            ccGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAg                                                                                                                                       >  1:886962/1‑150 (MQ=255)
            ccGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAg                                                                                                                                       >  2:1110822/1‑150 (MQ=255)
             cGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGa                                                                                                                                      >  1:525839/1‑150 (MQ=255)
                aGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGc                                                                                                                                   <  2:1551844/150‑1 (MQ=255)
                 taTTTGCCTCTGCCCGATAAGCCCATCATGACGTAAATGTCGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCa                                                                                                                                  <  2:548743/149‑1 (MQ=255)
                  aTTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCaa                                                                                                                                 <  1:150377/150‑1 (MQ=255)
                     tGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAatat                                                                                                                              <  1:49473/150‑1 (MQ=255)
                       ccGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTATTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATatcc                                                                                                                            <  1:1541324/150‑1 (MQ=255)
                           tGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATatccatcc                                                                                                                        <  2:702514/150‑1 (MQ=255)
                              ccGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATg                                                                                                                     >  2:1180955/1‑150 (MQ=255)
                              ccGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATg                                                                                                                     >  2:156537/1‑150 (MQ=255)
                               cGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGc                                                                                                                    >  2:1013929/1‑150 (MQ=255)
                                          tCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGtat                                                                                                         >  1:1365434/1‑150 (MQ=255)
                                            aaGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATaa                                                                                                       >  2:1479776/1‑150 (MQ=255)
                                              gACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCACCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCAAGCCTTAAGTGTATAAgt                                                                                                     >  1:579875/1‑150 (MQ=255)
                                              caCGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAgt                                                                                                     <  1:299448/149‑1 (MQ=255)
                                               aCGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAgtg                                                                                                    >  1:565936/1‑150 (MQ=255)
                                               aCGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAgtg                                                                                                    <  2:1319004/150‑1 (MQ=255)
                                                cGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAgtgt                                                                                                   >  1:369295/1‑150 (MQ=255)
                                                    aaTGTGGTGATCATGGTTAGTCTTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGc                                                                                               <  2:1180991/150‑1 (MQ=255)
                                                     aTGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAACACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCt                                                                                              <  1:420252/150‑1 (MQ=255)
                                                           tgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCa                                                                                        >  1:877893/1‑150 (MQ=255)
                                                           tgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCa                                                                                        >  2:181458/1‑150 (MQ=255)
                                                            gATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATCGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCaa                                                                                       >  2:987893/1‑150 (MQ=255)
                                                                 tGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtg                                                                                  <  1:661084/150‑1 (MQ=255)
                                                                 tGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtg                                                                                  >  1:1529456/1‑150 (MQ=255)
                                                                   gTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgtg                                                                                <  1:523935/150‑1 (MQ=255)
                                                                   gTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgtg                                                                                >  2:1277059/1‑150 (MQ=255)
                                                                    ttAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTgg                                                                               <  1:129250/150‑1 (MQ=255)
                                                                    ttAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTgg                                                                               <  2:565936/150‑1 (MQ=255)
                                                                         cGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTa                                                                          >  1:909896/1‑150 (MQ=255)
                                                                          gTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaa                                                                         >  2:1252481/1‑150 (MQ=255)
                                                                           tGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaaa                                                                        <  2:1246243/150‑1 (MQ=255)
                                                                            gCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAAt                                                                       >  1:1330739/1‑150 (MQ=255)
                                                                             cTCAAAGTGGTGCCCCCACAATCCAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATa                                                                      >  1:530716/1‑150 (MQ=255)
                                                                             cTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATa                                                                      <  1:307004/150‑1 (MQ=255)
                                                                               cAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGa                                                                    >  2:625048/1‑150 (MQ=255)
                                                                               cAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGa                                                                    >  1:1276305/1‑150 (MQ=255)
                                                                                aaaGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGAt                                                                   <  1:131337/150‑1 (MQ=255)
                                                                                 aaGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATc                                                                  <  1:705278/150‑1 (MQ=255)
                                                                                 aaGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATc                                                                  <  2:1415483/150‑1 (MQ=255)
                                                                                 aaGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATc                                                                  <  1:1014556/150‑1 (MQ=255)
                                                                                 aaGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCAAGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATATGTGGGCTTAAATAGATc                                                                  <  2:579875/150‑1 (MQ=255)
                                                                                       tgCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCg                                                            >  2:1203481/1‑150 (MQ=255)
                                                                                         cccccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGt                                                          >  2:25103/1‑150 (MQ=255)
                                                                                         cccccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGt                                                          >  2:1095198/1‑150 (MQ=255)
                                                                                          ccccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTTAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTg                                                         >  1:859645/1‑150 (MQ=255)
                                                                                          ccccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTg                                                         >  1:1042624/1‑150 (MQ=255)
                                                                                          ccccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTg                                                         <  1:300515/150‑1 (MQ=255)
                                                                                           cccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGa                                                        <  2:1507755/150‑1 (MQ=255)
                                                                                           cccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGa                                                        <  1:525186/150‑1 (MQ=255)
                                                                                           cccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGa                                                        <  1:1110822/150‑1 (MQ=255)
                                                                                             cacaATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACa                                                      >  2:470235/1‑150 (MQ=255)
                                                                                              acaATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACaa                                                     <  1:1188643/150‑1 (MQ=255)
                                                                                                aaTACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAgg                                                   <  2:877893/150‑1 (MQ=255)
                                                                                                  tACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTg                                                 <  2:525839/150‑1 (MQ=255)
                                                                                                  tACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTg                                                 >  1:420351/1‑150 (MQ=255)
                                                                                                    cAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGaa                                               <  1:1180955/150‑1 (MQ=255)
                                                                                                    cAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGaa                                               >  2:49871/1‑150 (MQ=255)
                                                                                                      aGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTTAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAAc                                             <  2:859645/150‑1 (MQ=255)
                                                                                                      aGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAAc                                             <  2:464591/150‑1 (MQ=255)
                                                                                                       gAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAAcc                                            <  1:1203481/150‑1 (MQ=255)
                                                                                                        aaTTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCt                                           >  2:1336465/1‑150 (MQ=255)
                                                                                                           tAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCGATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTaaa                                        >  1:126447/1‑150 (MQ=255)
                                                                                                             aaaaTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATc                                      <  2:1529456/150‑1 (MQ=255)
                                                                                                             aaaaTGTTACGGGTAACTGTTATCGGTAACATTGTACAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAACCCTaaagc                                      >  2:855461/1‑148 (MQ=255)
                                                                                                                aTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTaa                                   <  2:244730/150‑1 (MQ=255)
                                                                                                                 tGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACaa                                                     >  2:773702/1‑131 (MQ=255)
                                                                                                                 tGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACaa                                                     <  1:773702/131‑1 (MQ=255)
                                                                                                                  gTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACa                                 >  1:1512642/1‑150 (MQ=255)
                                                                                                                    tACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACAtt                               >  2:1018058/1‑150 (MQ=255)
                                                                                                                        ggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTg                                                 >  1:788051/1‑128 (MQ=255)
                                                                                                                        ggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTg                                                 <  2:788051/128‑1 (MQ=255)
                                                                                                                            aCTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGgtgt                       >  2:1432886/1‑150 (MQ=255)
                                                                                                                             cTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTc                      >  2:1094715/1‑150 (MQ=255)
                                                                                                                                  aTCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGa                 >  2:248291/1‑150 (MQ=255)
                                                                                                                                  aTCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGATACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGa                 >  1:695307/1‑150 (MQ=255)
                                                                                                                                     tctgataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAATCCTAAATCTAAGATTTTCGGTGTCACAGATTc              >  1:917806/7‑150 (MQ=255)
                                                                                                                                     ggTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCCTCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTc              >  1:1184248/1‑150 (MQ=255)
                                                                                                                                      gTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCg             <  2:420351/150‑1 (MQ=255)
                                                                                                                                      ctgataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCg             <  1:423512/145‑1 (MQ=255)
                                                                                                                                       tgataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTcac                    >  1:731466/5‑142 (MQ=255)
                                                                                                                                       tgataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTcac                    <  2:731466/138‑1 (MQ=255)
                                                                                                                                        aaCATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGcc           >  1:1336945/1‑150 (MQ=255)
                                                                                                                                        aaCATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGcc           <  2:909896/150‑1 (MQ=255)
                                                                                                                                          cATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCAc         >  1:1092991/1‑150 (MQ=255)
                                                                                                                                           aTTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGGCACAGATTCGCCACg        <  2:1541727/150‑1 (MQ=255)
                                                                                                                                             tGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGaa      <  2:1124199/150‑1 (MQ=255)
                                                                                                                                             tGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGaa      <  1:1252481/150‑1 (MQ=255)
                                                                                                                                               tCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATa    >  1:1421354/1‑150 (MQ=255)
                                                                                                                                                 cAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATAcg  >  1:637937/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACG  >  NC_000913/3577511‑3577805

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