Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,219,948 C→T R681H (CGC→CAC)  bglX ← beta‑D‑glucoside glucohydrolase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,219,9480CT97.6% 131.6 / ‑3.9 42R681H (CGC→CAC) bglXbeta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/19);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

AAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCATCTGTTGATTCCAGAACTTCAGCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCATTATACGGGCGACCGGTATTCAGATGGCTGTAGTACACCGGGATCTGCCCGACAGAACGCGG  >  NC_000913/2219704‑2220192
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
aaaCTCTCCTATCTTAACGCGTGCTGTATCAGTGCCGATAAAGACATTGAGCTTGCCAGGCTCGGCGTCATATTTCATCTGTTGATTCCAGATCTTCAGCGCCTCAATATCGATCGTGAAGCGGACAGTCTGAGTTTCGCCTGGTTTGAGAGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGACGACTCATGGAAGCCGTCATATCTTTCAAGGACATCTGCCCTATCGAAGCACCCTCGTGCttcc                                                                                                                                                                                                                                                <  2:416260/251‑3 (MQ=255)
                  gcgTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTAATCTGTTGATTCCAGAACTTCAGCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGAATCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGcc                                                                                                                                                                                                                              >  2:143444/1‑250 (MQ=255)
                            tCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATCTTTCATCTGTTGATTCCAGAACTTCAGCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTACAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACTCTGGCAGt                                                                                                                                                                                                                    >  2:308975/1‑251 (MQ=255)
                                     aTAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCATCTGTTGATTCCAGAACTTCAGCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCGCGCTGGCAGTCACTTTGcc                                                                                                                                                                                                           <  1:379549/251‑1 (MQ=255)
                                                gAACTTGCCAGGCGCGGCGTCATATTTCATCTGTTGATTCCAGAACTTCAGCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTAACGCTTCa                                                                                                                                                                                                <  1:278037/251‑1 (MQ=255)
                                                  aCTTGCCAGGCTCGGCGTCATATTTCATCTGTTGATTCCAGAACTTCAGCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATg                                                                                                                                                                                              <  1:62596/251‑1 (MQ=255)
                                                   cTTGCCAGGCTCGGCGTCATATTTCATCTGTTGATTCCAGAACTTCAGCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATgg                                                                                                                                                                                             >  2:434037/1‑251 (MQ=255)
                                                                gCGTCATATTTCATCTGTTGATTCCAGAACTTCAGCGCCTCAATATCGCTCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTATCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGaa                                                                                                                                                                                  >  1:253026/1‑249 (MQ=255)
                                                                    cATATTTCATCTGTTGATTCCAGAACTTCAGCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGttt                                                                                                                                                                             <  1:117950/250‑1 (MQ=255)
                                                                       atTTCATCTGTTGATTCCAGAACTTCAGCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTcac                                                                                                                                                                          >  1:93478/1‑250 (MQ=255)
                                                                                gTTGATTCCAGACCTTCAGCGCCTCCATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCACAGCCTTTCAGATGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTCCCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTCGCCGTCACGCTTCATGGCCGCCGCAGAAAGTTTCACACCAGAGACg                                                                                                                                                                >  2:275809/1‑251 (MQ=255)
                                                                                   gATTCCAGAACTTCAGCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAACGTTTCACATCAGAGACGGTg                                                                                                                                                             >  2:113078/1‑251 (MQ=255)
                                                                                        agacggtgctcttccgatctatCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCAGGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAgt                                                                                                                                                        <  1:73116/232‑1 (MQ=255)
                                                                                             cTTCAGCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGgtgt                                                                                                                                                   >  1:158973/1‑251 (MQ=255)
                                                                                                        cAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTTCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCa                                                                                                                                        <  1:228738/251‑1 (MQ=255)
                                                                                                                  cgGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCATCCCATAGCCGAACg                                                                                                                              >  1:128659/1‑251 (MQ=255)
                                                                                                                           taCAGTCTGAGTTTCGCCCTGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCGCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGATACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCTAACGGATACAAc                                                                                                                      <  1:288399/249‑1 (MQ=255)
                                                                                                                               gTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCTCATAGCCGAACGGATACAACGcccc                                                                                                                 <  2:93478/251‑1 (MQ=255)
                                                                                                                                         cGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATAAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTtc                                                                                                        >  1:322608/1‑250 (MQ=255)
                                                                                                                                                         gATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCGCCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCGTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTTTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCTAAATAAcgcgta                                                                                        <  2:24217/250‑3 (MQ=255)
                                                                                                                                                                   aaGCCTTTCAGCTGTTTCACAGGGCGGCTCATGGAAGCCTTCACAGCCGGGAAGGCCATCTGCACTACCGTGGCGCCCTCGTGCTTACCGGTGGGCGTCACCTGCACGCTGGCAGGTACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACAACAGAGACGGTGAGAGTGGTGTAGCTCAGCCCATAGCCGAACGCATACATCGCCCCGTTAGCTTCATCAAAGTAACGCGAAGTGTATTTGtt                                                                             <  2:160763/251‑1 (MQ=255)
                                                                                                                                                                      ccTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCgg                                                                          >  2:377709/1‑251 (MQ=255)
                                                                                                                                                                        tttCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCt                                                                        >  1:82763/1‑251 (MQ=255)
                                                                                                                                                                        tttCAGCTGTTTCACAGAGCAACTCATGGGAGCCGACACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGGGCTTACAGGAGATCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGTCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGGGTGTCTCAGCCCATAGCCGACCGGATAAAACGCCCCATTAGCTTAATCAAAATCACACGACGAGGATTTGTTCGGCt                                                                        >  2:93095/1‑251 (MQ=255)
                                                                                                                                                                         ttCAGCTTTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTGTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCt                                                                        <  2:381116/250‑1 (MQ=255)
                                                                                                                                                                                gTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCCAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACCCCCCGTTAGCTTCATCAAAATAACGCGAAGGGTATTTGTTCGGCTTGTCGGc                                                                 >  2:409812/1‑250 (MQ=255)
                                                                                                                                                                                            gACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCATTATACGGGCg                                                     >  1:94746/1‑250 (MQ=255)
                                                                                                                                                                                            gACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGGGTATTTGTTCGGCTTGTCGGCATTATACGGGCg                                                     >  2:310770/1‑250 (MQ=255)
                                                                                                                                                                                                cATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCATTATACGGGCGACCgg                                                <  2:94746/251‑1 (MQ=255)
                                                                                                                                                                                                             acaTCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCATTATACGGGCGACCGGTATTCAGATGGCt                                   <  1:409812/251‑1 (MQ=255)
                                                                                                                                                                                                               aTCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCATTATACGGGCGACCgg                                                <  2:318094/236‑1 (MQ=255)
                                                                                                                                                                                                               aTCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCATTATACGGGCGACCgg                                                >  1:318094/1‑236 (MQ=255)
                                                                                                                                                                                                                 ccTGCAAGTACATCTGCACTACCTTGGCACCCTCGTGCTTACCGGTGTTCGCCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCATTATACGGGCGACCGGTATTCAGATGGCTgta                                <  1:113078/250‑1 (MQ=255)
                                                                                                                                                                                                                 ccTGCAAGAACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTa                                                                                                             >  2:97004/1‑173 (MQ=255)
                                                                                                                                                                                                                 ccTGCAAGAACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTa                                                                                                             <  1:97004/173‑1 (MQ=255)
                                                                                                                                                                                                                     cAAGTACATCTGCACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCCTTCTACGGGCGACCGGTATTCAGATGGCTGTAGTac                            >  2:35485/1‑250 (MQ=255)
                                                                                                                                                                                                                         tACATCTACACTACCGTGGCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCATTATACGGGCGACCGGTATTCAGATGGCTGTAGTACACCgg                       >  2:377928/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tACCGTGGCACCCTCGTGCTTACCGGTGTGCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCATTATACGGGCGCCCGGTATTCAGATGGCTGTAGTACACCGGGATCt                  <  1:438484/245‑1 (MQ=255)
                                                                                                                                                                                                                                    tACCGTGGCACCCTCGTGCTTACCGGTGGTCGTAACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAGAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTAGGCATTATACGGGAGACCGGTATTCCGATGGCTGTAGTACACCGGGATCt                  >  2:438484/1‑245 (MQ=255)
                                                                                                                                                                                                                                           gCACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCATTATACGGGCGACCGGTATTCAGATGGCTGTAGTACACCGGGATCTGCCCGACAGAAcg     >  1:59204/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cACCCTCGTGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTGCGGCTTGTCGGCATTATACGGGCGACCGGTATTCAGATGGCTGTAGTACACCGGGATCTGCCCGACAGAAcgc    <  1:353037/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cccTCGTGCTTACCGGTGGTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCATTATACGGGCGACCGGTATTCAGATGGCTGTAGTACACCGGGATCTGCCCGACAGAACGCgg  <  2:127583/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
AAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCATCTGTTGATTCCAGAACTTCAGCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACGCTTCATGGTCGGCGCAGAAAGTTTCACATCAGAGACGGTGAAAGTGGTGTAGCTCAGCCCATAGCCGAACGGATACAACGCCCCGTTAGCTTCATCAAAATAACGCGAAGTGTATTTGTTCGGCTTGTCGGCATTATACGGGCGACCGGTATTCAGATGGCTGTAGTACACCGGGATCTGCCCGACAGAACGCGG  >  NC_000913/2219704‑2220192

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