Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,577,657 A→G 6.2% intergenic (‑65/+74) 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,6570AG 6.2% 454.2 / 12.6 142intergenic (‑65/+74)gntK/gntRgluconate kinase 2/d‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (59/74);  new base G (4/5);  total (63/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 145 total aligned reads displayed.

ACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACG  >  NC_000913/3577510‑3577805
                                                                                                                                                   |                                                                                                                                                    
aCTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAg                                                                                                                                                     >  2:860010/1‑149 (MQ=255)
aCTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAg                                                                                                                                                     >  3:1020169/1‑149 (MQ=255)
   tCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGctg                                                                                                                                                  <  3:238827/149‑3 (MQ=255)
         gCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAat                                                                                                                                            <  4:696123/149‑1 (MQ=255)
         gCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAat                                                                                                                                            >  2:878173/1‑149 (MQ=255)
             ccGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatgtta                                                                                                                                        <  1:489374/149‑5 (MQ=255)
             ccGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCa                                                                                                                                        <  1:236479/149‑1 (MQ=255)
              cGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAg                                                                                                                                       >  1:955392/1‑149 (MQ=255)
               gCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGa                                                                                                                                      >  3:552633/1‑149 (MQ=255)
                 aGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAg                                                                                                                                    <  3:725289/149‑1 (MQ=255)
                  gATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGc                                                                                                                                   <  2:969081/149‑1 (MQ=255)
                       gCCGCTGCACGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAGTTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAatat                                                                                                                              <  4:603507/149‑1 (MQ=255)
                          gCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATatcca                                                                                                                           >  2:390401/1‑149 (MQ=255)
                             gCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAata                                                                                                                               >  1:135627/1‑142 (MQ=255)
                             gCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAata                                                                                                                               <  2:135627/142‑1 (MQ=255)
                                 gATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGc                                                                                                                    >  3:759841/1‑149 (MQ=255)
                                  aTACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGcc                                                                                                                   <  4:151778/149‑1 (MQ=255)
                                   tACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatgttac                                                                                                                                       >  3:808462/1‑123 (MQ=255)
                                      gCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAAGATTGTCCAGCCGTACAATATCAGAAGCAATATCCATCCATGCCTTaa                                                                                                               <  1:762134/149‑1 (MQ=255)
                                        ccATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAgt                                                                                                             >  1:635138/1‑149 (MQ=255)
                                        ccATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAgt                                                                                                             <  4:20432/149‑1 (MQ=255)
                                         cATCAAGACGTAAATCTGGTGCTCATGGTTAGTCGTGCTCAAAGTGGTGCCCGCACAATACAAGAATTAAAAAGGTTACGGGTAACTGTTATCGGTAACTTTGTCCAGCCGGACAATATCAGAAGCAATACCCATCCATGCTTTAAgtg                                                                                                            <  2:787517/149‑1 (MQ=255)
                                          aTCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTaa                                                                                                               >  4:582856/1‑145 (MQ=255)
                                           tCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGc                                                                                                                    <  4:714829/139‑1 (MQ=255)
                                           tCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGta                                                                                                          <  3:185011/149‑1 (MQ=255)
                                           tCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCAAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGta                                                                                                          <  3:923199/149‑1 (MQ=255)
                                              aGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATaa                                                                                                       <  1:218801/149‑1 (MQ=255)
                                              aGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATaa                                                                                                       <  4:765692/149‑1 (MQ=255)
                                              aGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATaa                                                                                                       <  3:61941/149‑1 (MQ=255)
                                              aGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATaa                                                                                                       <  1:858789/149‑1 (MQ=255)
                                               gACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAg                                                                                                      <  3:165714/149‑1 (MQ=255)
                                               gACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAg                                                                                                      >  2:1051825/1‑149 (MQ=255)
                                                  gTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGc                                                                                                                                                    <  1:524925/100‑1 (MQ=255)
                                                  gTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGc                                                                                                                                                    >  2:524925/1‑100 (MQ=255)
                                                  gTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAgtgt                                                                                                   >  2:25926/1‑149 (MQ=255)
                                                  gTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAgtgt                                                                                                   <  2:516788/149‑1 (MQ=255)
                                                   tAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGtgtg                                                                                                  >  3:561401/1‑149 (MQ=255)
                                                    aaaTGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGa                                                                                                 <  2:324178/149‑1 (MQ=255)
                                                    aaaTGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGa                                                                                                 <  3:50191/149‑1 (MQ=255)
                                                    aaaTGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGa                                                                                                 <  4:52846/149‑1 (MQ=255)
                                                      aTGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCGTCCATGCCTTAAGTGTATAAGTGTGAGc                                                                                               <  1:74020/149‑1 (MQ=255)
                                                         tggtgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTAc                                                                                            <  3:871765/149‑1 (MQ=255)
                                                            tgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAat                                                                                                                                            <  3:275000/98‑1 (MQ=255)
                                                            tgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAat                                                                                                                                            >  4:275000/1‑98 (MQ=255)
                                                            tgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTc                                                                                         >  2:611772/1‑149 (MQ=255)
                                                               tCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCaaa                                                                                      <  2:344618/149‑1 (MQ=255)
                                                                    gTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgt                                                                                 >  1:136092/1‑149 (MQ=255)
                                                                    gTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgt                                                                                 <  1:878173/149‑1 (MQ=255)
                                                                     ttAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgtg                                                                                <  4:552633/149‑1 (MQ=255)
                                                                      tAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTgg                                                                               <  2:432237/149‑1 (MQ=255)
                                                                      tAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTgg                                                                               <  2:136092/149‑1 (MQ=255)
                                                                             gCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCAGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaaa                                                                        <  3:916728/149‑1 (MQ=255)
                                                                               tCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCaa                                                                                                                                 >  4:447527/1‑90 (MQ=255)
                                                                               tCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATa                                                                      >  1:258732/1‑149 (MQ=255)
                                                                               tCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATa                                                                      >  2:194809/1‑149 (MQ=255)
                                                                                 aaaGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgt                                                                                 <  3:41116/136‑1 (MQ=255)
                                                                                  aaGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGAt                                                                   <  4:826932/149‑1 (MQ=255)
                                                                                  aaGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGAt                                                                   <  1:57084/149‑1 (MQ=255)
                                                                                    gtggtgCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATccatccat                                                                                                                      <  3:38965/96‑1 (MQ=255)
                                                                                       gtgCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGcc                                                              >  1:833138/1‑149 (MQ=255)
                                                                                            cccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTg                                                         >  4:471128/1‑149 (MQ=255)
                                                                                            cccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTg                                                         >  3:777251/1‑149 (MQ=255)
                                                                                             ccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGa                                                        >  2:692238/1‑149 (MQ=255)
                                                                                                   tACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGt                                                  <  1:343652/149‑1 (MQ=255)
                                                                                                    aCAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTg                                                 >  1:536072/1‑149 (MQ=255)
                                                                                                      aaGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGaa                                               >  2:782986/1‑149 (MQ=255)
                                                                                                       aGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGaaa                                              <  4:19380/149‑1 (MQ=255)
                                                                                                        gAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAAc                                             >  3:871418/1‑149 (MQ=255)
                                                                                                           ttAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTa                                          >  1:1036626/1‑149 (MQ=255)
                                                                                                           ttAAAAATGTTACGGGTAACTGTTATCGGCAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGaaa                                              >  4:138468/1‑145 (MQ=255)
                                                                                                              aaaaTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAAt                                       >  2:342116/1‑149 (MQ=255)
                                                                                                              aaaaTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGCGAAACCTAAAt                                       >  2:708963/1‑149 (MQ=255)
                                                                                                                aaTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCt                                     <  3:365238/149‑1 (MQ=255)
                                                                                                                  tGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTaa                                   <  4:824674/149‑1 (MQ=255)
                                                                                                                  tGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTaa                                   <  3:471128/149‑1 (MQ=255)
                                                                                                                  tGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTaa                                   <  1:686433/149‑1 (MQ=255)
                                                                                                                        gggTAACTGTTATCGGTAACATTGTTCAGCCGGACAATATCAGACGCAATATCCATCCGTGACTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAAACTAAATCTTCCAtttt                             <  4:345146/149‑1 (MQ=255)
                                                                                                                        gggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaaa                                                                        >  4:405937/1‑106 (MQ=255)
                                                                                                                           tAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCgg                          <  4:682425/149‑1 (MQ=255)
                                                                                                                            aaCTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGgt                         <  4:561401/149‑1 (MQ=255)
                                                                                                                                gTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTca                     <  3:437748/149‑1 (MQ=255)
                                                                                                                                   aTCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAg                  >  1:984200/1‑149 (MQ=255)
                                                                                                                                    tCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGa                 <  1:342116/149‑1 (MQ=255)
                                                                                                                                       gTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGt                                                          <  3:655154/105‑1 (MQ=255)
                                                                                                                                       gTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGt                                                          >  4:655154/1‑105 (MQ=255)
                                                                                                                                       gTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTc              >  3:559431/1‑149 (MQ=255)
                                                                                                                                       gTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTc              <  1:692238/149‑1 (MQ=255)
                                                                                                                                       gTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTc              <  2:3410/149‑1 (MQ=255)
                                                                                                                                       gTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTc              <  2:984200/149‑1 (MQ=255)
                                                                                                                                        tgataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCg             >  2:591426/5‑149 (MQ=255)
                                                                                                                                        tAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCg             >  4:372278/1‑149 (MQ=255)
                                                                                                                                          ataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTcaca                   >  3:805993/3‑141 (MQ=255)
                                                                                                                                          ataTTGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTcaca                   <  4:805993/139‑1 (MQ=255)
                                                                                                                                          aCATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTaa                                         <  1:938215/119‑1 (MQ=255)
                                                                                                                                          aCATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTaa                                         >  2:938215/1‑119 (MQ=255)
                                                                                                                                           cATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCa          >  1:17114/1‑149 (MQ=255)
                                                                                                                                            aTTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCAc         >  2:955302/1‑149 (MQ=255)
                                                                                                                                            aTTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCAc         >  2:689649/1‑149 (MQ=255)
                                                                                                                                                   aGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATAcg  >  2:426443/1‑149 (MQ=255)
                                                                                                                                                   aGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATAcg  >  4:398996/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACG  >  NC_000913/3577510‑3577805

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