Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,139,288 C→T D141N (GAC→AAC)  fsaB ← fructose‑6‑phosphate aldolase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,139,2880CT97.8% 148.7 / NA 45D141N (GAC→AAC) fsaBfructose‑6‑phosphate aldolase 2
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (26/18);  total (27/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTC  >  NC_000913/4139056‑4139532
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
gtagtGCCAAATGCGGCATTCCAGGCGGTTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTAc                                                                                                                                                                                                                                    <  2:157012/251‑1 (MQ=255)
      ccAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGtt                                                                                                                                                                                                                              <  1:113952/251‑1 (MQ=255)
        aaaTGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTa                                                                                                                                                                                                                             >  1:347762/1‑250 (MQ=255)
        aaaTGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGGGTTGTGCATTTGTTGCGCTACATCTAAGGGCAGGGAGATGGATTCACTTCCTGCCAGTAACCAGTCCAGCTCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAGCCGAACGAATGCCGTTTCCGCCCTGGGCTTCTCTGCGGTTa                                                                                                                                                                                                                             >  2:26805/1‑250 (MQ=255)
         aaTGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAAc                                                                                                                                                                                                                           >  1:201810/1‑251 (MQ=255)
                 aTTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACgg                                                                                                                                                                                                                     >  2:82088/1‑249 (MQ=255)
                      aGTCGTGTTCGACCGTCTCTATAGCTGACTCTACCGCAGGAGAGTTGAGCATATGATGCGCTACATCTAAGGGCCGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGAGCCTGACGCGGCGTTTTCAAGCCGGCTGCCAGCACCATGATTTCTAGCGCGGGCATTTCTAACAGCAATTGCAACTCCTGAACCGTACGCATACCGATTCCTCCCTTGGCATCAACGCTATTAACATACGAAGCacct                                                                                                                                                                                                              >  2:287678/1‑248 (MQ=255)
                      aGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACg                                                                                                                                                                                                              >  1:182136/1‑251 (MQ=255)
                            gtTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGTGCATCTACGCGGTTAACATACGGAGCAACGTAtttt                                                                                                                                                                                                        <  1:137175/251‑1 (MQ=255)
                                        ctATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAgc                                                                                                                                                                                            >  2:288860/1‑250 (MQ=255)
                                                          cAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTtgt                                                                                                                                                                          >  1:154525/1‑251 (MQ=255)
                                                                  ttGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCGGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGgcgc                                                                                                                                                                   <  1:173430/250‑1 (MQ=255)
                                                                    gAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGATCCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCtat                                                                                                                                                                <  1:173184/251‑1 (MQ=255)
                                                                        aTTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCtatat                                                                                                                                                              <  1:306147/249‑1 (MQ=255)
                                                                        aTTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATAc                                                                                                                                                            >  2:306147/1‑251 (MQ=255)
                                                                                           aGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTgtagt                                                                                                                                          >  1:304966/1‑250 (MQ=255)
                                                                                            gggCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTgtagta                                                                                                                                         >  1:101023/1‑250 (MQ=255)
                                                                                              gCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAAt                                                                                                                                       >  2:119974/1‑250 (MQ=255)
                                                                                                  agggATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGCCGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGGTTGCAGCTCCTGAACCGTACGACTGCCGTTTACGACCTTGGCATCTACGCGGTTAACATACAGAGCACCGAATTTTGCCCCTTGCAGTACCGCTAATCACCCTTGTGCGTCGCTATATACAGCAGTGCCAAGTGTATTAATAACct                                                                                                                                  >  2:43339/4‑251 (MQ=255)
                                                                                                     gATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGAGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGCGCGGCTAACAACCCTTGTGCGGCGCTTTATACAGCAGTGCCAAGTGTAATCATACCCTCtt                                                                                                                               >  2:134944/1‑251 (MQ=255)
                                                                                                      aTGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCCAACAGCGTTTGCAGCTCCCGAACCGTACGAATGCCGTTTCAGCCCTGGGCATCTCCGCGGTAAACATACGGAGCAACGTCTTTTGCCCCTGCCAGTCCGGCTAATACCCATCGTGCGGCGGTAACTACAGCACTGCACAGTGTAGTAATACCctcgtt                                                                                                                              >  2:340665/1‑248 (MQ=255)
                                                                                                         gagtacaaTCCTGCCAGGAAACAGACCAGCGCCTGACGCGGCGTTGTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGACAAGTGTAGTAATACCCTCtttttt                                                                                                                           <  2:111404/244‑1 (MQ=255)
                                                                                                          aTTCACAGCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTc                                                                                                                          <  1:288860/251‑1 (MQ=255)
                                                                                                            tCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCa                                                                                                                         <  1:381549/250‑1 (MQ=255)
                                                                                                                  ccTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGATCCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTAtttt                                                                                                                  <  2:262726/251‑1 (MQ=255)
                                                                                                                  ccTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTAttt                                                                                                                   <  1:80044/250‑1 (MQ=255)
                                                                                                                        aGTAAACGGTCCAGAGCCTGACGCGGCGTTTTAAAGCTGGATGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTACCAGCGTTTGCAGCTCCTGAACCGTACGCATGACGTTTCCGCCCTGGGCATCTACGAGGTTAACATACGGCGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATCACCCTTGGGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATtgc                                                                                                            >  2:426667/1‑251 (MQ=255)
                                                                                                                                       gcCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGACTGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCCAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCtttgg                                                                                              >  2:174899/1‑247 (MQ=255)
                                                                                                                                            aCGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTc                                                                                         >  2:359798/1‑250 (MQ=255)
                                                                                                                                                ggcgTTTTAAAGATGGCTGCCAGCACCATGATTTCTGGCGCGTGCATTTCTACCAGAGTTTGCAGCTCCTGAACCGTACGCATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCACCGTCTTTTGCCCCTGCCCGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCCGTATTTTAATTGCTGCCAGACCTTCGGAAGACACCgg                                                                                    >  2:415902/1‑251 (MQ=255)
                                                                                                                                                     tttAAAGCTGGCTGCCAGCACCATGCTTGCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGAtt                                                                                <  2:201810/250‑1 (MQ=255)
                                                                                                                                                           gctggctgCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTcacc                                                                          >  2:355207/1‑250 (MQ=255)
                                                                                                                                                                  gCCAGCACCATGCTTTCTGGCGCGTGCATTTNTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCAGCTACGCGGTTAACATGCGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTGCGGAAGTCACCGGGATTTTCACCACAATAcc                                                                  <  2:236798/251‑1 (MQ=255)
                                                                                                                                                                         tccctacacgacgctcttccggatCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAAt                                                            <  2:95500/227‑1 (MQ=255)
                                                                                                                                                                           aTGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGACCCGTACGAATGCCGTTTACGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGAGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGCTTTTCACCACAATCCCCGGAATAGc                                                         >  2:132665/1‑251 (MQ=255)
                                                                                                                                                                                ttCTGGCGCGTGAATTTCTAACAGCGTTTTCAGCTCCTGAACCGTACGAATGCCGTTTACGCCCTGGGCATCTACCCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTACTAACCCTTCTGCGGCGCTATATACAGCAGTGCCAAGAGCAGTAATACCCTCTTTTTTCAGAATTTTAATTGCTGCCAGACCTACGGAACACACCGGGATTTTCACCACACTCCCCGGAATACCGTctc                                                    >  2:439573/1‑249 (MQ=255)
                                                                                                                                                                                 tCTGGCGCGTGCATTTCCAACAGCGTTTGCCGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGTAGCATCGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCTGAAGTCACCGGGATTTTCACCACAATACCCTGAATAGCGTcgc                                                    <  2:485780/250‑1 (MQ=255)
                                                                                                                                                                                             aTTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGACCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGGATTTTAATTGCTGCCAGACCTTCGGAAGTCCACGGGATTTTCACCACANTCCCCGGAATCGCGTCGCGCatgagcttagc                                       >  2:217971/1‑241 (MQ=255)
                                                                                                                                                                                                  tAACAGCGTTTGCAGCTCCTGAACCGTACGAATTCCGTTTCCGCCCTGGGCATCTACGCGGTTAACCTACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGcttct                                   <  2:488542/250‑1 (MQ=255)
                                                                                                                                                                                                   aaCAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAAt                                                            >  1:95500/1‑224 (MQ=255)
                                                                                                                                                                                                             gCAGCTCATGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATccct                       >  2:31511/1‑250 (MQ=255)
                                                                                                                                                                                                                gCTCCTGAACCGCACGAATCCCGTTTCCGCCCCGGGCCTCTACGCGGTCCACATACGCAGCAACGTATTTTGCCCCCGCCACTGCGGCTCATAACCCTTGTGCGGCGCTATATACAGCAGTGCCACGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAACTCCCCGGGCTTTTCACCACAATACCCGGCATCGCCTCGCGCAGCCGCTCCGCCTCTTCCACCCTccctgcc                    >  2:336271/1‑247 (MQ=255)
                                                                                                                                                                                                                             aCGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCa       <  1:82088/250‑1 (MQ=255)
                                                                                                                                                                                                                             aCGAATGCCGTTTCCGCCCGGGGCATCTACNCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTGGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCa       <  2:214432/250‑1 (MQ=255)
                                                                                                                                                                                                                                 aTGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTc  >  1:115236/1‑251 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
GTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTC  >  NC_000913/4139056‑4139532

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: