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,4700CA96.4% 150.8 / ‑5.7 55G19V (GGC→GTC) idnTL‑idonate/5‑ketogluconate/gluconate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (29/24);  minor base G (2/0);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TAGCGCCCATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGG  >  NC_000913/4492223‑4492713
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
gaGCGCCCAGTGCACGCGTTTTTTACCAAAAGGAGCAATCAGCGTAGTGGCGATACGTTGGGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGGGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGa                                                                                                                                                                                                                                                     <  2:196031/250‑1 (MQ=255)
   cgcCCATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGatc                                                                                                                                                                                                                                                   >  1:55980/1‑250 (MQ=255)
    gcCCATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGatcat                                                                                                                                                                                                                                                 >  1:184211/1‑251 (MQ=255)
    gcCCATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAAAGTGAAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAACCGTTAACTTTCAAGACGaaca                                                                                                                                                                                                                                                  >  2:478620/1‑247 (MQ=255)
                     tttACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAagcagc                                                                                                                                                                                                                                 <  1:601411/250‑1 (MQ=255)
                          cAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACg                                                                                                                                                                                                                           >  1:506807/1‑251 (MQ=255)
                          cAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGACGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACg                                                                                                                                                                                                                           >  1:261005/1‑251 (MQ=255)
                             aaGTAGCAACCAGCTTATTTGCGATACTTTGTGACGCACCCGTATAATAAATCAGCTTTCATAACATGGCCCCGAAACCGAGGATCATTGCAAGCCCGCCGATCGTGCTGCCGATCCCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACTATCATCAGGATAAGAAGCAGCGCTACGcat                                                                                                                                                                                                                        <  2:300820/251‑3 (MQ=255)
                              aGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTg                                                                                                                                                                                                                       >  1:623542/1‑251 (MQ=255)
                                              gTGGCGATACGTTGTGCCGCACCCGTATCAGAAACCAGCTTGCCGAATATGGCCCCGAACCCGAGGATCATTGCCAGCCCGCCGAGCTTGCTGCCGATTCCATTTTGTATAGAGTGCAGGACGGCCTGTTCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTTCCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAATCAGCGCGACGCCTGCTGCAATAATGATTaa                                                                                                                                                                                                       <  1:558865/251‑1 (MQ=255)
                                                       cGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGACGAGCGTGCTGCCGATACCATTTTGTATAGCGTGCAGGCCGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTAACCGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCCGGATAAGAAGCCGCGCGACGCCTGCCGCAATAAAGATTAAttgacaagt                                                                                                                                                                                              >  2:48820/1‑243 (MQ=255)
                                                           gtgCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGTCGCAATAATGATTAATGGCATAGTCtt                                                                                                                                                                                           <  1:587503/250‑1 (MQ=255)
                                                           gtgCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCttt                                                                                                                                                                                          >  1:251601/1‑251 (MQ=255)
                                                             gccgccCCCGTATCAGAAATCAGCTTCCCTAACATGGCCACGAACCCGAGGCTCATTGCTAGCCCGCCGAGCGTGCTGCCGATAAAATTTTATATAGAGTGCAGAACGGCCCGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTCCCAGAACGAGGGCAATACAGCCGTTAACTTTAAAGAAGATAATCAGGATAAGAAGCAGCGCGCCGACTTCCGCagtaatagttaatagcacaagtattt                                                                                                                                                                                         >  2:147209/1‑230 (MQ=255)
                                                              ccGCACCCGTATCAGAAATCAGCATGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGAGCAGGACGGCCTGTGCAGCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGGCGATCATCAGGATAAGAAGCAGCGTGACGCCTGCCGCAATAATGATTAATGGCATAGTCATTAtc                                                                                                                                                                                       >  2:47709/1‑251 (MQ=255)
                                                                gCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGCGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGAGCAGGACGGCCTGTGCATCCAGCCCTTCGGCAAATCCGACGCCGGCAGCTACCAGAACGAGGACAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCAAGACGCCGGCCGCAACAATGATTAATGGCATAGTATTTAtctc                                                                                                                                                                                     >  2:452239/1‑251 (MQ=255)
                                                                   cccGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCtt                                                                                                                                                                                   >  1:240689/1‑250 (MQ=255)
                                                                        aTCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATGGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAg                                                                                                                                                                              <  1:478620/250‑1 (MQ=255)
                                                                          caaaaaTCACCTTTCCTAACTGGGCCCCGAAACCGAGTATTATTGCCAGCCCGCCGAGCGTGCTGCCGATTCCATTTTGTATAGAGTGCAGGACGGCCTGTGCAGCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGc                                                                                                                                                                           <  1:82667/248‑1 (MQ=255)
                                                                            gAAATCAGCTTGCCAAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCGTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCa                                                                                                                                                                         <  1:388086/251‑1 (MQ=255)
                                                                                 cAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGaaa                                                                                                                                                                     <  1:63348/250‑1 (MQ=255)
                                                                                   ccTTGCATAACTTGCCCCCGAAACCTAGGATCATTGCCAGCCCGCCGGGCGGGTTGCCGTTACCATTTTGTATAGAGTGCCGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGTCGGCAGCTCCCATAACGAGGGCAATAAAGCCGTTAACGTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGaaaa                                                                                                                                                                    <  2:465844/248‑1 (MQ=255)
                                                                                    cTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAagc                                                                                                                                                                  >  1:486176/1‑250 (MQ=255)
                                                                                       gCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTGCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAgcg                                                                                                                                                              <  1:48820/251‑1 (MQ=255)
                                                                                            acttGGCCCCAATACCGTGGATAATTGCCACCCCGCCGAGCGTCCTGCCGATACCTTTTTGTATAGAGTGCAGGACGGCCTGTCAATCCATCCCTTCGGCAAATCCGACGACGGCATCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTc                                                                                                                                                         <  2:184211/248‑1 (MQ=255)
                                                                                              aTGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGAc                                                                                                                                                                                                                            <  1:349629/182‑1 (MQ=255)
                                                                                               gggCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTcgcg                                                                                                                                                      <  1:219722/250‑1 (MQ=255)
                                                                                                            gAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAgggg                                                                                                                                        <  1:85108/251‑1 (MQ=255)
                                                                                                                   aTTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGACCGAGGGCAATAAAGACGTTAACTTTAACGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAgggggggctgt                                                                                                                                >  2:385545/1‑246 (MQ=255)
                                                                                                                           cccGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCCCGTAGGGGGGGCTGTTAAACAGc                                                                                                                        >  1:588962/1‑251 (MQ=255)
                                                                                                                            ccgccggtcGTGCTGCCGATACCATTTTGTAAAGAGTGCAGGACGGCCTGTGCACCCATCCCTTCGGCAAATCCAACGACGGCAGCTACCAGAACGAGGGAAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGcc                                                                                                                       <  1:368941/243‑1 (MQ=255)
                                                                                                                                  aGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTTCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAgcg                                                                                                                                                              <  2:499514/208‑1 (MQ=255)
                                                                                                                                  aGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAgcg                                                                                                                                                              >  1:499514/1‑208 (MQ=255)
                                                                                                                                    cGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCATAACGAGGGAAATAAAGCCGTTAACTTAAAAGACGATCAGCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACAt                                                                                                               <  1:58100/251‑1 (MQ=255)
                                                                                                                                        ctgcAGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGAAGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTTGGGGGGGCTGTTAAACAGCCACTAAcctcccc                                                                                                           >  2:228123/1‑245 (MQ=255)
                                                                                                                                          gcCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATgccgcc                                                                                                         >  1:370919/1‑251 (MQ=255)
                                                                                                                                           ccGATACCATTTTGTATAGAGAGCAGGACGGCCTGTGCATCCATCCATTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCACCGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGAGGCGGTTAGACAGCCACTAACaggacgcc                                                                                                         >  2:370662/1‑243 (MQ=255)
                                                                                                                                            caataccattttagatggataccgggACGGCCTGTGCACCCACCCCTCCGGCAATTCCCGCGAGGGCAGATACCAAAAGGAGGGGAATAAGGCCGTTAACTTTAAAGACTATCATCAGTAAAGGAAGCAGCTCGCCGCCTGCCGCAAGAATGATTAATGGCATAGTCTTTATCTCTTATAAGTGCAAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCAt                                                                                                       <  1:310456/227‑1 (MQ=255)
                                                                                                                                                 ccATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCACGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATcaac                                                                                                   >  2:410741/1‑250 (MQ=255)
                                                                                                                                                   aTTTTGTATAGAGTGCAGGACGGCCTGTACATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCCTAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGACGCCATcaacaa                                                                                                 >  2:441152/1‑250 (MQ=255)
                                                                                                                                                     tttGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGGAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaaac                                                                                              >  2:418482/1‑249 (MQ=255)
                                                                                                                                                     tttGTATAGAGTGCAGGACGGCCGGTGCATCCATCCCTTCGGAAAATCCGACGCCGGCAGCTACCAGAACGAGGGAAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaac                                                                                               <  1:67617/250‑1 (MQ=255)
                                                                                                                                                                       cGGCCTGTGCATCCATCCCTTCGGCACATCCGACGCCGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGAAAATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGAATAGGCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATAAACAACCAACAGGGGGCCGtgta                                                                               >  2:43039/1‑245 (MQ=255)
                                                                                                                                                                                     atccCTTCGGCAAATCCCACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGGCGATCATCAGGATAAGAAGCCGCGCGACGCCTGCCGCAATCATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCGGTTCAACAGCCCCTAACATGACGCCATCACCCAACAACAGGTGGCCGTTTTAGCAATCAGAGGCtttt                                                              >  2:547600/1‑251 (MQ=255)
                                                                                                                                                                                                aaaTCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGa                                                    >  1:240282/1‑250 (MQ=255)
                                                                                                                                                                                                      gacgacGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCgc                                             >  1:662924/1‑251 (MQ=255)
                                                                                                                                                                                                         aaaggaagctccaaaaaCGAGGGTAATAAGGCCGTTAACTTTAAAGACGATCATCAGGATAAAAAATACCGCTACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTGATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCAAGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCgcagc                                          <  2:55980/236‑1 (MQ=255)
                                                                                                                                                                                                           cgcaagctcccaaaacaggggCAAAAAACCCGTTATCTTTAAAGACGATCACCAGGATAAAAAGCACCGCGTCCCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCAc                                        <  2:240282/233‑1 (MQ=255)
                                                                                                                                                                                                            gccagcgaccagcccGAGGACAATAAAGCCGTTAACTTTAAAGAAGAGCATCAGAAAAAGACGATGTTCGCCGCCGTCCTCATAAATGTTGAATGGCAAAGTCTTTACCTCTTATTAGTGCCAGAAAAGCAGTGCTTTCGCGCTGCACGTAGGGGGGTCTGTTAAACAGCCACTAACATGCCGTCATCAACAAACAACAGGTGGCCGTTTACGCAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCAcc                                       <  2:214528/237‑1 (MQ=255)
                                                                                                                                                                                                                  tACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGCTAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAACCAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGTAATCAGAGGCTTTTGAAGAAAGGAACACAGCAGCACCAATc                                   >  2:337628/1‑249 (MQ=255)
                                                                                                                                                                                                                                 aaaaaaGCCGTTAACTTTAAAGACGATCATCAGGATAATAACCAGCGCGACCCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGTAAAGCAGCGCTTTCGCGCTGCACg                                                                                                                                               <  2:377846/125‑1 (MQ=255)
                                                                                                                                                                                                                                 aaTAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACg                                                                                                                                               >  1:377846/1‑128 (MQ=255)
                                                                                                                                                                                                                                      aGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGc             >  1:210777/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gTTAACTTTAAAGACGAACATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCAGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGACAGGAACACCGCAGCACCAATCAGTTCATGCGGATCTCCCCAGCGTgc         >  2:281462/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 ttAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCCATCAGTTCCTGCGGATCTCCCCAGCGGGCGCCgggg  >  2:537648/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
TAGCGCCCATTGCACGCGTTTTTTACCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCATTGCCAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGG  >  NC_000913/4492223‑4492713

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