Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,577,660 C→T 9.1% intergenic (‑68/+71) gntK ← / ← gntR gluconate kinase 2/d‑gluconate inducible gluconate regulon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,577,6600CT 9.1% 545.3 / 32.7 177intergenic (‑68/+71)gntK/gntRgluconate kinase 2/d‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (78/82);  new base T (8/8);  total (87/90)
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.49e-01

Only 100 of 177 total aligned reads displayed.

CCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCACGCTCGCAAGTCGTGGCTCCATCACCTGACCAATGTCATGACCGTGGAAACCGGCAATCGCCATATCGTCAGG  >  NC_000913/3577391‑3577931
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                               
ccTGCAACACCGGGTGGCGATCGTAGTCATTCAGTGGTTCGCCACACGCCAGTTTTTCGATAGTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGAc                                                                                                                                                                                                                                                                             <  2:472960/274‑1 (MQ=255)
                 cGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATatc                                                                                                                                                                                                                                                           <  2:313769/275‑1 (MQ=255)
                     cgtcgtCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATCTCAGAAGCAATatccatc                                                                                                                                                                                                                                                       >  2:416270/1‑275 (MQ=255)
                            aTTCAGGGGGTCGCCAGACGCCATTTCTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATGAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGcc                                                                                                                                                                                                                                                 <  2:115871/274‑1 (MQ=255)
                             ttCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCtt                                                                                                                                                                                                                                               <  2:169908/275‑1 (MQ=255)
                                                              ttGCGCCGTGGATGGGGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgtg                                                                                                                                                                                                              <  4:272698/275‑1 (MQ=255)
                                                                                    cGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCg                                                                                                                                                                                          >  1:127292/1‑273 (MQ=255)
                                                                                         tCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCa                                                                                                                                                                                                                                                                      >  3:80191/1‑192 (MQ=255)
                                                                                             gAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAg                                                                                                                                                                                                                                                                  <  1:80191/192‑1 (MQ=255)
                                                                                                               gatcagcCACGTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACTGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACAtt                                                                                                                                                             <  4:177072/270‑1 (MQ=255)
                                                                                                                               ggCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACa                                                                                                                                                                                                                                                                            >  6:281957/1‑148 (MQ=255)
                                                                                                                                gCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACaa                                                                                                                                                                                                                                                                           <  5:281957/148‑1 (MQ=255)
                                                                                                                                gCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCgggttcaacaattc                                                                                                                                            >  4:394487/1‑263 (MQ=255)
                                                                                                                                 cgcnngCAGATTTGCCGCGGCCCGATACGCCCATGTAGACGGAAATGTGGTGANCAGGGGTAGTCGGACTCAAAGTGGTGCACCCACAATACAGTAAGTAAAAATGTTACGGGTGACTGTTATCGGTAACATTGTCCAGCCGCACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtg                                                                                                                                                                                                                <  4:545519/201‑1 (MQ=255)
                                                                                                                                           tttGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACaa                                                                                                                                                                                                                                                                           <  5:189421/137‑1 (MQ=255)
                                                                                                                                           tttGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACaa                                                                                                                                                                                                                                                                           >  6:189421/1‑137 (MQ=255)
                                                                                                                                                           aCGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTc                                                                                                                 <  2:102234/275‑1 (MQ=255)
                                                                                                                                                            cGCCAAGCAAGACGTAAATGGGGTGATCATGTTTAGTAGTGCTCAAATTGGTGCCCCCAAAATACAAGAATTAAAANTGTTACGGGTAAATGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGTAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTc                                                                                                                 <  3:127292/274‑1 (MQ=255)
                                                                                                                                                            ataaaagataGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCACCACAAAACAAAAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTc                                                                                                                 <  1:155589/265‑1 (MQ=255)
                                                                                                                                                              cccATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTaa                                                                                                                                                                                                                                             >  8:706710/1‑148 (MQ=255)
                                                                                                                                                                 aTCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACaa                                                                                                                                                                                                                                                                           >  7:647831/1‑115 (MQ=255)
                                                                                                                                                                  tCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAgtgt                                                                                                                                                                                                                                         <  7:706710/148‑1 (MQ=255)
                                                                                                                                                                     aGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTa                                                                                                                                                                                                                                              >  7:746980/1‑140 (MQ=255)
                                                                                                                                                                              tgtgGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCTGACAATATCAGAAGCAATATCCATCCATGCCTTAAg                                                                                                                                                                                                                                            >  6:896750/1‑133 (MQ=255)
                                                                                                                                                                              tgtgGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAg                                                                                                                                                                                                                                            <  5:896750/133‑1 (MQ=255)
                                                                                                                                                                               gtggtgATCATGGTTAGTCGTGCTCAAAGTGGGGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGcc                                                                                                                                                                                            <  1:12868/180‑1 (MQ=255)
                                                                                                                                                                                 ggtgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAat                                                                                                                                                                                                                                                                          <  8:737241/100‑1 (MQ=255)
                                                                                                                                                                                  gtgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACa                                                                                                                                                                                                                                                                            >  8:693102/1‑97 (MQ=255)
                                                                                                                                                                                  gtgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCAGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAAATTCGCCACGAATACGCGACACNNNGCGTTCAGCGCCAAAAcctttcctttct                                                                                           >  3:155589/1‑263 (MQ=255)
                                                                                                                                                                                   tgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGtgtg                                                                                                                                                                                                                                >  6:346639/1‑140 (MQ=255)
                                                                                                                                                                                   tgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGtgtg                                                                                                                                                                                                                                <  5:346639/140‑1 (MQ=255)
                                                                                                                                                                                     aTCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCa                                                                                                                                                                                                                      >  6:1026111/1‑148 (MQ=255)
                                                                                                                                                                                       cATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAata                                                                                                                                                                                                                                                                         <  6:623705/95‑1 (MQ=255)
                                                                                                                                                                                           gTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtg                                                                                                                                                                                                                <  5:1026111/148‑1 (MQ=255)
                                                                                                                                                                                            ttAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGtata                                                                                                                                                                                                                                      <  7:297191/125‑1 (MQ=255)
                                                                                                                                                                                            ttAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGtata                                                                                                                                                                                                                                      >  8:297191/1‑125 (MQ=255)
                                                                                                                                                                                            ttAGTCGTGCTCAAAGTGGTGCCCACACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGATAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgt                                                                                                                                                                                                               <  8:600909/148‑1 (MQ=255)
                                                                                                                                                                                               gTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCTTTCAGCGCCAATACTGCCCATCCGCCCACACGGCGTCa                                                                              >  4:484515/1‑274 (MQ=255)
                                                                                                                                                                                               gTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTCTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTAAGCGCCAATACTGCCCATCCGATCaacgggcgtaag                                                                             >  4:205702/1‑264 (MQ=255)
                                                                                                                                                                                                tCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGtgtg                                                                                                                                                                                                                                <  8:844458/127‑1 (MQ=255)
                                                                                                                                                                                                    gCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaa                                                                                                                                                                                                       <  7:700747/148‑1 (MQ=255)
                                                                                                                                                                                                    gCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCCTCAGCACg                                                                         >  2:478569/1‑274 (MQ=255)
                                                                                                                                                                                                           gtggtgCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATc                                                                                                                                                                                                >  6:581416/1‑148 (MQ=255)
                                                                                                                                                                                                                gCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTCCAAAtt                                                                                                                                                                                                                  >  8:127085/1‑125 (MQ=255)
                                                                                                                                                                                                                gCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTCCAAAtt                                                                                                                                                                                                                  <  7:127085/125‑1 (MQ=255)
                                                                                                                                                                                                                 cccccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTa                                                                                                                                                                                                                                              <  8:576916/96‑1 (MQ=255)
                                                                                                                                                                                                                 cccccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTa                                                                                                                                                                                                                                              >  7:576916/1‑96 (MQ=255)
                                                                                                                                                                                                                    ccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCaa                                                                                                                                                                                                                     <  5:510909/118‑1 (MQ=255)
                                                                                                                                                                                                                     cacaATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATg                                                                                                                                                                                                                                                   >  6:354034/1‑87 (MQ=255)
                                                                                                                                                                                                                      acaATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCCCCAGCAGGCTTTCAGCGCCAATCCTGCCCATCCGCTCACGCGGCCTCAGCCCGCTCCCAAGTCGTTGcccc                                                       >  4:471969/1‑274 (MQ=255)
                                                                                                                                                                                                                              aGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGtgtg                                                                                                                                                                                                                                <  7:612842/97‑1 (MQ=255)
                                                                                                                                                                                                                              aGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGt                                                                                                                                                                                >  7:1015599/1‑145 (MQ=255)
                                                                                                                                                                                                                              aGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGt                                                                                                                                                                                <  8:1015599/145‑1 (MQ=255)
                                                                                                                                                                                                                                  ttAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAgg                                                                                                                                                                                 <  5:815554/140‑1 (MQ=255)
                                                                                                                                                                                                                                  ttAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAgg                                                                                                                                                                                 >  6:815554/1‑140 (MQ=255)
                                                                                                                                                                                                                                  ttAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGAGTTCCGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCCCGCTAGCACGCCGTGGCTCCATACCATGACCa                                           >  4:240860/1‑274 (MQ=255)
                                                                                                                                                                                                                                    aaaaaTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCACGCCCGCAAGTCGTGGCTCCATCACCTGACCAATg                                        >  2:387251/1‑275 (MQ=255)
                                                                                                                                                                                                                                      aaaTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATatccatc                                                                                                                                                                                                                                                       >  7:1096705/1‑66 (MQ=255)
                                                                                                                                                                                                                                      aaaTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATatccatc                                                                                                                                                                                                                                                       <  8:1096705/66‑1 (MQ=255)
                                                                                                                                                                                                                                       aaTGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatg                                                                                                                                                                                                                                                                         >  7:918227/1‑46 (MQ=37)
                                                                                                                                                                                                                                       aaTGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatg                                                                                                                                                                                                                                                                         <  8:918227/47‑2 (MQ=37)
                                                                                                                                                                                                                                        aTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGcc                                                                                                                                                                                            >  5:685398/1‑123 (MQ=255)
                                                                                                                                                                                                                                           ttACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaa                                                                                                                                                                                                       <  8:791045/109‑1 (MQ=255)
                                                                                                                                                                                                                                                ggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCgg                                                                                                                                                                                         >  8:32596/1‑118 (MQ=255)
                                                                                                                                                                                                                                                ggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCgg                                                                                                                                                                                         <  7:32596/118‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAAt                                                                                                                                                                                                     <  8:631597/104‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTaa                                                                                                                                                                       >  8:149723/1‑134 (MQ=255)
                                                                                                                                                                                                                                                  tAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTaa                                                                                                                                                                       <  7:149723/134‑1 (MQ=255)
                                                                                                                                                                                                                                                       gTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgtg                                                                                                                                                                                                              <  6:563731/90‑1 (MQ=255)
                                                                                                                                                                                                                                                       gTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgtg                                                                                                                                                                                                              >  5:563731/1‑90 (MQ=255)
                                                                                                                                                                                                                                                        ttATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgtg                                                                                                                                                                                                              >  7:345349/1‑89 (MQ=255)
                                                                                                                                                                                                                                                         tATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTcac                                                                                                                                                  >  6:110593/1‑148 (MQ=255)
                                                                                                                                                                                                                                                         tATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATAcgcg                                                                                                                              >  4:329167/1‑168 (MQ=255)
                                                                                                                                                                                                                                                          aTCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCACGCTCGCAAGTCGTGGCGCCATCACCTGACCAATGTCATGACCGTTGAAACCGGCaa                  >  2:266856/1‑275 (MQ=255)
                                                                                                                                                                                                                                                            cGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAgtgt                                                                                                                                                                                                                                 <  7:828107/66‑1 (MQ=255)
                                                                                                                                                                                                                                                             ggTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCtt                                                                                                                                                                                                                                               <  2:558500/51‑1 (MQ=255)
                                                                                                                                                                                                                                                               tAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGa                                                                                                                                                                                                  <  5:9388/94‑1 (MQ=255)
                                                                                                                                                                                                                                                               tAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGa                                                                                                                                                                                                  >  6:9388/1‑94 (MQ=255)
                                                                                                                                                                                                                                                                aaCATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATc                                                                                                                                                                                                >  7:15970/1‑95 (MQ=255)
                                                                                                                                                                                                                                                                aaCATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATc                                                                                                                                                                                                <  8:15970/95‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaCATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCACGCTCGCCAGTCGTGCCCCCATCACCTGACCAATGTCCTGACCGTGGAAACCGGCAATCGc              >  4:440086/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                  cATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCaaa                                                                                                                                                                                                                    >  3:26217/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                   aTTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaa                                                                                                                                                                                                       <  5:925345/85‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aTTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaa                                                                                                                                                                                                       >  6:925345/1‑85 (MQ=255)
                                                                                                                                                                                                                                                                    ttGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACCTTTTCGGTGTCACAGAtt                                                                                                                                             <  6:233684/142‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ttGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGAtt                                                                                                                                             >  5:233684/1‑142 (MQ=255)
                                                                                                                                                                                                                                                                    ttGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTaaa                                                                                                                                                                      >  5:689788/1‑117 (MQ=255)
                                                                                                                                                                                                                                                                    ttGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTaaa                                                                                                                                                                      <  6:689788/117‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ttGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTa                                                                                                                                                                  <  5:195619/121‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ttGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTa                                                                                                                                                                  >  6:195619/1‑121 (MQ=255)
                                                                                                                                                                                                                                                                     tGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACg                                                                                                                                      <  5:581416/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtg                                                                                                                                                                                                                <  1:26217/73‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGa                                                                                                                                                                                                  >  7:889219/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                       tCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGa                                                                                                                                                                                                  <  6:150612/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGCGAGCTACTTCAAATTTGTGGGCTTAAATAGa                                                                                                                                                                                                  >  5:150612/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                        ccAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCa                                                                                                                                                                                                                      <  4:567732/65‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATa                                                                                                                                  <  5:110593/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCACGCTCGCAAGTCGTGGCTCCATCACCTGACCAATGTCATGACCGTGGAAACCGGCAATCGCCATATCGTCAgg  <  2:205703/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATc                                                                                                                                                                                                >  5:529510/1‑84 (MQ=255)
                                                                                                                                                                                                                                                                           gCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATc                                                                                                                                                                                                <  6:529510/84‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                               
CCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCACGCTCGCAAGTCGTGGCTCCATCACCTGACCAATGTCATGACCGTGGAAACCGGCAATCGCCATATCGTCAGG  >  NC_000913/3577391‑3577931

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