Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,635 C→T intergenic (‑85/+200) araD ← / ← araA L‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,6350CT96.8% 177.4 / ‑4.4 62intergenic (‑85/+200)araD/araAL‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (39/21);  total (40/22)
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

GTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGT  >  NC_000913/66389‑66878
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
gTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCa                                                                                                                                                                                                                                                 <  1:117712/251‑1 (MQ=255)
  cAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGACCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCTTGACCATGTTGTGTTTTTGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTTAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAgg                                                                                                                                                                                                                                               <  2:254985/251‑1 (MQ=255)
      ggTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGAGGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGCCGTATCAGGCGtt                                                                                                                                                                                                                                           >  2:286598/1‑251 (MQ=255)
       gTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCAGCCTACAAAATCGCAGCGTGTAGGCCTGTTAAGAGGCGCCATCCTCGTATCAGGCGTTg                                                                                                                                                                                                                                          >  2:299659/1‑251 (MQ=255)
         cATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGa                                                                                                                                                                                                                                         >  1:186665/1‑250 (MQ=255)
                    aagcGACGCCGGAAGGTTTGACCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTGTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCGTCTTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATgcg                                                                                                                                                                                                                             <  2:230660/248‑1 (MQ=255)
                      gcGACGCCGGAAGGGTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCt                                                                                                                                                                                                                           <  2:148259/250‑1 (MQ=255)
                       cGACGCCGGAAGGTTGGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCGTGATAAGACGCGCCAGCGTAGCATCAGGCGTTGAATGCCGGATGCGCtt                                                                                                                                                                                                                          <  2:54997/251‑1 (MQ=255)
                           gCCGGCAGGTTTTATCACAAAGTCGCCGCGCTCGCGATCAACGGCGCTGACGTGGCCCCATGTGAGCGTTACCAGGTTGTTTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCt                                                                                                                                                                                                                      <  2:244064/251‑1 (MQ=37)
                             cGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCtt                                                                                                                                                                                                                     >  1:88462/1‑250 (MQ=255)
                                aaGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATcc                                                                                                                                                                                                                 <  1:329919/251‑1 (MQ=255)
                                            caAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCTGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATTCGCTTTGCTTATCCGTCCTACAAAAt                                                                                                                                                                                                     >  2:357129/1‑251 (MQ=255)
                                             aaaGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCCACAAAATc                                                                                                                                                                                                    >  2:258731/1‑251 (MQ=255)
                                                   gccgcgttccggtgcAACGGCGCTGGCTTTGCCCCATGTGAGCGTGACCAGTTTGGGTTTGTGCAGCGCCAGGTTGGCTTCTAATACCTTGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGTAGGCGCTTTGCTTATCCGGCCTATAAATTCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGTCCTACAAAATCGCAgcg                                                                                                                                                                                              <  2:80549/237‑1 (MQ=21)
                                                      gcgcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGt                                                                                                                                                                                           >  1:297435/1‑251 (MQ=255)
                                                        gcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGt                                                                                                                                                                                           >  2:470369/1‑249 (MQ=255)
                                                           cgcgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGcc                                                                                                                                                                                      >  2:471454/1‑251 (MQ=255)
                                                              gATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTg                                                                                                                                                                                    >  1:326144/1‑250 (MQ=255)
                                                               aTCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTggt                                                                                                                                                                                  >  2:28067/1‑249 (MQ=255)
                                                                  aaCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGa                                                                                                                                                                                                                                 >  2:97052/1‑201 (MQ=255)
                                                                  aaCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGa                                                                                                                                                                                                                                 <  1:97052/201‑1 (MQ=255)
                                                                                tGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCAAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGACTGATAAGACGCGCCAGCGTCg                                                                                                                                                                 >  1:415908/1‑251 (MQ=255)
                                                                                 gCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGgcgc                                                                                                                                                                >  2:294545/1‑251 (MQ=255)
                                                                                     cATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATc                                                                                                                                                             >  1:415926/1‑250 (MQ=255)
                                                                                         tgAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCg                                                                                                                                                        >  1:72231/1‑251 (MQ=255)
                                                                                         tgAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACATGACGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATCAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGCCTGCGCTTTCCTTATCCGGCCTACCAACTCGCATCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCACCAGGCg                                                                                                                                                        >  2:266439/1‑251 (MQ=21)
                                                                                               tGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAAt                                                                                                                                                  >  2:475547/1‑251 (MQ=255)
                                                                                                gACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAAt                                                                                                                                                  <  1:28067/250‑1 (MQ=255)
                                                                                                    aGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCgg                                                                                                                                             >  1:351260/1‑251 (MQ=255)
                                                                                                      gttgtgttgtGTCGGCGCCAGGGTGCCTTCTTATACCTGGCTGTAGATATCTGCTAACATGTTGGCTCCTGCGGGGCGGATGCGCGTTGCTTTTCCGGCGCTCAAAATCGGATGCTGAAGGCCTGATAAGACGTGCCTGCGGCGTATCAGCCCTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGAt                                                                                                                                           <  2:219167/251‑1 (MQ=9)
                                                                                                       ttGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATg                                                                                                                                          >  1:359158/1‑251 (MQ=255)
                                                                                                             tttGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTAAGATCTTCTAACGTGTTATCTCCTTCGTGCCGGTTGCGCTTTGCTTATCCGCACTACAAAATCGAAGCGTGTAGACCTGATACGACACGCCAGCGTCGCATCAGGCGATAGATGCCGGATGCCCTTTGCTTCTCCGGCCTGCAAAATCGCCGCCCGTAGCCCTATTAGGACGACCCAGCTTCGCATAAGGCGTTGCGTGCCGCGTGCGCttt                                                                                                                                    >  1:441273/1‑251 (MQ=255)
                                                                                                                       ccAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGAGGCGCTTTGATTATCCgg                                                                                                                          >  1:458199/1‑251 (MQ=255)
                                                                                                                         aGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGcc                                                                                                                        <  1:475547/251‑1 (MQ=255)
                                                                                                                             tGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGATGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTAc                                                                                                                     >  1:424685/1‑250 (MQ=255)
                                                                                                                              ggCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTaaa                                                                                                                    >  1:312770/1‑248 (MQ=255)
                                                                                                                                     aaTACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGCCGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGCCCTACAAAATCGCAGCGCGGTGGCCAGATACGACGCGCCAGCGTCGCACCAGTAGTTGAAGGCCGGATTCCCTTTGCTTGTCCGGACTAccaacactcc                                                                                                            >  2:147830/1‑242 (MQ=18)
                                                                                                                                              cGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAgg                                                                                                   <  1:470369/251‑1 (MQ=255)
                                                                                                                                                tttGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATa                                                                                                                                                                                 >  1:150228/1‑171 (MQ=255)
                                                                                                                                                tttGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATa                                                                                                                                                                                 <  2:150228/171‑1 (MQ=255)
                                                                                                                                                tttGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGcc                                                                                                 >  1:247584/1‑251 (MQ=255)
                                                                                                                                                 ttGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGACTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGcc                                                                                                 >  1:478197/1‑250 (MQ=255)
                                                                                                                                                        cttctAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATa                                                                                            >  1:279803/1‑248 (MQ=255)
                                                                                                                                                         ttctAACATGTAGACTCCTTCGTGCCGGATGAGCTTTACTTATCCGACCTACAAAATCGCAGCGAGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGACCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGAc                                                                                        >  1:137357/1‑251 (MQ=21)
                                                                                                                                                             aaCATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCATCGTATCAGGCGTTGCATGCAGGATGCGCTTTGCTTACCCGGCCTACAAAATCGCAGCGCGTAGGACTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGt                                                                                    >  1:199723/1‑251 (MQ=255)
                                                                                                                                                               cATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTc                                                                                   >  1:93901/1‑250 (MQ=255)
                                                                                                                                                                  gTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCa                                                                                                >  1:446148/1‑234 (MQ=255)
                                                                                                                                                                  gTTGACGCCTGCGGGCCGGAGGCGCTGTTCGTACCCGGCCTACAAAATCGCAGCGTTTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGGTGAAGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCa                                                                                                <  2:446148/234‑1 (MQ=18)
                                                                                                                                                                     gACTCCGTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCg                                                                            <  1:294545/251‑1 (MQ=255)
                                                                                                                                                                       ccccTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCGGATAAGGCGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCa                                                                          <  2:359158/249‑1 (MQ=255)
                                                                                                                                                                         ccgtCGTGCCGGATGCGCTTGTCTCTTCCGGCCCACAAAACCGTAGCGTGTAGGCCGGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCAt                                                                         <  2:93901/247‑1 (MQ=18)
                                                                                                                                                                         ccTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCAt                                                                         <  1:473699/250‑1 (MQ=255)
                                                                                                                                                                         ccTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATc                                                                        <  1:225790/251‑1 (MQ=255)
                                                                                                                                                                                        cgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCCTTGCTTATCCCGCCTACAAAATCGCCGCGTGTCGGCCAGATAAGGCGCGTCacatgcgcctcaggcgttgaataccg                                                         >  2:169634/1‑226 (MQ=18)
                                                                                                                                                                                               cTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGCCGCGTCACCGGCGCATCAGGCGTTACATCCCGGATGCgg                                                  >  2:257927/1‑251 (MQ=255)
                                                                                                                                                                                                             cAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGAGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAACTCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTCGcgacg                                     >  1:218623/1‑251 (MQ=255)
                                                                                                                                                                                                                 aTCGCAGCGTGTAGGCCTGATAATACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAAc                                 <  2:72231/251‑1 (MQ=255)
                                                                                                                                                                                                                      cgtgggttgggctggatcaGACGCGCCAGCGTCTTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGt                             <  2:186436/233‑1 (MQ=18)
                                                                                                                                                                                                                               ggCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGCTAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGCCGAACCCCGTAATACACtt                    >  1:458031/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aaGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTcc              >  1:294582/1‑249 (MQ=255)
                                                                                                                                                                                                                                             cgcCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACCAAAACGCAGCGTGTAGGCCAGAAAAGCCGCGTCAGCGTCGCCTCAGGCGTTACATACCGGATGCGGCTCCTTATCAACGCAACCCGTAATACCGTTCGTTCcaccgcagg      >  2:291629/1‑250 (MQ=18)
                                                                                                                                                                                                                                                 aGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGATCCAGCGCAGcgcgg  >  1:11203/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
GTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGT  >  NC_000913/66389‑66878

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: