Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,784 A→T intergenic (‑234/+51) 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,7840AT98.5% 180.4 / ‑4.5 68intergenic (‑234/+51)araD/araAL‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (33/34);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

TCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGT  >  NC_000913/66543‑67029
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
tctAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGAc                                                                                                                                                                                                                                                    <  1:396057/250‑1 (MQ=255)
           gACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCg                                                                                                                                                                                                                                        <  1:334751/251‑1 (MQ=255)
              tCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACCAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCCGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCAt                                                                                                                                                                                                                                     >  1:238413/1‑251 (MQ=255)
               ccTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCa                                                                                                                                                                                                                                      <  1:560252/249‑1 (MQ=255)
                cTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGGGTCGGCCTGATAAGACGCGACAACGTCGccta                                                                                                                                                                                                                                    >  2:702729/1‑247 (MQ=18)
                      gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGtt                                                                                                                                                                                                                             <  1:692805/251‑1 (MQ=255)
                        cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTGTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGAGTCGCATCAGGCGTTAc                                                                                                                                                                                                                           >  2:485226/1‑251 (MQ=255)
                            tGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACacac                                                                                                                                                                                                                       >  1:234386/1‑248 (MQ=255)
                                       tagccgggcgaCAAAATCGCAGCGTGTAGGCCTGAATAGACGCGCCAGCGGCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCga                                                                                                                                                                                                              <  2:8182/239‑2 (MQ=18)
                                          ccgtccTATAAAATCGCGGAGTGTATGCCTGATAAGACGCGCCAGCGTCGCAGCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGGTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCGGCGTCGCATCAGGCGTTACACACCTGATGCGactact                                                                                                                                                                                                         <  2:469276/247‑1 (MQ=17)
                                             gtcTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAGCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACAACTGATGCGACTACTTAg                                                                                                                                                                                                      <  2:393852/249‑1 (MQ=18)
                                              ccTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGc                                                                                                                                                                                                     >  1:502055/1‑251 (MQ=255)
                                                 ccAAAACCGCATCGTGCAGGCCCGCTAAGACGCGCCAGCGTAGCATCAGGCGTTGAACGCCGGATGCGCTTTGCTTATCCGGCCTACCAAAGCGGAACGCGTAGGCCCGATAAGATGCGTTAGCGTCGCATCAGGCGTTGAATGTCGGATGCGCTTTGTTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAATACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGcga                                                                                                                                                                                                   <  1:445710/249‑1 (MQ=21)
                                                      accGCAGCGGGTAGGCCTGAAAAGACGCGCCGGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGaaat                                                                                                                                                                                             <  2:112382/249‑2 (MQ=18)
                                                         cgagcttgtggtccggaaaaGACGCGCTATCGGCGCAACATGCGTAGAAGGCCGGGTTCGCTTTGCTTATCCGGCCTTCAAAAACGGACCGCGTAGGCCTGATAAGACGCGTCAGCGGCGCATCAGGCGTTGAATGCCGGATGCGCATTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCg                                                                                                                                                                                          <  2:6456/233‑1 (MQ=16)
                                                              gtgtAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAAGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATa                                                                                                                                                                                     <  1:360681/251‑1 (MQ=18)
                                                                gttggCCTGATAAGACGCGCTAGCGGCGCATCAGGCGTTGTATGCCGGTTGCGCTTTGCTTATACGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGGGTCAGCGTCGCATCAGGCGTTTAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCTGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATac                                                                                                                                                                                    <  2:408807/247‑1 (MQ=18)
                                                                    gCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTc                                                                                                                                                                               <  1:145343/251‑1 (MQ=18)
                                                                         aTAAGACGCGCCAGCGTCGAATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACTCCTGATGCGACTACTGAGCGACGAAATCCGTAATACACTGCGTTcc                                                                                                                                                                          <  2:640660/251‑1 (MQ=18)
                                                                                gcgcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTTTCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCcagcgcag                                                                                                                                                                   <  2:429084/251‑1 (MQ=18)
                                                                                            cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGGCTTTaa                                                                                                                                                        >  1:555911/1‑250 (MQ=18)
                                                                                             aTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAAc                                                                                                                                                      <  1:224452/251‑1 (MQ=18)
                                                                                               cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAACTCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGc                                                                                                                                                    >  1:630660/1‑251 (MQ=18)
                                                                                                  gCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTgg                                                                                                                                                 >  1:355268/1‑251 (MQ=18)
                                                                                                   cGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGc                                                                                                                                                >  1:328889/1‑251 (MQ=18)
                                                                                                         attccGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgg                                                                                                                                          <  2:234386/248‑2 (MQ=18)
                                                                                                         aTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgg                                                                                                                                          <  1:467564/251‑2 (MQ=18)
                                                                                                           tccgtattcgcTTTGCTTATCCGGCCTACAAAATCTGAACGCGTAGTCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTTAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGTCTACTTAGCGACGAAATCCGTAATACTCTTCTTTCCAGCGCAGCGCTTCTTTAAACGCTGGCAGGCggg                                                                                                                                         <  1:609674/242‑3 (MQ=18)
                                                                                                             cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATAAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggt                                                                                                                                       >  2:512529/1‑246 (MQ=18)
                                                                                                              ggATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgggtatc                                                                                                                                     <  1:305698/251‑7 (MQ=18)
                                                                                                                aTGCGCGTTGCTTATCCGGCCTAAAAAATCGGAACGCGTATGCCGGTTAAGACGCGCATGCGTCGCTTCTGGCGTTTAATGCCGGATTCTCTTTGCTTTTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCTTCATCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTATCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGTCAGGCGGGTATCGt                                                                                                                                  <  2:447001/251‑1 (MQ=19)
                                                                                                                  gcgcTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTCGGCCTGATAAGACGCGTCAGAGTCACATCAGGCGTTACACACCAGATGCGACTACTTAGCGACGAAATCCGAAATACACTACAGTCCAGCGCAGCGCGTCCTTAAACGCaggcaggcaggttgcgtt                                                                                                                                  >  2:116528/1‑240 (MQ=17)
                                                                                                                     gttgtcttagcggggcTACAAAATCGGAACGCGTAGGCCGTGTAAGACGCGCCCGCGTCGCATCGGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAACGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCGTCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTGAAACGCTGGCAGGCGGGTATCGTTatca                                                                                                                             <  2:247892/236‑1 (MQ=18)
                                                                                                                            tATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAACGTGCAGGCCTGGTAAGACGCGTCCGCGTCGCATAAGGCGCTACACACCTGACGCGACGACTTCGCGTCGAAATCCGGAATACGCTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggactcgctacaacaccac                                                                                                                       >  2:169781/1‑231 (MQ=12)
                                                                                                                            tATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTACTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCAc                                                                                                                       <  1:234613/249‑1 (MQ=21)
                                                                                                                             aTCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACCCCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAgggggggatagttatccatccccgg                                                                                                                      >  2:24034/1‑228 (MQ=18)
                                                                                                                             aTCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGt                                                                                                                    <  1:24034/251‑1 (MQ=21)
                                                                                                                                cGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTTTCAATCACCGTggt                                                                                                                 >  1:295393/1‑249 (MQ=21)
                                                                                                                                               gaacgcgTAGGCCTGATAAGACGCGCCAGCGTCGCGTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAACATCGCAGCGTGTAGGCCTGATAAGACACGTCAGCGGCGCATCAGACGTTACACACCAGATGCGCCTACTTAGCAACGACATCCGTTATACACTTCGTTCCAGCACAGCGCGTCTTTAACAGCTGGCCGGCGGGGATCGGTATCAATCAccacgagtgtcatgttgtgagt                                                                                                  >  2:432537/4‑231 (MQ=12)
                                                                                                                                                 acgcgTAGACCTGATAAGACGCGCCAGCGTCGCATGTGGCGTTGGATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGAGATATCCTGCTAAGATGCGTCAGCGACGAATCATGCGTTACACACCTGATGCGAAAAATTATGGACGAACTCCGTCATACACTTCAGTCCAGCGCAGCGCGTCTTTAACCGCTTGCAGGAGGGTGGTGTTATACATCACCGTggtgtctatgcagtgacatc                                                                                                >  2:73815/2‑232 (MQ=9)
                                                                                                                                                     gTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGc                                                                                           >  1:280445/1‑251 (MQ=25)
                                                                                                                                                              gATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgc                                                                                   <  1:485226/251‑1 (MQ=25)
                                                                                                                                                                      gcgcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTGCAATGTCGTGCATCTCGGCGAACTGGCGCATATCg                                                                           <  1:203741/251‑1 (MQ=25)
                                                                                                                                                                                cGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCAACCTGATTTTCATGTCGTACATCTCGGCGacccggtgcatacacttgtggttc                                                                  >  2:417907/1‑227 (MQ=21)
                                                                                                                                                                                 gaaTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGGCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTAGCAATCACCGGGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAg                                                                <  2:295393/249‑1 (MQ=25)
                                                                                                                                                                                  cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgc                                                               <  1:512529/251‑2 (MQ=37)
                                                                                                                                                                                      aGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCCTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTTGGGTTCAgcgcct                                                           >  1:173148/1‑246 (MQ=37)
                                                                                                                                                                                       ggCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGGGATTTCAATGTCGGGCATCTCGGCGAACAGGCGCATATCGTTGAGGGTCAGCGCATg                                                          >  1:637094/1‑251 (MQ=25)
                                                                                                                                                                                       ggCGTTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGCAGCTTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATg                                                          <  1:699048/251‑1 (MQ=255)
                                                                                                                                                                                           ttGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGAGATCGTTATCAATCACCGTGATTTCAATGTCGGGAATCTTAGCGAACCGGCGCataacggttaggatcagcgccaggctg                                                      >  1:380731/1‑227 (MQ=21)
                                                                                                                                                                                                   cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGGTCCGCGCATCGCTGAAGACggg                                            >  1:428156/1‑250 (MQ=37)
                                                                                                                                                                                                    ggATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGGACCGGCGCATATCGTTGGGGGTCCGCGCATGGCTGAAGACggg                                            >  1:71588/1‑249 (MQ=25)
                                                                                                                                                                                                          gcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTgc                                     <  1:31632/251‑1 (MQ=38)
                                                                                                                                                                                                                cTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTGATTTCCATGTCGTGCATTTCGGCGAAATGGCGCATATCGTTGAGGGTCAGCGCATGGCTGAAGACGGTATGGTGCGCGccc                               >  1:234604/1‑250 (MQ=37)
                                                                                                                                                                                                                 ttATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGaactg                                                                                     <  1:80130/198‑4 (MQ=37)
                                                                                                                                                                                                                   aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGGTTTTAATCAACGTGTTTTCAATTTCTGTCATCTGGGCGAACTGGCGCCTATCGTTGAGGTTCAGCGCCTGGCTTGCAGCGGTCTGGTGCGCGccaccc                            >  2:298902/1‑250 (MQ=19)
                                                                                                                                                                                                                      cGGCCTACAAAATCGCAACGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATTAATACACGGTATTTCAATGTCGTGCATCTTGGCGTACTCGCGCATATCGGTGGGGGTCAGCGCATGGCTGGAGACGGGATTGGGCGGGCCACCcgcg                         >  2:347379/1‑251 (MQ=255)
                                                                                                                                                                                                                        gCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACAAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGAGTCGGTATCAAATAACGAGATTTCATTATCGTTCATTTCGGCGAACTGGCCCCTATAGGTGAGGGTTACGGCATGGCTGAAGACgggagaggtggtgccacccgca                         >  1:370894/1‑229 (MQ=255)
                                                                                                                                                                                                                            aCAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTAGTTCCAGCGCAGCGCGTCTTTAAACGCTGGAAGGCGGGTATCGTTATCCCACACCGTGTTTTCCATATAGTGCATATCGGCGAATTGGCGCATATCGTTTGGGGTCAGCGCCTGTCATGAGACGGTATGGGGGGGGCCACCCGCGAGGATc                   >  2:534223/1‑251 (MQ=255)
                                                                                                                                                                                                                               aaagcGCAGCGTGTAGGCCAGATAAGACGCGGCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATTGTGCGCGCCACCAGCGAGGATCCAc                <  2:301142/247‑1 (MQ=255)
                                                                                                                                                                                                                               aaaacGCAGCGTGTAGGCCTGATAAGACGCGTAAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCAc                <  2:328889/247‑1 (MQ=255)
                                                                                                                                                                                                                               aaaTCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGGATGGTGCGCGCCACCAGCGAGGATCCAc                >  1:484964/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCg          <  2:355268/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aGCGTGTAGGCCTGACAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCCTGGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGAACCACGCTTCgg         >  1:432032/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gtgtAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGGAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGGTTCGGa        >  1:478823/1‑249 (MQ=255)
                                                                                                                                                                                                                                         gtgtAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGAATGGTGCGCGCCACCAGCCAGGATCCCCGCTTCGGa        >  1:66449/1‑249 (MQ=255)
                                                                                                                                                                                                                                              aGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGCTTTCAATGTCGTGCATCTCGGCGAACTGGCGCCTATCGTTGTGGTTCAGCGCATGGCTGAAGACGGTATGGTGGGCGCCCCCAGCGGGGAGCCCCGCTTCGGAAGCAgg  >  2:381247/1‑250 (MQ=255)
                                                                                                                                                                                                                                               ggCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgcg                                                              <  1:586858/192‑3 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
TCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGT  >  NC_000913/66543‑67029

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