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

TCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGT  >  NC_000913/66543‑67029
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
tctAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGAc                                                                                                                                                                                                                                                    <  1:396057/250‑1 (MQ=255)
           gACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCg                                                                                                                                                                                                                                        <  1:334751/251‑1 (MQ=255)
              tCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACCAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCCGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCAt                                                                                                                                                                                                                                     >  1:238413/1‑251 (MQ=255)
               ccTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCa                                                                                                                                                                                                                                      <  1:560252/249‑1 (MQ=255)
                cTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGGG‑TCGGCCTGATAAGACGCGACAACGTCGccta                                                                                                                                                                                                                                    >  2:702729/1‑247 (MQ=18)
                      gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGtt                                                                                                                                                                                                                             <  1:692805/251‑1 (MQ=255)
                        cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTGTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGAGTCGCATCAGGCGTTAc                                                                                                                                                                                                                           >  2:485226/1‑251 (MQ=255)
                            tGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACacac                                                                                                                                                                                                                       >  1:234386/1‑248 (MQ=255)
                                       tagccgggcgaCAAAATCGCAGCGTGTAGGCCTGAATAGACGCGCCAGCGGCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCga                                                                                                                                                                                                              <  2:8182/239‑2 (MQ=18)
                                          ccgtccTATAAAATCGCGGAGTGTATGCCTGATAAGACGCGCCAGCGTCGCAGCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGGTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCGGCGTCGCATCAGGCGTTACACACCTGATGCGactact                                                                                                                                                                                                         <  2:469276/247‑1 (MQ=17)
                                             gtcTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAGCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACAACTGATGCGACTACTTAg                                                                                                                                                                                                      <  2:393852/249‑1 (MQ=18)
                                              ccTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGc                                                                                                                                                                                                     >  1:502055/1‑251 (MQ=255)
                                                 ccAAAACCGCATCGTGCAGGCCCGCTAAGACGCGCCAGCGTAGCATCAGGCGTTGAACGCCGGATGCGCTTTGCTTATCCGGCCTACCAAAGCGGAACGCGTAGGCC‑CGATAAGATGCGTTAGCGTCGCATCAGGCGTTGAATGTCGGATGCGCTTTGTTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAATACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGcga                                                                                                                                                                                                   <  1:445710/249‑1 (MQ=21)
                                                      accGCAGCGGGTAGGCCTGAAAAGACGCGCCGGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGaaat                                                                                                                                                                                             <  2:112382/249‑2 (MQ=18)
                                                         cgagcttgtggtccggaaaaGACGCGCTATCGGCGCAACATGCGTAGAAGGCCGGGTTCGCTTTGCTTATCCGGCCTTCAAAAACGGACCGCGTAGGCC‑TGATAAGACGCGTCAGCGGCGCATCAGGCGTTGAATGCCGGATGCGCATTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCg                                                                                                                                                                                          <  2:6456/233‑1 (MQ=16)
                                                              gtgtAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAAGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATa                                                                                                                                                                                     <  1:360681/251‑1 (MQ=18)
                                                                gttggCCTGATAAGACGCGCTAGCGGCGCATCAGGCGTTGTATGCCGGTTGCGCTTTGCTTATACGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGGGTCAGCGTCGCATCAGGCGTTTAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCTGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATac                                                                                                                                                                                    <  2:408807/247‑1 (MQ=18)
                                                                    gCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTc                                                                                                                                                                               <  1:145343/251‑1 (MQ=18)
                                                                         aTAAGACGCGCCAGCGTCGAATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACTCCTGATGCGACTACTGAGCGACGAAATCCGTAATACACTGCGTTcc                                                                                                                                                                          <  2:640660/251‑1 (MQ=18)
                                                                                gcgcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTTTCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCcagcgcag                                                                                                                                                                   <  2:429084/251‑1 (MQ=18)
                                                                                            cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGGCTTTaa                                                                                                                                                        >  1:555911/1‑250 (MQ=18)
                                                                                             aTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAAc                                                                                                                                                      <  1:224452/251‑1 (MQ=18)
                                                                                               cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAACTCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGc                                                                                                                                                    >  1:630660/1‑251 (MQ=18)
                                                                                                  gCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTgg                                                                                                                                                 >  1:355268/1‑251 (MQ=18)
                                                                                                   cGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGc                                                                                                                                                >  1:328889/1‑251 (MQ=18)
                                                                                                         attccGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgg                                                                                                                                          <  2:234386/248‑2 (MQ=18)
                                                                                                         aTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgg                                                                                                                                          <  1:467564/251‑2 (MQ=18)
                                                                                                           tccgtattcgcTTTGCTTATCCGGCCTACAAAATCTGAACGCGTAGTCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTTAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGTCTACTTAGCGACGAAATCCGTAATACTCTTCTTTCCAGCGCAGCGCTTCTTTAAACGCTGGCAGGCggg                                                                                                                                         <  1:609674/242‑3 (MQ=18)
                                                                                                             cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATAAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggt                                                                                                                                       >  2:512529/1‑246 (MQ=18)
                                                                                                              ggATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgggtatc                                                                                                                                     <  1:305698/251‑7 (MQ=18)
                                                                                                                aTGCGCGTTGCTTATCCGGCCTAAAAAATCGGAACGCGTATGCC‑GGTTAAGACGCGCATGCGTCGCTTCTGGCGTTTAATGCCGGATTCTCTTTGCTTTTCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCTTCATCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTATCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGTCAGGCGGGTA‑TCGt                                                                                                                                  <  2:447001/251‑1 (MQ=19)
                                                                                                                  gcgcTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCG‑TCGGCCTGATAAGACGCGTCAGAGTCACATCAGGCGTTACACACCAGATGCGACTACTTAGCGACGAAATCCGAAATACACTACAGTCCAGCGCAGCGCGTCCTTAAACGCaggcaggcaggttgcgtt                                                                                                                                  >  2:116528/1‑240 (MQ=17)
                                                                                                                     gttgtcttagcggggcTACAAAATCGGAACGCGTAGGCCGTG‑TAAGACGCGCCCGCGTCGCATCGGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAACGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCGTCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTGAAACGCTGGCAGGCGGGTA‑TCGTTatca                                                                                                                             <  2:247892/236‑1 (MQ=18)
                                                                                                                            tATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAACGTG‑CAGGCCTGGTAAGACGCGTCCGCGTCGCATAAGGCGCTACACACCTGACGCGACGACTTCGCGTCGAAATCCGGAATACGCTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggactcgctacaacaccac                                                                                                                       >  2:169781/1‑231 (MQ=12)
                                                                                                                            tATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTACTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CAc                                                                                                                       <  1:234613/249‑1 (MQ=21)
                                                                                                                             aTCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACCCCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAgggggggatagttatccatccccgg                                                                                                                      >  2:24034/1‑228 (MQ=18)
                                                                                                                             aTCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGt                                                                                                                    <  1:24034/251‑1 (MQ=21)
                                                                                                                                cGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTTTCAAT‑CACCGTggt                                                                                                                 >  1:295393/1‑249 (MQ=21)
                                                                                                                                               gaacgcgTAGGCC‑TGATAAGACGCGCCAGCGTCGCGTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAACATCGCAGCGTG‑TAGGCCTGATAAGACACGTCAGCGGCGCATCAGACGTTACACACCAGATGCGCCTACTTAGCAACGACATCCGTTATACACTTCGTTCCAGCACAGCGCGTCTTTAACAGCTGGCCGGCGGGGA‑TCGGTATCAAT‑CAccacgagtgtcatgttgtgagt                                                                                                  >  2:432537/4‑231 (MQ=12)
                                                                                                                                                 acgcgTAGACC‑TGATAAGACGCGCCAGCGTCGCATGTGGCGTTGGATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGAGATAT‑CCTGCTAAGATGCGTCAGCGACGAATCATGCGTTACACACCTGATGCGAAAAATTATGGACGAACTCCGTCATACACTTCAGTCCAGCGCAGCGCGTCTTTAACCGCTTGCAGGAGGGTG‑GTGTTATACAT‑CACCGTggtgtctatgcagtgacatc                                                                                                >  2:73815/2‑232 (MQ=9)
                                                                                                                                                     gTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGc                                                                                           >  1:280445/1‑251 (MQ=25)
                                                                                                                                                              gATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGgc                                                                                   <  1:485226/251‑1 (MQ=25)
                                                                                                                                                                      gcgcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTGCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCg                                                                           <  1:203741/251‑1 (MQ=25)
                                                                                                                                                                                cGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CAACCTGATTTTCATGTCGTA‑CATCTCGGCGacccggtgcatacacttgtggttc                                                                  >  2:417907/1‑227 (MQ=21)
                                                                                                                                                                                 gaaTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGGCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTAGCAAT‑CACCGGGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAg                                                                <  2:295393/249‑1 (MQ=25)
                                                                                                                                                                                  cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgc                                                               <  1:512529/251‑2 (MQ=37)
                                                                                                                                                                                      aGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCCTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTTGGGTTCAgcgcct                                                           >  1:173148/1‑246 (MQ=37)
                                                                                                                                                                                       ggCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGGGATTTCAATGTCGGG‑CATCTCGGCGAACAGGCGCATATCGTTGAGGGTCAGCGCATg                                                          >  1:637094/1‑251 (MQ=25)
                                                                                                                                                                                       ggCGTTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGCAGCTTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATg                                                          <  1:699048/251‑1 (MQ=255)
                                                                                                                                                                                           ttGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCG‑GAGATCGTTATCAAT‑CACCGTGATTTCAATGTCGGG‑AATCTTAGCGAACCGGCGCataacggttaggatcagcgccaggctg                                                      >  1:380731/1‑227 (MQ=21)
                                                                                                                                                                                                   cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGGTCCGCGCATCGC‑T‑GAAGACggg                                            >  1:428156/1‑250 (MQ=37)
                                                                                                                                                                                                    ggATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGGACCGGCGCATATCGTTGGGGGTCCGCGCATGGC‑T‑GAAGACggg                                            >  1:71588/1‑249 (MQ=25)
                                                                                                                                                                                                          gcTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTgc                                     <  1:31632/251‑1 (MQ=38)
                                                                                                                                                                                                                cTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CACCGTGATTTCCATGTCGTG‑CATTTCGGCGAAATGGCGCATATCGTTGAGGGTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGccc                               >  1:234604/1‑250 (MQ=37)
                                                                                                                                                                                                                 ttATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGaactg                                                                                     <  1:80130/198‑4 (MQ=37)
                                                                                                                                                                                                                   aTCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGGTTTTAAT‑CAACGTGTTTTCAATTTC‑TGTCATCTGGGCGAACTGGCGCCTATCGTTGAGGTTCAGCGCCTGGC‑TTGCAG‑CGGTCTGGTGCGCGccaccc                            >  2:298902/1‑250 (MQ=19)
                                                                                                                                                                                                                      cGGCCTACAAAATCGCAACGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATTAATACACGGT‑ATTTCAATGTCGTG‑CATCTTGGCGTACTCGCGCATATCGGTGGGGGTCAGCGCATGGC‑T‑GGAGACGGGATTGGGCGGGCCACCcgcg                         >  2:347379/1‑251 (MQ=255)
                                                                                                                                                                                                                        gCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACAAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGAG‑TCGGTATCAAA‑TAACGAGATTTCATTATCGTT‑CATTTCGGCGAACTGGCCCCTATAGGTGAGGGTTACGGCATGGC‑T‑GAAGACgggagaggtggtgccacccgca                         >  1:370894/1‑229 (MQ=255)
                                                                                                                                                                                                                            aCAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTAGTTCCAGCGCAGCGCGTCTTTAAACGCTGGAAGGCGGGTA‑TCGTTATCCCA‑CACCGTGTTTTCCATATAGTG‑CATATCGGCGAATTGGCGCATATCGTTTGGGGTCAGCGCCTGTCAT‑GA‑GACGGTATGGGGGGGGCCACCCGCGAGGATc                   >  2:534223/1‑251 (MQ=255)
                                                                                                                                                                                                                               aaagcGCAGCGTG‑TAGGCCAGATAAGACGCGGCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATTGTGCGCGCCACCAGCGAGGATCCAc                <  2:301142/247‑1 (MQ=255)
                                                                                                                                                                                                                               aaaacGCAGCGTG‑TAGGCCTGATAAGACGCGTAAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCAc                <  2:328889/247‑1 (MQ=255)
                                                                                                                                                                                                                               aaaTCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGGATGGTGCGCGCCACCAGCGAGGATCCAc                >  1:484964/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCg          <  2:355268/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aGCGTG‑TAGGCCTGACAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCCTGGC‑T‑GAAGACGGTATGGGGCGCGCCACCAGCGAGGAACCACGCTTCgg         >  1:432032/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gtg‑tAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GGAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGGTTCGGa        >  1:478823/1‑249 (MQ=255)
                                                                                                                                                                                                                                         gtg‑tAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGAATGGTGCGCGCCACCAGCCAGGATCCCCGCTTCGGa        >  1:66449/1‑249 (MQ=255)
                                                                                                                                                                                                                                              aGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGCTTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCCTATCGTTGTGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGGGCGCCCCCAGCGGGGAGCCCCGCTTCGGAAGCAgg  >  2:381247/1‑250 (MQ=255)
                                                                                                                                                                                                                                               ggCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgcg                                                              <  1:586858/192‑3 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
TCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGT  >  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: ‑