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,6350CT94.9% 158.0 / ‑3.2 60intergenic (‑85/+200)araD/araAL‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (25/31);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

CGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAG  >  NC_000913/66388‑66868
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
cGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATc                                                                                                                                                                                                                                        <  1:430089/251‑1 (MQ=255)
       ggTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCTCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGtt                                                                                                                                                                                                                                 <  1:61932/251‑1 (MQ=255)
       ggTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCAGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGtt                                                                                                                                                                                                                                 >  1:259741/1‑251 (MQ=255)
        gTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTg                                                                                                                                                                                                                                >  2:248620/1‑251 (MQ=255)
          cATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTCGGCCTGATAAGACGCGCCCGCGTCGTATCAGGCGTTGa                                                                                                                                                                                                                               >  2:346080/1‑250 (MQ=255)
          cATGACGCTGTAAGCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTGCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGaa                                                                                                                                                                                                                              <  2:501502/251‑1 (MQ=255)
           aTGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGATGaa                                                                                                                                                                                                                              >  2:430240/1‑250 (MQ=255)
              aCGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGCCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACACAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGcc                                                                                                                                                                                                                          >  1:141194/1‑251 (MQ=255)
                 cTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCgg                                                                                                                                                                                                                        >  1:405245/1‑250 (MQ=255)
                         gACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCttt                                                                                                                                                                                                               <  1:68766/251‑1 (MQ=255)
                               ggAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTAt                                                                                                                                                                                                         >  1:386434/1‑251 (MQ=255)
                                   ggttttgtcatcaaGACTCCGCGCGCTCGATCAACGGCGCTGACGTTGCCACATGTGAGCGTTCCCAGGTTGTGTTTTGGTAGCGCCAGGTTGGCTTCTAATTCCTGGCGTTTGAGATCTTCTAACATGTTGTCTCCTTCGTGCCAGATGCGGTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCTCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCgg                                                                                                                                                                                                     <  2:329595/239‑1 (MQ=21)
                                         aTCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTAc                                                                                                                                                                                                >  2:101023/1‑250 (MQ=255)
                                            acaAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGATGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACaaa                                                                                                                                                                                             >  1:128656/1‑250 (MQ=255)
                                            acaAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACaaa                                                                                                                                                                                             >  1:106437/1‑250 (MQ=255)
                                              ccaGACGCCGAGATCGCGATCAACGGCGCTGCCGTTGCCCCATGAGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCTGATGCGCTTTTCTTATCCGACCTACCAAATCGCAGCTAGTAGGCCTCGTAAGACTCGCCAACGTAGAATCCCGCGTAGAATGCCGGATCCGCTTTGCTTTTCCGGCCTCCAAAATc                                                                                                                                                                                          >  2:384034/3‑251 (MQ=255)
                                              aaaGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATc                                                                                                                                                                                          <  1:53220/251‑1 (MQ=255)
                                                aGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAtagc                                                                                                                                                                                        >  1:33722/1‑248 (MQ=255)
                                                aGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGATTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGc                                                                                                                                                                                        >  2:279993/1‑251 (MQ=255)
                                                   cgctgcccgcgctattaccGGCTCTGACGTTGCCCCATGTTATCGTGTCCATGTTGTGTTTTGGCAGCGCGTGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCGGCGTGCCGGATGCGCTTTGCTTATCCGGACTCCAAAATCGCAGCGTGGAGGCCTGATAAGACGCGCCAGGGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAgc                                                                                                                                                                                     <  2:268638/233‑1 (MQ=18)
                                                   cgccgcGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTGTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATTAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAgc                                                                                                                                                                                     <  1:72724/251‑1 (MQ=255)
                                                       gcgcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGCGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGt                                                                                                                                                                                 >  2:457496/1‑251 (MQ=255)
                                                                aTCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGAt                                                                                                                                                                        <  1:506287/251‑1 (MQ=255)
                                                                       gcgcTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGAcgc                                                                                                                                                                 <  1:248620/251‑1 (MQ=255)
                                                                           tGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGccc                                                                                                                                                             >  1:236893/1‑250 (MQ=255)
                                                                            gACGTTGCCCCATGTGCGCGTGACCAGGTTGTGTTTTGCCGGCGCCAGGTTGGCTTCAAATACCTGGCGTTCGATATCTTCTAACATGTTGACTCCTTCGGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCa                                                                                                                                                             <  2:386434/250‑1 (MQ=21)
                                                                               gTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGt                                                                                                                                                         <  1:126942/251‑1 (MQ=255)
                                                                                 tGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTc                                                                                                                                                        <  1:109129/250‑1 (MQ=255)
                                                                                 tGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTc                                                                                                                                                        <  1:226400/250‑1 (MQ=255)
                                                                                            ggCGTGACCAGGGTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCGGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGGAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGtt                                                                                                                                            <  2:310213/250‑1 (MQ=255)
                                                                                                gggccAGGTGGTGTTCTGGCAGCGGCAGGTTTGCTGCCAATGCCTGGCGTTTGATATCTTCTGACAGGTTGACTCCTTGGTGCCGGATGCGGTTTGCGTGTCGGGTCTACAAAATCGCAGGGTGGAGGCCTGATAAGACGCGCTAGTGTCGTATCAGGCGCTGAGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAAt                                                                                                                                        <  2:42030/248‑1 (MQ=17)
                                                                                                    cAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCg                                                                                                                                    <  1:279993/251‑1 (MQ=255)
                                                                                                             ttttGGCAGCGCCCGGTTGGCTTCTAATACCTGGCGTTTGGGATCTTCTAACATGTTGACTCCTTCGTGCCTGATGCGCTTTGCTTATTCTGCCTACAAAATCGCAGCGTGTAGGCCTGATAAAACGCGCCAGCGTCGCATCAGGCGATGAATGCCGGATGCGCTTTGCTTATCCGGCCCACACAATCACAGCGCGTAGGCCTGCTAAGACGCGCCAGCGTCGTATCCGGCGTTGCATGCCggctggcctt                                                                                                                           >  2:343663/1‑243 (MQ=21)
                                                                                                               gtGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGAGCTTCTAACATGTTGACTCCTTCGTGCCGGATTCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTg                                                                                                                         <  1:414712/250‑1 (MQ=255)
                                                                                                                tGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTg                                                                                                                         <  1:132518/250‑1 (MQ=255)
                                                                                                                   cAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCACCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCtt                                                                                                                      >  1:312052/1‑250 (MQ=255)
                                                                                                                    aGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTAt                                                                                                                    >  1:325211/1‑251 (MQ=255)
                                                                                                                     gcgcCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATc                                                                                                                   <  1:346080/251‑1 (MQ=255)
                                                                                                                      cgcCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTTACTCCTTCGTGCCGGATGCGCGTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCTCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATcc                                                                                                                  <  1:417396/251‑1 (MQ=255)
                                                                                                                            gTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGACTACAAAATCGCAGCGCGTCGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCt                                                                                                             >  2:114088/1‑250 (MQ=255)
                                                                                                                              tGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTAc                                                                                                           <  1:430240/250‑1 (MQ=255)
                                                                                                                                  ttCTAGTACCTGGCGTTTGAGCTCTTCTAACATGTTGCCTAATTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTGGGCCTGATACAACGCGCCAGCGTCTCACCGAGCGTTGAATGCCGGATGCGGTTTGCTTATCCGGCCGACAAAAt                                                                                                      >  2:317325/1‑251 (MQ=18)
                                                                                                                                                   tGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGcc                                                                                       >  1:503525/1‑249 (MQ=255)
                                                                                                                                                       aTCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCtt                                                                                                                           <  1:301729/209‑1 (MQ=255)
                                                                                                                                                         cttctAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATa                                                                                  >  1:170699/1‑248 (MQ=255)
                                                                                                                                                             caactagctGACTCCTTAGTGCGGGATGCGATTTGCTTGTCCGGCTTGCGATAGAGCAGAGTGTAGGCCCGATAAGTCTCGCCGGCGTCGTAGTAGCCGTTGGATGCCGTGTGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCTCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCTGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGAcgcg                                                                           <  2:376476/243‑1 (MQ=9)
                                                                                                                                                                      gACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTc                                                                   <  1:156815/250‑1 (MQ=255)
                                                                                                                                                                                        gcgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATAc                                                 <  1:343663/250‑1 (MQ=255)
                                                                                                                                                                                        gcgcTGTGCTTATCCGGCCTACAAAGTCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTAGCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATAcc                                                <  2:106437/251‑1 (MQ=255)
                                                                                                                                                                                           cTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTCCTTATCCGGCCTCCAAAATCACAGCGTGTAGGCACGATAATCCGCGTCAGcctctaccccggcgttgcctaccgg                                              >  2:263534/1‑226 (MQ=17)
                                                                                                                                                                                            tttGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGAt                                            >  1:105166/1‑251 (MQ=255)
                                                                                                                                                                                                   aTCCGGCCTACAAAATCGCAGCGTGAAAGCCTGATAAGACGCCCCAGCGTCGCATCCGGAGTTGAATGCCGGAAGCGCTTTGCTTATACGGTCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATTATGGATGCGCTTTGCTTATCCGCACTTCAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGACGTTACATACCGGaaccggct                                      >  1:395835/1‑243 (MQ=18)
                                                                                                                                                                                                       ggCCTACAAAATCGCAGCGTGTAGGCCTGATAGGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTa                                 <  1:79802/251‑1 (MQ=255)
                                                                                                                                                                                                        gCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACtt                                  >  1:266980/1‑249 (MQ=255)
                                                                                                                                                                                                         ccTACAAAATCGCAGCGTGTGGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGTATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTGCAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGc                               <  1:208739/251‑1 (MQ=255)
                                                                                                                                                                                                            ccAAAATCGCAGCGTGTAGGCCTGATAAGCCGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGcgac                            <  2:414575/250‑1 (MQ=255)
                                                                                                                                                                                                                         gtgtAGGCCTGATAAGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAAt                >  1:301683/1‑250 (MQ=255)
                                                                                                                                                                                                                           gtatgCCTGATAAGACGCGCCAGCGGCGTAGCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGGCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAATACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAATCCCGTAATaca             <  2:485815/247‑1 (MQ=18)
                                                                                                                                                                                                                           gtAGGCCGGATAAGACGCGCCAGCGCCGTATCAGGCGTTGAATGCCGGATGCGCTTGGCTTAGCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATac              <  2:144751/250‑1 (MQ=255)
                                                                                                                                                                                                                                      aaGACGCGCCAGCGTCGTATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTCCAAAATCGCAGCGTGTAGGCCAGATAAGCCGCGTCACCGTCGCATCAGGCGTTACATACCTGATGCGCCTCGTTATCGAGGAAACCCGTATTACCCTTCTTTcccc  >  2:150858/1‑249 (MQ=17)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
CGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAG  >  NC_000913/66388‑66868

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