Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,117,047 G→T 62.3% S59Y (TCC→TAC)  glpK ← glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,117,0470GT62.3% 157.2 / 278.2 204S59Y (TCC→TAC) glpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base G (33/44);  new base T (64/63);  total (97/107)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

Only 100 of 204 total aligned reads displayed.

CACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGA  >  NC_000913/4116797‑4117303
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                
cacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGTCTTTCGCCAGCCCTTCTACCAGCTTGGAGCTTTGTGTGGCCCAGATTTCCa                                                                                                                                                                                                                                          >  1:15559/1‑275 (MQ=255)
                                ccAGTCCGGTTTTGCTGCGGATTTCATCTTCTAAACCGTCCCGGTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTggg                                                                                                                                                                                                          <  4:142206/275‑1 (MQ=255)
                                         tATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGTCAATGGCGTTATAGATAGGTTTGCCGGTATCCTTTTCCCAGACAATCGTGGTTTCACGCGGGTTCGTAGTACCTATATCTGCAATTTGATCGGAACTGCTACCGGCTTTCGCAAGCCATTCTACCAGCGTATAGCTTTGGGTCGCCCAGATTTCCATTGGTTCCGTTTCTACCCAACCTGGGTTTGGTAAGAtttcc                                                                                                                                                                                                 >  4:273884/1‑273 (MQ=255)
                                                          tcttcgAAACGTTAACGTTTTAAATCCTCGACGGTTTCTGCGGTCGGACGGCACTGCCAGACAATGGCGTTATTGTTAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGTTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGa                                                                                                                                                                                <  4:237253/275‑1 (MQ=255)
                                                                                                       cgacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                                                                                                                                                                                                                                                 <  5:938242/149‑1 (MQ=255)
                                                                                                              aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTgg                                                                                                                                                                                                                                                          <  5:734507/149‑1 (MQ=255)
                                                                                                              aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                                                                                                                                                                                                                                          <  5:526142/149‑1 (MQ=255)
                                                                                                               cTGCCAGACAATTGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTggg                                                                                                                                                                                                                                                         <  8:754829/149‑1 (MQ=255)
                                                                                                                 gCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTg                                                                                                                                                                                                                                                       <  8:315180/149‑1 (MQ=255)
                                                                                                                    aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGcc                                                                                                                                                                                                                                                    >  7:173565/1‑149 (MQ=255)
                                                                                                                              gTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCa                                                                                                                                                                                                                                          >  6:1133022/1‑149 (MQ=255)
                                                                                                                               ttATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGCTATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTccgt                                                                                                                                                                                                                                         >  8:1252642/1‑147 (MQ=255)
                                                                                                                               ttATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGt                                                                                                                                                                                                                                                        <  7:1059796/134‑1 (MQ=255)
                                                                                                                               ttATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGt                                                                                                                                                                                                                                                        >  8:1059796/1‑134 (MQ=255)
                                                                                                                                     atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                                                                                                                                                                                   >  5:1036219/1‑149 (MQ=255)
                                                                                                                                     atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                                                                                                                                                                                   <  6:986234/149‑1 (MQ=255)
                                                                                                                                     atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                                                                                                                                                                                   <  7:509738/149‑1 (MQ=255)
                                                                                                                                     atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                                                                                                                                                                                   >  5:1224487/1‑149 (MQ=255)
                                                                                                                                           ttGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTc                                                                                                                                                                                                                             >  8:418694/1‑149 (MQ=255)
                                                                                                                                            tGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCt                                                                                                                                                                                                                            <  7:138012/149‑1 (MQ=255)
                                                                                                                                             gCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAttttt                                                                                             >  2:274844/1‑275 (MQ=255)
                                                                                                                                              ccGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAc                                                                                                                                                                                                                          <  6:473387/149‑1 (MQ=255)
                                                                                                                                               cGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAcc                                                                                                                                                                                                                         >  5:909662/1‑149 (MQ=255)
                                                                                                                                                   ttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                                                                                                                                                                                                                                           >  7:1065303/1‑127 (MQ=255)
                                                                                                                                                      tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTgg                                                                                                                                                                                                                                                          >  7:1132455/1‑109 (MQ=255)
                                                                                                                                                      tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                                                                                                                                                  >  6:646143/1‑149 (MQ=255)
                                                                                                                                                      tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                                                                                                                                                  >  5:1409784/1‑149 (MQ=255)
                                                                                                                                                      tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                                                                                                                                                  >  6:843620/1‑149 (MQ=255)
                                                                                                                                                      tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                                                                                                                                                  >  8:763018/1‑149 (MQ=255)
                                                                                                                                                           cccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                                                                                                                                                                             <  6:728682/149‑1 (MQ=255)
                                                                                                                                                           cccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                                                                                                                                                                             <  7:1346089/149‑1 (MQ=255)
                                                                                                                                                            ccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGccc                                                                                                                                                                                                                                                   <  5:1159197/110‑1 (MQ=255)
                                                                                                                                                            ccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                                                                                                                                                                            <  8:1228234/149‑1 (MQ=255)
                                                                                                                                                             cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                                                                                                                           >  5:453064/1‑149 (MQ=255)
                                                                                                                                                             cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                                                                                                                           >  5:364388/1‑149 (MQ=255)
                                                                                                                                                             cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                                                                                                                           <  7:84769/149‑1 (MQ=255)
                                                                                                                                                             cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                                                                                                                           <  6:197349/149‑1 (MQ=255)
                                                                                                                                                             cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                                                                                                                           <  5:50879/149‑1 (MQ=255)
                                                                                                                                                             cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                                                                                                                           <  5:814645/149‑1 (MQ=255)
                                                                                                                                                               gACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                                                                                                                           <  7:680815/147‑1 (MQ=255)
                                                                                                                                                               gACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                                                                                                                                                                           >  8:680815/1‑147 (MQ=255)
                                                                                                                                                               gACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTAGCCAGCACTTCTACCAGCGTGTAGCTTTgg                                                                                                                                                                                                                                                          <  6:895803/100‑1 (MQ=255)
                                                                                                                                                                    tAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAtt                                                                                                                                                                                                    <  7:216160/149‑1 (MQ=255)
                                                                                                                                                                      gTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTg                                                                                                                                                                                                  <  8:1404599/149‑1 (MQ=255)
                                                                                                                                                                        ggTTTAACGCTGGTTCGTAATACCGATAGGTGCAATTTGAACGGAACTGATATCGGCTTTTGCCAGCACTTCTACCAGCGTGGGGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAg                                                                  <  4:243201/275‑1 (MQ=255)
                                                                                                                                                                         gTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCNGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGNCCTCCCGGTAGCTCGTGGTGCGTTTGCCTAGCGCGCCGAGATATTNTNTNTCAGTCATTTTTTATTGTccctttgg                                                                 >  3:41137/1‑270 (MQ=255)
                                                                                                                                                                          tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTggg                                                                                                                                                                                                                                     <  8:906760/110‑1 (MQ=255)
                                                                                                                                                                          tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                                                                                                                                                                             >  7:825753/1‑134 (MQ=255)
                                                                                                                                                                          tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                                                                                                                                                                             <  8:825753/134‑1 (MQ=255)
                                                                                                                                                                                gCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCt                                                                                                                                                                                                                   <  7:361636/122‑1 (MQ=255)
                                                                                                                                                                                  tGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgcg                                                                                                                                                                                      <  6:122819/149‑1 (MQ=255)
                                                                                                                                                                                           aTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTggg                                                                                                                                                                                                          >  8:360661/1‑120 (MQ=255)
                                                                                                                                                                                           aTACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGacac                                                                                                                                                                             >  5:578042/1‑149 (MQ=255)
                                                                                                                                                                                           aTACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGacac                                                                                                                                                                             <  5:1432837/149‑1 (MQ=255)
                                                                                                                                                                                            tACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACg                                                                                                                                                                            >  7:706448/1‑149 (MQ=255)
                                                                                                                                                                                                gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTaa                                                                                                                                                                        >  5:1350098/1‑149 (MQ=255)
                                                                                                                                                                                                gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTaa                                                                                                                                                                        <  6:1350098/149‑1 (MQ=255)
                                                                                                                                                                                                 aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCAt                                                                                                                                                                                                                                         >  2:51816/1‑83 (MQ=255)
                                                                                                                                                                                                 aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAAt                                                                                                                                                                       <  6:1371528/149‑1 (MQ=255)
                                                                                                                                                                                                 aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAAt                                                                                                                                                                       <  5:646143/149‑1 (MQ=255)
                                                                                                                                                                                                     cTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTggg                                                                                                                                                                                                                                     <  4:51816/83‑1 (MQ=255)
                                                                                                                                                                                                      tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat                                                                                                                                                                  <  5:947720/149‑1 (MQ=255)
                                                                                                                                                                                                      tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat                                                                                                                                                                  >  5:618699/1‑149 (MQ=255)
                                                                                                                                                                                                           tttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCAtt                                                                                                                                                                                                                                        <  5:1354352/74‑1 (MQ=255)
                                                                                                                                                                                                             tGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATc                                                                                                                                                           <  6:443061/149‑1 (MQ=255)
                                                                                                                                                                                                               aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGGTTTGGGTAGATTTTCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTTATCCATTACAACCGCGCGGGGGCTGTTGGTGCCCTGTTCGAGCGCGACGATATATTTGTTTTTAGTTCTGTTTTATTGTCCCGTAGTAATATTACATCGATGCTTTTTGTTTTTAATGATTtgt                            >  1:60364/1‑274 (MQ=255)
                                                                                                                                                                                                                acGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGt                                                                                                                                                         >  6:333625/2‑148 (MQ=255)
                                                                                                                                                                                                                acGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGt                                                                                                                                                         <  5:333625/147‑1 (MQ=255)
                                                                                                                                                                                                                  ggAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGAt                                                                                                                                                      >  8:1285012/1‑149 (MQ=255)
                                                                                                                                                                                                                   gAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATc                                                                                                                                                     >  5:1470187/1‑149 (MQ=255)
                                                                                                                                                                                                                   gAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATc                                                                                                                                                     <  8:185304/149‑1 (MQ=255)
                                                                                                                                                                                                                    aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATcc                                                                                                                                                    <  5:972390/149‑1 (MQ=255)
                                                                                                                                                                                                                     acggATATCGGCTTGCTCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCATATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAGTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTTCGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTcc                     <  2:200595/272‑1 (MQ=255)
                                                                                                                                                                                                                       tGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCAtt                                                                                                                                                 <  7:263438/149‑1 (MQ=255)
                                                                                                                                                                                                                        gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAACTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGCCACGCTAATGATATTGGCATCGTGATCCATTa                                                                                                                                                <  8:459844/149‑1 (MQ=255)
                                                                                                                                                                                                                        gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTa                                                                                                                                                <  8:1125260/149‑1 (MQ=255)
                                                                                                                                                                                                                           atCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATc                                                                                                                                                           <  6:404130/135‑1 (MQ=255)
                                                                                                                                                                                                                           atCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGa                                                                                                                                             <  5:113859/149‑1 (MQ=255)
                                                                                                                                                                                                                            tCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTg                                                                                                                                                        >  6:510847/1‑137 (MQ=255)
                                                                                                                                                                                                                            tCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATc                                                                                                                                                     >  6:562973/1‑140 (MQ=255)
                                                                                                                                                                                                                            tCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAAt                                                                                                                                                                                           <  5:877804/102‑1 (MQ=255)
                                                                                                                                                                                                                            tCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGAc                                                                                                                                            >  8:1447845/1‑149 (MQ=255)
                                                                                                                                                                                                                              ggCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcg                                                                                                                                          >  7:621166/1‑149 (MQ=255)
                                                                                                                                                                                                                              ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCt                                                                                                                     >  2:337379/1‑170 (MQ=255)
                                                                                                                                                                                                                                  ttCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTc                                                                                                                 <  4:337379/170‑1 (MQ=255)
                                                                                                                                                                                                                                   tCGCCAGCATGTCTACCAGCGGTTAGCTTTGGGGGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGAAACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCacaac       <  4:170034/275‑1 (MQ=255)
                                                                                                                                                                                                                                    cGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCggg                                                                                                                                    <  7:88104/149‑1 (MQ=255)
                                                                                                                                                                                                                                      ccAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                                                                                  >  5:542653/1‑149 (MQ=255)
                                                                                                                                                                                                                                        aGCACTTTGACCAGAGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGGCGTGTTCTACCCAACCTGGGTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGa  <  2:43569/275‑1 (MQ=255)
                                                                                                                                                                                                                                            cTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt                                                                                                                                                                 <  8:22337/112‑1 (MQ=255)
                                                                                                                                                                                                                                            cTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt                                                                                                                                                                 >  7:22337/1‑112 (MQ=255)
                                                                                                                                                                                                                                            cTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                                                                                                                            <  5:226353/149‑1 (MQ=255)
                                                                                                                                                                                                                                            cTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                                                                                                                            <  5:395075/149‑1 (MQ=255)
                                                                                                                                                                                                                                            cTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                                                                                                                            <  8:1165216/149‑1 (MQ=255)
                                                                                                                                                                                                                                             ttCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtgg                                                                                                                           >  7:354574/1‑149 (MQ=255)
                                                                                                                                                                                                                                             ttCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtgg                                                                                                                           <  6:820900/149‑1 (MQ=255)
                                                                                                                                                                                                                                             ttCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtgg                                                                                                                           <  6:1431354/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGAt                                                                                                                                                      >  6:104419/1‑117 (MQ=255)
                                                                                                                                                                                                                                                   cAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGCTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCt                                                                                                                     <  6:1199836/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        tGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCg                                                                                                                >  7:923983/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                
CACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGA  >  NC_000913/4116797‑4117303

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