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→GGGsgbUputative 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

AATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTC  >  NC_000913/3749206‑3749651
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                
aaacgacgaTTTCTGGGGAGTTGCGCGTGCGTAATCATCATGTAGGGAGTGATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTCTGATTTTGTCGAAATGACGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTgg                                                                                                                                                                                                                             <  1:583298/226‑1 (MQ=255)
 aTCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGCCTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGGTggg                                                                                                                                                                                                                            >  1:658128/1‑229 (MQ=255)
        tccgcTGGGGGGGTGGGTGTGCGTAATAAGCGGTGTGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGt                                                                                                                                                                                                                     <  2:315655/225‑1 (MQ=255)
         tttgtGGGGGGGTGGGCGTGGGTAATCCTCAGTTAGGGATTTATGAAAAAGCGGTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTTCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCGCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGt                                                                                                                                                                                                                     <  2:51853/224‑1 (MQ=255)
           tCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGAGCccgccg                                                                                                                                                                                                                   >  1:400179/1‑228 (MQ=255)
                ggAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTCGATGTGc                                                                                                                                                                                                             >  1:369919/1‑229 (MQ=255)
                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGGAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGGGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGt                                                                                                                                                                                                         >  1:237827/1‑228 (MQ=255)
                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGTCCGGAGCGGCTGGGGCTGGCACAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAACCCGACGCACGCGCCTCACGTCTTGACTGGAGCGCCGGCCCAAGGACTTCCTCGGGTGCCCGAGTGATCGAAACAGGGGCTGGGAGCCCGTCCTTGGGCCttg                                                                                                                                                                                                         >  2:614260/1‑226 (MQ=255)
                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGGTGAAACCGACGAACGGGTCTCACGTCTTGACTGGAGGGCCGCACAAAGGACTTCCCTGGTTGCCGCGTGTTTCGAAACAGGGGTTTGGAACCCGAGGAGGGGCCTGtt                                                                                                                                                                                                        >  2:667663/1‑228 (MQ=255)
                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTGt                                                                                                                                                                                                         >  1:188302/1‑228 (MQ=255)
                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTGt                                                                                                                                                                                                         >  1:626675/1‑228 (MQ=255)
                       cgcgTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGGGCCTgcggt                                                                                                                                                                                                      >  1:315655/1‑225 (MQ=255)
                       cgcgTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTCGAGGTGCCTGTCCg                                                                                                                                                                                                      >  1:105269/1‑229 (MQ=255)
                       cgcgTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTATCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTTTcc                                                                                                                                                                                                       >  1:419014/1‑228 (MQ=255)
                                        gTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGAGTTGGGATCCCGTCCATTTGCCTTTCCCCACATCCTCGCTTTccc                                                                                                                                                                                     >  1:364190/1‑229 (MQ=255)
                                                tttATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGGTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCCCGCTGGTCGAAACAGGGGTTGGGAGCCCGTCGATGGGCCTGGCCGCCAATTGTCGCTTTCCCTTTGGGAg                                                                                                                                                                             >  2:548458/1‑229 (MQ=255)
                                                 ttATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGTTTGGGAGCCCGTCGAGTGGCCTGTCGGCACATCGTCGCTTTCCCTTTgtaat                                                                                                                                                                             >  1:413773/1‑224 (MQ=255)
                                                 ttATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCCGGTTGCCCCGAGTACTGAAACCGGTGTTGGGAGCCCGTCGATGTGCCCTGCCGCACATCCGCGCTTTCCCCTTTGGAGc                                                                                                                                                                            >  2:503520/1‑229 (MQ=255)
                                                   aTGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGCTCGAAACAGGGGTTGGGATCCCGTCGATGGGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCg                                                                                                                                                                          >  1:730792/1‑229 (MQ=255)
                                                           gcgcGGGCGGAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGGTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAgg                                                                                                                                                                 <  1:657124/229‑1 (MQ=255)
                                                               gggCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCGGGCAAAAAGCTGCGGTTTTGATTTTGGCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGt                                                                                                                                                             <  1:502672/228‑1 (MQ=255)
                                                                                  gCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCATGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATccc                                                                                                                                                                                                                       <  2:317872/153‑1 (MQ=255)
                                                                                  gCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGCAAACGGGGTTGGGATCCCGTCGATGTGCCTGTCCGCACCACGTCCCTTTCCCTTTGGTGGCCCGGGAGGGGCAGGGGGCGAACGGGCGCGGGaa                                                                                                                                          >  2:133128/1‑229 (MQ=255)
                                                                                  gCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATccc                                                                                                                                                                                                                       >  1:317872/1‑153 (MQ=255)
                                                                                  gCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGGGGCAGTGCGCGAACGGGCGCGGGaa                                                                                                                                          >  1:651654/1‑229 (MQ=255)
                                                                                       ggtggctggggctggCAAAAAGCCGCGGTGTTGATTTGGGCGGAGTGGCGGGGGATGAAAGCGACGAACGGCTCTCGCCTCTTGACTGGAGCGCCGCTCAAAGGACTTCGCTGGTTGCCGCGATCATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCAt                                                                                                                                     <  2:140970/226‑1 (MQ=255)
                                                                                       aGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGGAATCAt                                                                                                                                     >  1:231657/1‑229 (MQ=255)
                                                                                       aGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCCTCGATGTGCCTGTCCGCCCATCGTCGCTTTCCCTTTGGTAGCCGTGACCAGCGCATGGCCCAAACGGCCGGGGAAACCAt                                                                                                                                     >  2:163112/1‑229 (MQ=255)
                                                                                            cTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGCGTTGGGATCCCGTCGATGTGCCTGTCCGACCATCGTCGCTTTCCCTTTGGGAGCCGTGACGAGGCAGTGCGCGAACGGGGGCGGGAAATCATggga                                                                                                                                 >  2:175271/1‑225 (MQ=255)
                                                                                              ggTTCTGGCAAAAAGCTGCCGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGGCGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTaaa                                                                                                                              <  1:667664/229‑1 (MQ=255)
                                                                                              ggTGTTGGCAAAAAGGTTGGGTTTTGATTTTTGCGAAATGTCGGTGAATGAAACCGACGAACGGCTCTCGCGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTaa                                                                                                                               <  2:237827/228‑1 (MQ=255)
                                                                                                       aaaagggttcGGTTTTGATTTTGTCGAAATGTCGGTGGATGAACACGACGAACGGCTCTCACGGCTTGACCGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGc                                                                                                                     <  2:105269/220‑1 (MQ=255)
                                                                                                          aaGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCCTTACCTTTGGTCGCCCTGCCGAGGCAGTGCGCGACCGGGCGCGGGACATCATGAGGAAAACGATTGGCGTg                                                                                                                  >  2:508316/1‑229 (MQ=255)
                                                                                                              tGCGGTTTTGTGTTTGTCGAAATTGGGGTGGAGGAAACCGACGAATGGCTCTCCCGTCTTGACTGGAGCGGCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGgcgc                                                                                                              <  2:730791/229‑1 (MQ=255)
                                                                                                               gggtgtttGTTTTGATCGAAATTGGGGTGGGTGAAACCGACTAACGGCCGTGACGTCTTGTCGGGCGCGCCGCTGAAAGGACTTCGCTGGTTGCCGCGAATAGCGAAAGAGGTGTGTGGATCCCGTGGATGTGCCTGTCCGCACATCGTCGTTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCTCGGGAAATCATGAGTAAAGCGATTCGCCTGgcgcg                                                                                                             <  2:247584/224‑1 (MQ=255)
                                                                                                               gCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGCACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGgcg                                                                                                               >  1:598139/1‑227 (MQ=255)
                                                                                                               gCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCCCGAACGGGCGCGGGAAAACATGAGTAAAGCGATTCGCCTGgcggg                                                                                                             >  2:707338/1‑227 (MQ=255)
                                                                                                                        aTTGTGGCGAAATGTCGGTGGATGAAACCGACGAACGGTTCTCACGTCTTGACTGGATCGCCGCACAAAGGACGTCGCTGGTTGCCGCGATGATCGAAACAGGTGGTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGAtctc                                                                                                     <  1:66796/227‑1 (MQ=255)
                                                                                                                                     gtggtgggctGAAACCGACGAAGGGCTCCCCCGTCTTGACTGGGGCGCCGCACAAGGGCCTTCGCGGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATc                                                                                      <  2:704707/220‑1 (MQ=255)
                                                                                                                                        ggTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGCTGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTCGCCGTGACTAGGCATTCCGCAGACGGGCCCGGGAAACCCTCGGAAAACGGATTCGCCACGCCGCGCATCCCGGCATTAGcccacaccga                                                                                   >  2:709971/1‑220 (MQ=255)
                                                                                                                                            gATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGGGACGAGGCAGTCGGCGAACGTGCGCGGGAAATAATGAGCAAAGCGAGTCGCCTGGCGCGCGATCTCGGGCTTCGGACCATCCCACTgg                                                                               >  2:659350/1‑229 (MQ=255)
                                                                                                                                             aTGAAACGGACGGACGGCTTTGACGTCTTGACTGGGTCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACAGCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGc                                                                              <  2:626674/229‑1 (MQ=255)
                                                                                                                                                 aaCCGACGAACGGCTCTCACGTCTTTACTGGCGCTCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAgg                                                                           <  1:614261/228‑1 (MQ=255)
                                                                                                                                                     gacgaACGGCTCTGACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTAt                                                                       <  1:139314/228‑1 (MQ=255)
                                                                                                                                                     gacgaACGGCTCTCCCGTCTGGACTGGAGCGCCGCACAAAGGCCTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTAt                                                                       <  1:556133/228‑1 (MQ=255)
                                                                                                                                                      acgaACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGCAATCATGAGTAAAGCGATTCGCCTGGCGCGCGGTCTCGGCATTCGCCCCGTCCCGCTTGACGGTTTTGa                                                                     >  1:285103/1‑229 (MQ=255)
                                                                                                                                                      acgaACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGCTATGa                                                                     >  1:602772/1‑229 (MQ=255)
                                                                                                                                                       ctagcGGCTCTCCCGTTTTGACTGGGTCGCTGCCCAAGGGACTTCTCTGGTTTGCGGCATGAACGGAACGGGTGTTGGGAGCCCGTCGATGTGGCTGGCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGa                                                                     <  2:457961/224‑1 (MQ=255)
                                                                                                                                                       cgaACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTTTGAc                                                                    >  1:659715/1‑229 (MQ=255)
                                                                                                                                                        ggagggCTCTCCCGTCTTGGAGGGAGGGGCGGCAAAGTGACTTCGCTGGGTGCCGCTTGGATCGGAACAGGTGTTGGGATCCCGTCGAGGTGCCTGTCCGCGCATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGGGGGAAATCATGAGTAAAGCGATTCGACTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGAc                                                                    <  2:364190/224‑1 (MQ=255)
                                                                                                                                                        gaACGGCCTTCACGGTTTGACTGGAGCGCCGGACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGGCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGAc                                                                    <  2:651653/228‑1 (MQ=255)
                                                                                                                                                                              ggAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACCGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCAcgacga                                             <  1:133128/229‑1 (MQ=255)
                                                                                                                                                                                       cacaAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCCCGCCGCAGGCACCCg                                    >  1:207091/1‑229 (MQ=255)
                                                                                                                                                                                          aaaGGCCTTGGCGTGTTGCCGCGATGATCGAAACAGGTGTGGGGATCCCGTCGATGTGCCTGGCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGc                                  <  2:658127/228‑1 (MQ=255)
                                                                                                                                                                                                  ttCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGCGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCGCCATCCAGCTGGCAGGTTATGACGTCTGTTCCGAAGACCACGACGGAGGCACCCGGCAACGtttt                          >  1:302190/1‑229 (MQ=255)
                                                                                                                                                                                                     gCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCt                       >  2:529964/1‑229 (MQ=255)
                                                                                                                                                                                                     gCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTGTTGCt                       >  1:712488/1‑229 (MQ=255)
                                                                                                                                                                                                      ctgggcgCCGGGATGATCGAAACGGGGGTTGGGATCCCGTCGATGTGCCTGTGCGCACATCGTCGGTGTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTg                      <  1:691169/223‑1 (MQ=255)
                                                                                                                                                                                                        ggTTGCCGCGATGATCGGAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTg                      <  1:509178/227‑1 (MQ=255)
                                                                                                                                                                                                                      atCGAAACAGGGTTTGGGTTCCGGTCGGGGGGCCCGTCCCGCCATGGGCGCTTTCCCTTTGGTAGCCGGGCCGGGGGCGGGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTggg       <  2:188302/229‑1 (MQ=255)
                                                                                                                                                                                                                           aaCAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGACTATTACGAAGACCACGACGAAGGCACCCTGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGtt  >  1:84888/1‑228 (MQ=255)
                                                                                                                                                                                                                                gtgtTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCgg                                   >  1:668095/1‑191 (MQ=255)
                                                                                                                                                                                                                                gtgtTGGGATCCCGTCGATGTGCCTGTCCGCACAGCGTCGCTTCCCCATTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCgg                                   <  2:668094/191‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                
AATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTC  >  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: