Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,464,859 A→C 15.6% V441G (GTG→GGG)  treB ← trehalose‑specific PTS enzyme: IIB and IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,464,8590AC15.6% 149.8 / 4.3 79V441G (GTG→GGG) treBtrehalose‑specific PTS enzyme: IIB and IIC component
Reads supporting (aligned to +/- strand):  ref base A (41/18);  new base C (2/17);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 7.59e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGC  >  NC_000913/4464652‑4465075
                                                                                                                                                                                                                |                                                                                                                                                                                                                        
tctTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCCCGGTTGaa                                                                                                                                                                                                     >  1:78375/1‑229 (MQ=255)
tctTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAACCCAAAACCCTGCCAGTAGCTCGGttta                                                                                                                                                                                                      >  2:549551/1‑226 (MQ=255)
tctTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACAAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGTCAGGTTGaa                                                                                                                                                                                                     >  1:483484/1‑229 (MQ=255)
  tttGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAAt                                                                                                                                                                                                    >  1:216530/1‑228 (MQ=255)
  tttGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGCGAGTACAATCGGGATGATGATGGCGATAGCCCTTGCCAGCGCAAACCCCTGCCAGTACCCTGGTTGAAt                                                                                                                                                                                                    >  1:137394/1‑228 (MQ=255)
  tttGGATAAATCTGGGAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAAt                                                                                                                                                                                                    <  1:497548/228‑1 (MQ=255)
        ttaaTTTGGTAAATAACCGTGTTTTGCCACCAGGGGGGAAATAGGTTCATTACGTTAATCCTGCGAATGCGAGGGGGCGCAATTGCGACCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCgagaga                                                                                                                                                                                             <  2:300438/227‑1 (MQ=255)
                gTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTCCCAGTAGCCCGGTTGAATCGAGAGAATTCCCgg                                                                                                                                                                                     >  1:74197/1‑229 (MQ=255)
                gTAGATAACGCCGTTTTGCCACCAGAGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGCGTACAATCGGGATGACCATGGCGATAGCCATTGCCAGCGCAAAAACCTTCCAGACGCTCGGTGGCATCGAGAGAGTTCCCgg                                                                                                                                                                                     >  2:389055/1‑229 (MQ=255)
                      aaCGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGCGAATTCCAggacggcg                                                                                                                                                                               >  1:158133/1‑223 (MQ=255)
                       aCGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGgcgg                                                                                                                                                                              >  1:534381/1‑227 (MQ=255)
                          ccGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGCAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGGTAAAAGAGGTGAGGACAATCGGGATGATGAGTGCGATAGCCATTGCCAGCGCAACAACCTCCCAGTAGCTCGGTTGAACCGCGGAAGATGCCGGCAGgccgcc                                                                                                                                                                            >  2:274527/1‑228 (MQ=255)
                                 gCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGa                                                                                                                                                                    <  1:27963/229‑1 (MQ=255)
                                     ccaGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATTGAGAGAATTCCCGGCCGGCCGCCTCAGCCGAGGcc                                                                                                                                                                >  1:361870/1‑229 (MQ=255)
                                      caGGGGGGAAGGAGGGTCATTAGGTTATTCAGGGGAATGGCGGGGGGGGCAATTGGCCCCCGAAGAAAATTAACAAATGTCCAGTGGTCCCAGGGGGTATTTCCGCTGAAAGATAAACGAGGTGAGTACAATCGGGATGATTATGGCGATAGCCATTGCCAGCGCAAACCCCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCa                                                                                                                                                               <  2:261300/229‑1 (MQ=255)
                                         tGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGACGCTCGGGTGGATCGAGGAAGTTACCGGAAGGCCGCCTACGCCGCTGCCATTc                                                                                                                                                            >  2:512407/1‑229 (MQ=255)
                                         gggggAAGATGGATCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCTGGCGGGATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCAtt                                                                                                                                                             <  2:216530/227‑1 (MQ=255)
                                                gATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCCCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGcccaac                                                                                                                                                     >  1:469822/1‑225 (MQ=255)
                                                gATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCGTCCAGCAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTCCCCCCATGCCATTCCCCaaaa                                                                                                                                                     >  2:724137/1‑229 (MQ=255)
                                                                      gCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGCTGCTGATGGCGATAGCCATTTCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGTGAGAATTCCCGGCAGGCCGCCTCCGCCCATGCAATTCGCCATCACCCCGTTCAGGCCGCaaacca                                                                                                                                >  2:720811/1‑223 (MQ=255)
                                                                      gCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATCGCa                                                                                                                                >  1:19665/1‑228 (MQ=255)
                                                                      gCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCCCATAGCa                                                                                                                                >  1:349267/1‑228 (MQ=255)
                                                                       ggAATGCGGGGGGCGGAAATTGTCCCCCGAAAAAAATTAAACAAGGGCCAGCGGGCCCAGGCGTTATTTCCGTTGAATGATAAACGAGGTGAGTAAAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAAt                                                                                                                              <  2:282922/228‑1 (MQ=255)
                                                                       cGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGGATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAAAGCCGTTCAcgccgccagccaa                                                                                                                               >  2:643889/1‑221 (MQ=255)
                                                                                 ggggccaatttGCGCCCGTAAAAAAATTAAACAATGGCCAGCGTCCCCGGCGAGTATTTCCGCTGATAGATAAACGGGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCCAAAACCCCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGccagac                                                                                                                    <  2:60810/220‑1 (MQ=255)
                                                                                 ggggCGAATTTGGCCCCGGAAGAAAATTAAACAAGGTCCAGGGTGCCCGGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTCCAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGccagac                                                                                                                    <  2:534381/229‑1 (MQ=255)
                                                                                 ggggCCAAATTGCGCCCCGAAGAAAATTAAACAATGCCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGccagac                                                                                                                    <  1:17500/229‑1 (MQ=255)
                                                                                   ggCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGccagcac                                                                                                                   >  1:290924/1‑225 (MQ=255)
                                                                                       cAATTGGGCGCCATAAAAAATTAAAACATGTCGAGCGTGCCCAGGAGAGATTACCGCTGATAGACAAACCAGTTGAGTCCAATCGGGATGATGATGGCGGTAGCAGTTGCACGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGTCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCCTAGCAAACCTGccagagcaga                                                                                                                <  2:728211/227‑1 (MQ=255)
                                                                                        aaTTGCGCCCCAAAGAAAATAAAACAAAGTCCAAGATGCCCGGGCGGTATTTCCGGTGAAAGAAAAACGAGGTGAGTACAATCGGGATGAGGATGGCGATAGCCATGGCCAGCGCAAACCCCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAAc                                                                                                              <  2:557312/228‑1 (MQ=255)
                                                                                               ccccGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGGTGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTAAGGCCGCATAGCAATCCTGCCAGCCACGAACCcatcat                                                                                                       >  1:310705/1‑228 (MQ=255)
                                                                                                  agaagaaAATTAAACAATGGCCAGGGTGCCCAGGCGGTATTTCCGCTGAAAGATAAACCAGTTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATcgcg                                                                                                   <  1:512407/228‑1 (MQ=255)
                                                                                                    agaaaatTTAAAAAATGTCCAGCGTGACCGGGGGGTATTTCCGCGGAGAGATAAACGAGGGGAGTCCAATCGGGTTGAGGAGGGCTATAGCCATGGCCGGCGCAACCACCTGCCAGTAGCTCGGTTGAAGCGAGAGAATGCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGca                                                                                                 <  2:158133/226‑1 (MQ=255)
                                                                                                    aagaaAATTAACAACTGTCCGGCGTGCCCAGGGGTTATTTCCGCTAATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCCGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGca                                                                                                 <  2:25805/229‑1 (MQ=255)
                                                                                                    aagaaAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCGGACCCGAACCacacaccccgcc                                                                                                 >  1:97327/1‑218 (MQ=255)
                                                                                                    aagaaAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCCGACAAGAACCAATCATCCCGca                                                                                                 >  1:329354/1‑229 (MQ=255)
                                                                                                    aagaaAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGAAAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCCTTGCCAGCGCAAACACCTGCCAGTACCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACCCCGCTGCCATTCGCCATAAATCCAGTCCGGCCGCATAGCAACCCTGCCAGACCCGAACCAAACATCGCCca                                                                                                 >  2:670645/1‑229 (MQ=255)
                                                                                                    aaaaaaaTTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCGGCTGATTGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACTTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGca                                                                                                 <  1:447202/226‑1 (MQ=255)
                                                                                                    aaaaaaaTTAAAAAAGGGCCGGTGTGCCGGGGCGTGATTGCCGCTGATAGAGAAACGGGGTGGGTACATCCGGGATGAGGATGGCAATAGCCATTGCCAGCGCAAACTCCTGCCAGTAGCTCGGTTGAATCGAGAGAAGTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGca                                                                                                 <  2:121510/226‑1 (MQ=255)
                                                                                                      gataaTTAAAGAGTGTACAGCGTGGCCGGGGGGTATTCCCGCTAGTAGATAAACGGGGTGAGTACAATCGGGATGTTGATGGCGATGGCCTGGGCCGGCGAAACCCCCTGCCAGTAGCTCGGTTGCATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGAAGCCTATCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGcaca                                                                                               <  2:137394/226‑1 (MQ=255)
                                                                                                         atttAAACAAGTTCCAGTGTGCCCGGGCGGTTTTTCCGCTGATAGAAAAACGAGGTGAGTACAATCGGGATGATGATGGGGATGGCCATTGCCAGCGCAAACCCCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCa                                                                                            <  2:662707/227‑1 (MQ=255)
                                                                                                                   tgcccAGGGGGCCCGGGCGGTATTTCCGCTGATAGATAAACGGGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCg                                                                                  <  1:600950/226‑1 (MQ=255)
                                                                                                                   ggTCCAGGGTGCCCGGGGGGTTTTTCCGCTGAAAGATAAACGGGGTGAGTACAATGGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACCCCTGCCAGTAGCCCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCg                                                                                  <  2:183914/228‑1 (MQ=255)
                                                                                                                     tCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCGGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGAGTCCATTACCCATAAGTCCGTTCAGGCCGCCTAGCAATCCTGCCCACCCCGAACCAACCCTCCCCGACAGCCTCGGGaaaggat                                                                                >  2:595692/1‑224 (MQ=255)
                                                                                                                          cGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTAAGGCCGCCTAGCAATCCCGCCAGACCAGAAACAATCATCGCGCACAGCCACGGGAACGGAGAttt                                                                            >  2:450961/1‑228 (MQ=255)
                                                                                                                              cccAGGGGGTATTTCCGCTGAAAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGGAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCg                                                                                  <  1:311759/218‑1 (MQ=255)
                                                                                                                              cccAGGCGGTATTTCCGGTGATAGAAAAACGAGGTTAGTACAAGCGGGATTATGATGGCGATAGCCAGTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGt                                                                       <  1:643890/229‑1 (MQ=255)
                                                                                                                              cccAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGACTCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACGCCATCGGGAAGCGATATTTCAgg                                                                        >  1:411410/1‑228 (MQ=255)
                                                                                                                              cccAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGAGAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCCCCATAACGCCGTTCAGGCCGCATAGCCATCCCGCCAGACCAGACCCACACACCGCGCACAGCACCCGGAAACGACATTTccagg                                                                       >  2:311759/1‑225 (MQ=255)
                                                                                                                              ccaggggggTTTTTCCGCGAATAGATAAAAGGGGTGAGGAAAATCGGGATGGAGATGGCGATAGCCATTGCCAGGGAAAACCCCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGt                                                                       <  2:427754/222‑1 (MQ=255)
                                                                                                                               ccAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAACGGAGATTTCCGGtt                                                                      >  1:223758/1‑229 (MQ=255)
                                                                                                                                    gggatttttCCGTTTGTAGATAAACGGGGTGGATCCAATCGGTATGATGTTGGCGTAAGCCTGTGCCAGGGCAAGCCCCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCAGTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATg                                                                  <  2:535325/223‑1 (MQ=255)
                                                                                                                                    gctgTATGTCCGCAGAAAGATAAACGAGGTGAGTACAACCGGGATTATTATGGCGATAGCCAGTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCCTCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATg                                                                  <  1:687205/226‑1 (MQ=255)
                                                                                                                                          tttCCGGTGTTAGATAAACGGGGTGATGCCAATCGGGTGGATGGTGGCGATAGCCATTCACCGCGCAAACCCCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTAc                                                            <  2:349931/229‑1 (MQ=255)
                                                                                                                                          tttCCACTGATAGGTAAACGAGGTTAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACCCCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTa                                                             <  1:546874/228‑1 (MQ=255)
                                                                                                                                                  gatagataAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCCTAGCAATCCTGCCAGACCAGAACCAACCATCGCGCACAGCCTCGGGAAACGATACTTCAGGATGATGCCGTACATTGAAg                                                     >  2:464105/1‑228 (MQ=255)
                                                                                                                                                  gatagataAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAg                                                     >  1:391799/1‑228 (MQ=255)
                                                                                                                                                  gatagataAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCCCGACCAGCATCGGGAAGCGATATTTCAGGTTGATGCCGAACATTGCAg                                                     >  2:200407/1‑228 (MQ=255)
                                                                                                                                                   aaagagaAACGAGGGGAATACAACCGGGATGATGATGGCGATAGCAATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGGCCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAgg                                                    <  1:74060/226‑1 (MQ=255)
                                                                                                                                                                   aGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGTAGCGATATTTCAGGTTGCGACCGTACATTGCAGGCTCAGTGACACCAAGCCa                                  >  1:284394/1‑229 (MQ=255)
                                                                                                                                                                     tACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCCTAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGACTCCGGAAGCGATATTTCAGGGTGCTGCCGCACATTGCCGGCCCAGTGCACCCAAGcccg                                 >  2:645430/1‑227 (MQ=255)
                                                                                                                                                                     tACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGGAGCTAGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCAATTCGCCGTACAGCCGTTCAGCCCGCCAAGCCAACCGCCCAGACCAGAACCAATCATCTGTCACAGCACCGGGGAGCGATATTTCAGTTTGTTGCCGAACCTTGCGGGCCCAAGGCCCCCAGACCAg                                 >  2:598999/1‑229 (MQ=255)
                                                                                                                                                                      aaaagcGGGGTGAGGAGGGCGAAGACCAGTGCCCGCGCAACCCCCGGCCAGTAGCTCGGTTGAATCGAGGGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAgg                                <  2:78375/224‑1 (MQ=255)
                                                                                                                                                                          tCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCCTAGCAATCCTGCCAGCCCAAAACCAATCATCGAGCACAGCACCGGGCAGCGAAATTTCAGGTTTATGCCGTACATCACAGGCTACGGGCCCCCAAAcccgcggga                            >  2:589943/1‑222 (MQ=255)
                                                                                                                                                                          tCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCCCAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCACAGCCAGGCGga                            >  1:225735/1‑229 (MQ=255)
                                                                                                                                                                          tCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTCGCAGACCCGAACCACTCATCGGCCACACACTCGGGCAACACTATTTACGGTGGTACCCGGACATTACAGGCCCGGGGACCCCCAACCAGggggg                            >  2:735570/1‑228 (MQ=255)
                                                                                                                                                                          tCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACACCGATCCCATTCGCCATAACGCCGTTCAGGCCGCGTAGCAATCCTGCCAGACCAGAACCAACCATCGCGCCCACACCCGGGTAGCGAACTTTCAGGTTGATGCCGGACATTGAAGACCCCGTGACCCCAGGcacgggggg                            >  2:225765/1‑221 (MQ=255)
                                                                                                                                                                             ggATGATGATGGCGATAGCCATTGCCAGCGCAGACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGAGgaga                          >  1:31752/1‑228 (MQ=255)
                                                                                                                                                                                   tgtgggCGATAGCCATTGCCCGCCAACACCCCTGCCGGGAGCTCGGTTGAAGCGAGAGAATCCCCGGCAGGCCTCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCggc                   <  2:310705/225‑1 (MQ=255)
                                                                                                                                                                                   tgtggcggaaggCCTTTGCCATCGCAAACACCTGCCAGTTGTCCGGTTGAACTGAGAGAATTCCCGGCCGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCggc                   <  2:168618/218‑1 (MQ=255)
                                                                                                                                                                                       ggCGGTGGCCATGTCCAGGGGAAACCCCTGGCAGTATCTCGGTTGAATCGAAAGAATTCCCGGAAGGCCGCATCCGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCggcagg                <  2:349267/228‑1 (MQ=255)
                                                                                                                                                                                          gATGGCCATTGCCAGCGCAAACCCCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACg            <  1:531365/229‑1 (MQ=255)
                                                                                                                                                                                             aGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGCCCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCCGGTTGATGCCGTACATTGAAGGACCATGGACCCCAAGCCAGGCGGGGATAGCGGAAAGCAcagg          >  1:208793/1‑225 (MQ=255)
                                                                                                                                                                                                cgtggCCAGGCCAAACCCCGGCCAGTAAGCCGGGTGAATCGAGAGAATTCCGGGCAGGCCGGCCACGGCGATGCCATTCGACATAACGCCGTTCAGGTCGCATAGCAATCCGGCCAGAACAGAACCAATCATCGCGGACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGTCCCCAAGCCAGGCGGAGATAGCGGCAGGCACggctat       <  2:603916/224‑5 (MQ=255)
                                                                                                                                                                                                 atctcccgccaaaaCCCCTTCCAGTCGCTGGGTTGAATCGAGACAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTTAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGAtct     <  2:329354/218‑1 (MQ=255)
                                                                                                                                                                                                   tGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCCCGCACAGCATCGGGAAGCGATATTTCAGGGTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCCGAGATAGCGGCAGGCACCgaaaacacg   >  2:656306/1‑222 (MQ=255)
                                                                                                                                                                                                    gCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGAACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGGTAGCGGCCGGGACgggagacacgc  >  2:531365/1‑220 (MQ=255)
                                                                                                                                                                                                         cgcAAACCCCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGca                                                                                                 <  2:26209/129‑1 (MQ=255)
                                                                                                                                                                                                         cgcAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGca                                                                                                 >  1:26209/1‑129 (MQ=255)
                                                                                                                                                                                                                |                                                                                                                                                                                                                        
TCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGCGCCCCGAAGAAAATTAAACAATGTCCAGCGTGCCCAGGCGGTATTTCCGCTGATAGATAAACGAGGTGAGTACAATCGGGATGATGATGGCGATAGCCATTGCCAGCGCAAACACCTGCCAGTAGCTCGGTTGAATCGAGAGAATTCCCGGCAGGCCGCCTACGCCGATGCCATTCGCCATAACGCCGTTCAGGCCGCATAGCAATCCTGCCAGACCAGAACCAATCATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGC  >  NC_000913/4464652‑4465075

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