Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,749,429 T→G 19.8% G66G (GGT→GGG)  sgbU → putative L‑xylulose 5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,749,4290TG19.8% 93.4 / 2.3 61G66G (GGT→GGG) sgbUputative L‑xylulose 5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (20/24);  new base G (13/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 6.70e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTC  >  NC_000913/3749206‑3749651
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                
aaacgacgaTTTCTGGGGAGTTGCGCGTGCGTAATCATCATGTAGGGAGTGATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTCTGATTTTGTCGAAATGACGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTgg                                                                                                                                                                                                                             <  1:583298/226‑1 (MQ=255)
 aTCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGCCTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGGGGTggg                                                                                                                                                                                                                            >  1:658128/1‑229 (MQ=255)
        tccgcTGGGGGGGTGGGTGTGCGTAATAAGCGGTGTGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGt                                                                                                                                                                                                                     <  2:315655/225‑1 (MQ=255)
         tttgtGGGGGGGTGGGCGTGGGTAATCCTCAGTTAGGGATTTATGAAAAAGCGGTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTTCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCGCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGt                                                                                                                                                                                                                     <  2:51853/224‑1 (MQ=255)
           tCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGAGCccgccg                                                                                                                                                                                                                   >  1:400179/1‑228 (MQ=255)
                ggAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGGGTTGGGATCCCGTCGATGTGc                                                                                                                                                                                                             >  1:369919/1‑229 (MQ=255)
                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGGAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGGGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGt                                                                                                                                                                                                         >  1:237827/1‑228 (MQ=255)
                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGTCCGGAGCGGCTGGGGCTGGCACAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAACCCGACGCACGCGCCTCACGTCTTGACTGGAGCGCCGGCCCAAG‑GACTTCCTCGGGTGCC‑CGAGTGATCGAAACAGGGGCTGGGAGCCCGTCCTTGGGCCttg                                                                                                                                                                                                         >  2:614260/1‑226 (MQ=255)
                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGGTGAAACCGACGAACGGGTCTCACGTCTTGACTGGAGGGCCGCACAAAG‑GACTTCCCTGGTTGCCGCGT‑GTTTCGAAACAGGGGTTTGGAACCCGAGGAGGGGCCTGtt                                                                                                                                                                                                        >  2:667663/1‑228 (MQ=255)
                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTGt                                                                                                                                                                                                         >  1:188302/1‑228 (MQ=255)
                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTGt                                                                                                                                                                                                         >  1:626675/1‑228 (MQ=255)
                       cgcgTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGGGCCTgcggt                                                                                                                                                                                                      >  1:315655/1‑225 (MQ=255)
                       cgcgTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGGGTTGGGATCCCGTCGAGGTGCCTGTCCg                                                                                                                                                                                                      >  1:105269/1‑229 (MQ=255)
                       cgcgTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTATCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTTTcc                                                                                                                                                                                                       >  1:419014/1‑228 (MQ=255)
                                        gTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGAGTTGGGATCCCGTCCATTTGCCTTTCCCCACATCCTCGCTTTccc                                                                                                                                                                                     >  1:364190/1‑229 (MQ=255)
                                                tttATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGGTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCCCGC‑TGGTCGAAACAGGGGTTGGGAGCCCGTCGATGGGCCTGGCCGCCAATTGTCGCTTTCCCTTTGGGAg                                                                                                                                                                             >  2:548458/1‑229 (MQ=255)
                                                 ttATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGGTTTGGGAGCCCGTCGAGTGGCCTGTCGGCACATCGTCGCTTTCCCTTTgtaat                                                                                                                                                                             >  1:413773/1‑224 (MQ=255)
                                                 ttATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCCGGTTGCCCCGA‑GTACTGAAACCGGTGTTGGGAGCCCGTCGATGTGCCCTGCCGCACATCCGCGCTTTCCCCTTTGGAGc                                                                                                                                                                            >  2:503520/1‑229 (MQ=255)
                                                   aTGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGCTCGAAACAGGGGTTGGGATCCCGTCGATGGGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCg                                                                                                                                                                          >  1:730792/1‑229 (MQ=255)
                                                           gcgcGGGCGGAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGGTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAg‑g                                                                                                                                                                 <  1:657124/229‑1 (MQ=255)
                                                               gggCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCGGGCAAAAAGCTGCGGTTTTGATTTTGGCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGt                                                                                                                                                             <  1:502672/228‑1 (MQ=255)
                                                                                  gCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCATGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATccc                                                                                                                                                                                                                       <  2:317872/153‑1 (MQ=255)
                                                                                  gCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGCAAACGGGGTTGGGATCCCGTCGATGTGCCTGTCCGCACCACGTCCCTTTCCCTTTGGTGGCCCGGGAGGG‑GCAGGGGGCGAACGGGCGCGGGaa                                                                                                                                          >  2:133128/1‑229 (MQ=255)
                                                                                  gCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATccc                                                                                                                                                                                                                       >  1:317872/1‑153 (MQ=255)
                                                                                  gCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGGG‑GCAGTGCGCGAACGGGCGCGGGaa                                                                                                                                          >  1:651654/1‑229 (MQ=255)
                                                                                       ggtggctggggctggCAAAAAGCCGCGGTGTTGATTTGGGCGGAGTGGCGGGGGATGAAAGCGACGAACGGCTCTCGCCTCTTGACTGGAGCGCCGCTCAAAG‑GACTTCGCTGGTTGCCGCGA‑TCATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCAt                                                                                                                                     <  2:140970/226‑1 (MQ=255)
                                                                                       aGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGGAATCAt                                                                                                                                     >  1:231657/1‑229 (MQ=255)
                                                                                       aGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCCTCGATGTGCCTGTCCGCCCATCGTCGCTTTCCCTTTGGTAGCCGTGACCAGCGCATGGCCCAAACGGCCG‑GGGAAACCAt                                                                                                                                     >  2:163112/1‑229 (MQ=255)
                                                                                            cTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGCGTTGGGATCCCGTCGATGTGCCTGTCCGACCATCGTCGCTTTCCCTTTGGGAGCCGTGACGAG‑GCAGTGCGCGAACGGGGGCGGGAAATCATggga                                                                                                                                 >  2:175271/1‑225 (MQ=255)
                                                                                              ggTTCTGGCAAAAAGCTGCCGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGGCGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTaaa                                                                                                                              <  1:667664/229‑1 (MQ=255)
                                                                                              ggTGTTGGCAAAAAGGTTGGGTTTTGATTTTTGCGAAATGTCGGTGAATGAAACCGACGAACGGCTCTCGCGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTaa                                                                                                                               <  2:237827/228‑1 (MQ=255)
                                                                                                       aaaagggttcGGTTTTGATTTTGTCGAAATGTCGGTGGATGAACACGACGAACGGCTCTCACGGCTTGACCGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGc                                                                                                                     <  2:105269/220‑1 (MQ=255)
                                                                                                          aaGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCCTTACCTTTGGTCGCCCTGCCGAG‑GCAGTGCGCGACCGGGCGCGGGACATCATGAGGAAAACGATTGGCGTg                                                                                                                  >  2:508316/1‑229 (MQ=255)
                                                                                                              tGCGGTTTTGTGTTTGTCGAAATTGGGGTGGAGGAAACCGACGAATGGCTCTCCCGTCTTGACTGGAGCGGCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGgcgc                                                                                                              <  2:730791/229‑1 (MQ=255)
                                                                                                               gggtgtttGTTTTGATCGAAATTGGGGTGGGTGAAACCGACTAACGGCCGTGACGTCTTGTCGGGCGCGCCGCTGAAAG‑GACTTCGCTGGTTGCCGCGA‑ATAGCGAAAGAGGTGTGTGGATCCCGTGGATGTGCCTGTCCGCACATCGTCGTTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCTCGGGAAATCATGAGTAAAGCGATTCGCCTGgcgcg                                                                                                             <  2:247584/224‑1 (MQ=255)
                                                                                                               gCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGCACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGgcg                                                                                                               >  1:598139/1‑227 (MQ=255)
                                                                                                               gCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCCCGAACGGGCGCGGGAAAACATGAGTAAAGCGATTCGCCTGgcggg                                                                                                             >  2:707338/1‑227 (MQ=255)
                                                                                                                        aTTGTGGCGAAATGTCGGTGGATGAAACCGACGAACGGTTCTCACGTCTTGACTGGATCGCCGCACAAAG‑GACGTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGGTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑Atctc                                                                                                     <  1:66796/227‑1 (MQ=255)
                                                                                                                                     gtggtgggctGAAACCGACGAAGGGCTCCCCCGTCTTGACTGGGGCGCCGCACAAGG‑GCCTTCGCGGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATc                                                                                      <  2:704707/220‑1 (MQ=255)
                                                                                                                                        ggTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGC‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTCGCCGTGACTAG‑GCATTCCGCAGACGGGCCCGGGAAACCCTCGGAAAACGGATTCGCCACGC‑CGCGCATCCCGGCATTAGcccacaccga                                                                                   >  2:709971/1‑220 (MQ=255)
                                                                                                                                            gATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGGGACGAG‑GCAGTCGGCGAACGTGCGCGGGAAATAATGAGCAAAGCGAGTCGCCTGGCGCGCG‑ATCTCGGGCTTCGGACCATCCCACTgg                                                                               >  2:659350/1‑229 (MQ=255)
                                                                                                                                             aTGAAACGGACGGACGGCTTTGACGTCTTGACTGGGTCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACAGCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGc                                                                              <  2:626674/229‑1 (MQ=255)
                                                                                                                                                 aaCCGACGAACGGCTCTCACGTCTTTACTGGCGCTCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAgg                                                                           <  1:614261/228‑1 (MQ=255)
                                                                                                                                                     gacgaACGGCTCTGACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTAt                                                                       <  1:139314/228‑1 (MQ=255)
                                                                                                                                                     gacgaACGGCTCTCCCGTCTGGACTGGAGCGCCGCACAAAG‑GCCTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTAt                                                                       <  1:556133/228‑1 (MQ=255)
                                                                                                                                                      acgaACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGCAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑GTCTCGGCATTCGCCCCGTCCCGCTTGACGGTTTTGa                                                                     >  1:285103/1‑229 (MQ=255)
                                                                                                                                                      acgaACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGCTATGa                                                                     >  1:602772/1‑229 (MQ=255)
                                                                                                                                                       ctagcGGCTCTCCCGTTTTGACTGGGTCGCTGCCCAAGG‑GACTTCTCTGGTTTGCGGCA‑TGAACGGAACGGGTGTTGGGAGCCCGTCGATGTGGCTGGCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGa                                                                     <  2:457961/224‑1 (MQ=255)
                                                                                                                                                       cgaACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTTTGAc                                                                    >  1:659715/1‑229 (MQ=255)
                                                                                                                                                        ggagggCTCTCCCGTCTTGGAGGGAGGGGCG‑GCAAAGTGACTTCGCTGGGTGCCGCTT‑GGATCGGAACAGGTGTTGGGATCCCGTCGAGGTGCCTGTCCGCGCATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGGGGGAAATCATGAGTAAAGCGATTCGACTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGAc                                                                    <  2:364190/224‑1 (MQ=255)
                                                                                                                                                        gaACGGCCTTCACGGTTTGACTGGAGCGCCGGACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGGCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGAc                                                                    <  2:651653/228‑1 (MQ=255)
                                                                                                                                                                              ggAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACCGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCAcgacga                                             <  1:133128/229‑1 (MQ=255)
                                                                                                                                                                                       cacaAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCCCGCCGCAGGCACCCg                                    >  1:207091/1‑229 (MQ=255)
                                                                                                                                                                                          aaaG‑GCCTTGGCGTGTTGCCGCGA‑TGATCGAAACAGGTGTGGGGATCCCGTCGATGTGCCTGGCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGc                                  <  2:658127/228‑1 (MQ=255)
                                                                                                                                                                                                  ttCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGCGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCGCCATCCAGCTGGCAGGTTATGACGTCTGTTCCGAAGACCACGACGGAGGCACCCGGCAACGtttt                          >  1:302190/1‑229 (MQ=255)
                                                                                                                                                                                                     gCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCt                       >  2:529964/1‑229 (MQ=255)
                                                                                                                                                                                                     gCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTGTTGCt                       >  1:712488/1‑229 (MQ=255)
                                                                                                                                                                                                      ctgggcgCCGGGA‑TGATCGAAACGGGGGTTGGGATCCCGTCGATGTGCCTGTGCGCACATCGTCGGTGTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTg                      <  1:691169/223‑1 (MQ=255)
                                                                                                                                                                                                        ggTTGCCGCGA‑TGATCGGAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTg                      <  1:509178/227‑1 (MQ=255)
                                                                                                                                                                                                                      atCGAAACAGGGTTTGGGTTCCGGTCGGGGGGCCCGTCCCGCCATGGGCGCTTTCCCTTTGGTAGCCGGGCCGGG‑GGCGGGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTggg       <  2:188302/229‑1 (MQ=255)
                                                                                                                                                                                                                           aaCAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGACTATTACGAAGACCACGACGAAGGCACCCTGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGtt  >  1:84888/1‑228 (MQ=255)
                                                                                                                                                                                                                                gtgtTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCgg                                   >  1:668095/1‑191 (MQ=255)
                                                                                                                                                                                                                                gtgtTGGGATCCCGTCGATGTGCCTGTCCGCACAGCGTCGCTTCCCCATTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCgg                                   <  2:668094/191‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                
AATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAG‑GACTTCGCTGGTTGCCGCGA‑TGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAG‑GCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCG‑ATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTC  >  NC_000913/3749206‑3749651

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