Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,688 G→A intergenic (‑138/+147) 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,6880GA96.4% 141.9 / ‑4.0 57intergenic (‑138/+147)araD/araAL‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (21/33);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

CGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGAT  >  NC_000913/66444‑66921
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          
cgcTCGCGAACAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAAGACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTa                                                                                                                                                                                                                                      <  2:306600/251‑1 (MQ=255)
   ccgcgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGc                                                                                                                                                                                                                                   <  1:353550/250‑1 (MQ=255)
      cgATCAACGGCTCTGACGTTGTCCCAGGTGAGCGGGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATCCCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTTCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑Tg                                                                                                                                                                                                                               <  2:463710/251‑1 (MQ=21)
             cGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGAc                                                                                                                                                                                                                        <  1:247820/251‑1 (MQ=255)
                ggctgtctgtgccctaagtgAGCGGGACCATGTTATGGTTTGGCAGCGCCAGGTGTGCCTCTAATACCTGGCGTTGGAGATCTTCTAACATGTTGCCTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGGGTAGGCTTGATAAGACGCGCCAGCGTCGCAGCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGAcgcg                                                                                                                                                                                                                     <  2:451433/234‑1 (MQ=18)
                       gTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGt                                                                                                                                                                                                              >  1:429084/1‑251 (MQ=255)
                                         gACCAGGTTATGGTTTGGCAGCGCCAGTTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAAt                                                                                                                                                                                             <  1:627843/250‑1 (MQ=21)
                                                  atgtTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATgcg                                                                                                                                                                                   >  1:285378/2‑251 (MQ=255)
                                                       ttGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATCAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGCTGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGACGCATCAGGCGTTGAATGCCGGCTGCGCtttt                                                                                                                                                                              >  2:369377/1‑250 (MQ=21)
                                                        tGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGGCGCATCAGGCGTTGAATGCCGGATGCGCTTTg                                                                                                                                                                              <  1:565022/250‑1 (MQ=255)
                                                               gcCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCg                                                                                                                                                                      >  1:663465/1‑251 (MQ=255)
                                                                 cAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGc                                                                                                                                                                    >  1:162451/1‑251 (MQ=255)
                                                                                   cTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTACCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCgtg                                                                                                                                                  >  1:381160/1‑251 (MQ=255)
                                                                                     ggCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTCCAAAATCGCAGCGCGTAGGCC‑TGATAAGACGCGCCAGCGACGCATCAGGCGTTGAACGCCGGATGCGCTTTGCTTATCCGCCCTACAAAATCGCAGCgtgt                                                                                                                                                 >  1:327372/1‑250 (MQ=255)
                                                                                      gCGTTTGAGATCTTCTAATATGTTGACTGCTTCGTGCCGGATGCGCGTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAg                                                                                                                                               <  1:680267/251‑1 (MQ=255)
                                                                                           tGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGcct                                                                                                                                           <  1:382460/250‑2 (MQ=255)
                                                                                                  ttctAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACCAAATCGGAACGCGGAGGCC‑TGCTACGACGCGCCAGCGTCGCACCAGTCGTTGAATGCCGGAAGCGCTTTGCTTATCCGGCCTACAAAAACGCAgagggtcggcccgaaacgac                                                                                                                                   >  2:514594/1‑232 (MQ=18)
                                                                                                   tctAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGAc                                                                                                                                   <  1:396057/250‑1 (MQ=255)
                                                                                                              gACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCg                                                                                                                       <  1:334751/251‑1 (MQ=255)
                                                                                                                 tCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACCAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCCGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCAt                                                                                                                    >  1:238413/1‑251 (MQ=255)
                                                                                                                  ccTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCa                                                                                                                     <  1:560252/249‑1 (MQ=255)
                                                                                                                   cTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGGGTCGGCCTGATAAGACGCGACAACGTCGccta                                                                                                                   >  2:702729/1‑247 (MQ=18)
                                                                                                                         gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGtt                                                                                                            <  1:692805/251‑1 (MQ=255)
                                                                                                                           cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTGTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGAGTCGCATCAGGCGTTAc                                                                                                          >  2:485226/1‑251 (MQ=255)
                                                                                                                               tGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACacac                                                                                                      >  1:234386/1‑248 (MQ=255)
                                                                                                                                    tttGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCg                                                                                                                                                                      >  1:106798/1‑182 (MQ=255)
                                                                                                                                    ttggcctATCCGGCCGACAAAATCGCAGCGTGTAGGCCGGATAAGGCGCGCCAGCGTCGCATCAGGCGCTGAATGCCGGAGGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCg                                                                                                                                                                      <  2:106798/176‑1 (MQ=17)
                                                                                                                                          tagccgggcgaCAAAATCGCAGCGTGTAGGCCTGAATAGACGCGCCAGCGGCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCga                                                                                             <  2:8182/239‑2 (MQ=18)
                                                                                                                                             ccgtccTATAAAATCGCGGAGTGTATGCCTGATAAGACGCGCCAGCGTCGCAGCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGGTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCGGCGTCGCATCAGGCGTTACACACCTGATGCGactact                                                                                        <  2:469276/247‑1 (MQ=17)
                                                                                                                                                gtcTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAGCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACAACTGATGCGACTACTTAg                                                                                     <  2:393852/249‑1 (MQ=18)
                                                                                                                                                 ccTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGc                                                                                    >  1:502055/1‑251 (MQ=255)
                                                                                                                                                    ccAAAACCGCATCGTGCAGGCCCGCTAAGACGCGCCAGCGTAGCATCAGGCGTTGAACGCCGGATGCGCTTTGCTTATCCGGCCTACCAAAGCGGAACGCGTAGGCC‑CGATAAGATGCGTTAGCGTCGCATCAGGCGTTGAATGTCGGATGCGCTTTGTTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAATACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGcga                                                                                  <  1:445710/249‑1 (MQ=21)
                                                                                                                                                         accGCAGCGGGTAGGCCTGAAAAGACGCGCCGGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGaaat                                                                            <  2:112382/249‑2 (MQ=18)
                                                                                                                                                            cgagcttgtggtccggaaaaGACGCGCTATCGGCGCAACATGCGTAGAAGGCCGGGTTCGCTTTGCTTATCCGGCCTTCAAAAACGGACCGCGTAGGCC‑TGATAAGACGCGTCAGCGGCGCATCAGGCGTTGAATGCCGGATGCGCATTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCg                                                                         <  2:6456/233‑1 (MQ=16)
                                                                                                                                                                 gtgtAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAAGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATa                                                                    <  1:360681/251‑1 (MQ=18)
                                                                                                                                                                   gttggCCTGATAAGACGCGCTAGCGGCGCATCAGGCGTTGTATGCCGGTTGCGCTTTGCTTATACGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGGGTCAGCGTCGCATCAGGCGTTTAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCTGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATac                                                                   <  2:408807/247‑1 (MQ=18)
                                                                                                                                                                       gCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTc                                                              <  1:145343/251‑1 (MQ=18)
                                                                                                                                                                            aTAAGACGCGCCAGCGTCGAATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACTCCTGATGCGACTACTGAGCGACGAAATCCGTAATACACTGCGTTcc                                                         <  2:640660/251‑1 (MQ=18)
                                                                                                                                                                                   gcgcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTTTCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCcagcgcag                                                  <  2:429084/251‑1 (MQ=18)
                                                                                                                                                                                               cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGGCTTTaa                                       >  1:555911/1‑250 (MQ=18)
                                                                                                                                                                                                aTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAAc                                     <  1:224452/251‑1 (MQ=18)
                                                                                                                                                                                                  cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAACTCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGc                                   >  1:630660/1‑251 (MQ=18)
                                                                                                                                                                                                     gCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTgg                                >  1:355268/1‑251 (MQ=18)
                                                                                                                                                                                                      cGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGc                               >  1:328889/1‑251 (MQ=18)
                                                                                                                                                                                                            attccGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgg                         <  2:234386/248‑2 (MQ=18)
                                                                                                                                                                                                            aTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgg                         <  1:467564/251‑2 (MQ=18)
                                                                                                                                                                                                              tccgtattcgcTTTGCTTATCCGGCCTACAAAATCTGAACGCGTAGTCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTTAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGTCTACTTAGCGACGAAATCCGTAATACTCTTCTTTCCAGCGCAGCGCTTCTTTAAACGCTGGCAGGCggg                        <  1:609674/242‑3 (MQ=18)
                                                                                                                                                                                                                cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATAAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggt                      >  2:512529/1‑246 (MQ=18)
                                                                                                                                                                                                                 ggATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgggtatc                    <  1:305698/251‑7 (MQ=18)
                                                                                                                                                                                                                   aTGCGCGTTGCTTATCCGGCCTAAAAAATCGGAACGCGTATGCC‑GGTTAAGACGCGCATGCGTCGCTTCTGGCGTTTAATGCCGGATTCTCTTTGCTTTTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCTTCATCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTATCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGTCAGGCGGGTATCGt                  <  2:447001/251‑1 (MQ=19)
                                                                                                                                                                                                                     gcgcTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTCGGCCTGATAAGACGCGTCAGAGTCACATCAGGCGTTACACACCAGATGCGACTACTTAGCGACGAAATCCGAAATACACTACAGTCCAGCGCAGCGCGTCCTTAAACGCaggcaggcaggttgcgtt                 >  2:116528/1‑240 (MQ=17)
                                                                                                                                                                                                                        gttgtcttagcggggcTACAAAATCGGAACGCGTAGGCCGTG‑TAAGACGCGCCCGCGTCGCATCGGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAACGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCGTCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTGAAACGCTGGCAGGCGGGTATCGTTatca             <  2:247892/236‑1 (MQ=18)
                                                                                                                                                                                                                               tATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAACGTGCAGGCCTGGTAAGACGCGTCCGCGTCGCATAAGGCGCTACACACCTGACGCGACGACTTCGCGTCGAAATCCGGAATACGCTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggactcgctacaacaccac      >  2:169781/1‑231 (MQ=12)
                                                                                                                                                                                                                               tATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTACTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCAc        <  1:234613/249‑1 (MQ=21)
                                                                                                                                                                                                                                aTCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACCCCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAgggggggatagttatccatccccgg     >  2:24034/1‑228 (MQ=18)
                                                                                                                                                                                                                                aTCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGt     <  1:24034/251‑1 (MQ=21)
                                                                                                                                                                                                                                   cGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTTTCAATCACCGTggt  >  1:295393/1‑249 (MQ=21)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          
CGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGAT  >  NC_000913/66444‑66921

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: ‑