Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,492,470 C→A G19V (GGC→GTC)  idnT ← L‑idonate/5‑ketogluconate/gluconate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,492,4700CA97.7% 128.0 / NA 43G19V (GGC→GTC) idnTL‑idonate/5‑ketogluconate/gluconate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (20/22);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCCATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGG‑‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCC  >  NC_000913/4492226‑4492699
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
cgcCCATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGTCCCCGAAACCGAGGATCATGGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGGATATAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCTACGACGTCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGatca                                                                                                                                                                                                                                   <  2:626633/251‑1 (MQ=255)
  cccATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGatcatc                                                                                                                                                                                                                                 <  1:579609/251‑1 (MQ=255)
   ccATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGatcatc                                                                                                                                                                                                                                 >  1:362285/1‑250 (MQ=255)
       tGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCAGCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGa                                                                                                                                                                                                                             <  1:587060/250‑1 (MQ=255)
             cgTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATaagaag                                                                                                                                                                                                                      <  2:464981/251‑1 (MQ=255)
             cgTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATaaga                                                                                                                                                                                                                        >  1:603644/1‑249 (MQ=255)
               ttttttACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAagc                                                                                                                                                                                                                     >  1:105125/1‑250 (MQ=255)
                      ccAAAAGTAGCAATCAGCGTCGTGGCGATACGTTGTGCCGCCCCCGTATCAGAAATCAGCTTGCCTAACAGGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTGTAAAGACGATCATCAGGATAAGAAGCAGCGCGAc                                                                                                                                                                                                             <  2:431949/251‑1 (MQ=255)
                      ccAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGAc                                                                                                                                                                                                             <  1:435980/251‑1 (MQ=255)
                           aGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTg                                                                                                                                                                                                        <  1:248931/251‑1 (MQ=255)
                              aGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGcc                                                                                                                                                                                                      <  1:462170/250‑1 (MQ=255)
                                cactCAGCGTAGTGGAGATACGTTGTGCGGCACCCGTTACAGAAATCAGCTGGCCTAACAGTGCCCCGACACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGGGCAGGACGGCCTGTGCAGCCATCCCTCCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCa                                                                                                                                                                                                   <  2:146217/248‑1 (MQ=255)
                                 aaTCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCTCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACTCCTGCCGCaa                                                                                                                                                                                                  <  1:541854/251‑1 (MQ=255)
                                 aaTCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGCAGCAGCGCGACGCCTGCCGCa                                                                                                                                                                                                   >  2:54044/1‑250 (MQ=255)
                                   tCAGCGTAGTGGCGCTACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGCCGGCCTGTGCATCCATCCCTTCGACAAATCCGACGCCGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCaat                                                                                                                                                                                                 >  2:247966/1‑250 (MQ=255)
                                          aGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCCAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGAtt                                                                                                                                                                                          <  2:472475/250‑1 (MQ=255)
                                               cGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGc                                                                                                                                                                                    >  1:226052/1‑251 (MQ=255)
                                               cGATACGTTGAGCCGCACCCGTAACAAAAACCAGCTTGCCAAACATGGCCCCGAAACCGAGGAGCATTGCCAGCCCGCCGAGCGTGCAGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTAAACTTTAAAGACAATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGc                                                                                                                                                                                    <  2:64148/251‑1 (MQ=255)
                                                 aTACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCAt                                                                                                                                                                                  >  1:361306/1‑251 (MQ=255)
                                                              cgcccGTATCAGAAATCATCTTGCCTAACATGGCCCCGAAACCGAGGACCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATctct                                                                                                                                                                     <  1:145310/249‑1 (MQ=255)
                                                                         aaaaTCAGCTTGCCTAACATGGCCCCTAAACCGAGGATCATTGCCAGCCCGCCGAGCGGTCTGCCGATTCCATTTTGTATAGAGTGCAGGTCGGCCTGTGCATCCATCCCGTCGGCAAATCCGACGACGGCATCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATGAATGGCATAGTCTTTATCGCTTATTAGTGCCa                                                                                                                                                          <  1:532366/250‑1 (MQ=255)
                                                                                   tGCCTAACATGGCCCCGAAACCGAGAATCATGGCCAGCCCCCCGAGCGTGCTGCCGATACCATTTTGTAAAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCGGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAagc                                                                                                                                                   <  2:81822/248‑1 (MQ=255)
                                                                                        aaCATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGCGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCACTAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCttt                                                                                                                                           >  1:618590/1‑251 (MQ=255)
                                                                                        aaCATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCttt                                                                                                                                           <  1:54044/251‑1 (MQ=255)
                                                                                                    aaaCCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATCATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGt                                                                                                                               >  1:629039/1‑251 (MQ=255)
                                                                                                                     cAGCCCGCCGAGCGTGCTGCCAGTACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGCCGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCGCTTATGAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAAc                                                                                                             <  2:204694/251‑1 (MQ=255)
                                                                                                                         ccgccgAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGcc                                                                                                         >  1:149714/1‑251 (MQ=255)
                                                                                                                               ggCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCACCCCTTCGGCAAATCCGACGACGGCAGCTACAAGAACGAGGGAAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAAc                                                                                                   <  2:437508/250‑1 (MQ=255)
                                                                                                                                          gATACTATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGCCGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGAAAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATc                                                                                        <  1:2053/251‑1 (MQ=255)
                                                                                                                                           aTACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTGCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATca                                                                                       <  1:296241/251‑1 (MQ=255)
                                                                                                                                             cccATTTTGTAAGGAGGGCAGGACGGCCTGGGCACCCATCCCTTCGGCAAATCCGACGACGGCAGCTCCCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATcaac                                                                                     <  2:140841/250‑1 (MQ=255)
                                                                                                                                                       ataGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGAACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAgg                                                                           >  1:325310/1‑251 (MQ=255)
                                                                                                                                                        taGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGt                                                                          >  1:59093/1‑251 (MQ=255)
                                                                                                                                                           agggCAGGACGGCTGGTGCATCCACCCCCTCGGAAAATCCGACGACGGCAGCTACCAAAACGAGGGCAATAAAGCCTTTAATTTTAAAGACGTTCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGAGTAATGGCAGAGTCTTTATCTCTTATTAGGGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGc                                                                       <  2:534747/248‑1 (MQ=255)
                                                                                                                                                                  gACGGCCTGTGCACCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTAc                                                                <  1:240076/251‑1 (MQ=255)
                                                                                                                                                                     ggCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAAAAGGTGGCCGTTTACGa                                                              >  1:591587/1‑250 (MQ=255)
                                                                                                                                                                                             aaaTCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCCGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGATTTCGCGCTGCACGGAGGGGGGGCTGTAAAACAGCCACTAACATGCCTCCATCAACAAACAACAGTTGGCCGTTTACGAAATCAGAGGCTTTTGAAAAAAGGa                                      >  1:1320/1‑250 (MQ=255)
                                                                                                                                                                                                         ggCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCAcc                         >  2:162806/1‑251 (MQ=255)
                                                                                                                                                                                                          gCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCAcc                         >  1:16498/1‑250 (MQ=255)
                                                                                                                                                                                                           cAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCAGCCGCAATAATGATTAATGGCATAGTCCTCATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCAACAACAAACAACCGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCgcagcaaca                        >  1:380126/1‑247 (MQ=255)
                                                                                                                                                                                                                              aaTAAAGCCGTTAACTTTTAAGACGATCATCAGGATAAGAAGCAACGCGACGCCTGCCGCAATAATTATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGCGGCCGTTTACGAAATCAGAGTCTTTTGCAGAAAGGAACACCGCATCCCCAATCAGTTCCTCCGGATCTcc    >  2:73470/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaaGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTcccc  >  1:66814/1‑250 (MQ=255)
                                                                                                                                                                                                                                  aaGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGAGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTcccc  >  2:358413/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
CGCCCATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGG‑‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCC  >  NC_000913/4492226‑4492699

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