Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,239 G→T A74E (GCG→GAG)  yaaA ← DNA binding and peroxide stress response protein YaaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009136,2390GT89.8% 110.6 / ‑2.5 52A74E (GCG→GAG) yaaADNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base T (29/15);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 6.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTT  >  NC_000913/6003‑6481
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
tgGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGGGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCCTTCTCCGGCTTGGAATCt                                                                                                                                                                                                                                         >  2:609999/1‑251 (MQ=255)
  gTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTgg                                                                                                                                                                                                                                       >  1:580282/1‑251 (MQ=255)
    gctgatATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTTGCt                                                                                                                                                                                                                                     >  2:45848/3‑251 (MQ=255)
        atatCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCa                                                                                                                                                                                                                                <  1:78756/251‑1 (MQ=255)
           tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTc                                                                                                                                                                                                                             >  1:443623/1‑251 (MQ=255)
           tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTc                                                                                                                                                                                                                             >  1:83315/1‑251 (MQ=255)
           tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGACGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCAGGCTGCCCGTc                                                                                                                                                                                                                             >  2:385562/1‑251 (MQ=255)
             cccccAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCAt                                                                                                                                                                                                                           <  1:497737/250‑1 (MQ=255)
              ccccAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCa                                                                                                                                                                                                                            <  1:30001/249‑1 (MQ=255)
                     aTTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAagcg                                                                                                                                                                                                                    >  1:463931/1‑250 (MQ=255)
                       tGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAagcgagc                                                                                                                                                                                                                 >  1:145433/1‑251 (MQ=255)
                         aTACAGATCTTTCCCTCGGGCATTCGTAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCTGGAAAGCATTCGCAAATGCGGTTGGGCAAAATCGAAATCGTCTTCGCGGAAGGTTTCGGTCTGCAAGCCGGTGGAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCGGGCTGCCAGTCATGAAAGCGAgcgg                                                                                                                                                                                                               <  2:618537/251‑1 (MQ=255)
                           aCAGATCTTTCCCTCGGGCATTCTCAAGACGGCTCCCCATTTCCAGACGATATGGCTGCCTTAAATCGAGCGGGCGCAGTCCGCCATCCACGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGCAATCGTCTTCGATGAAGGTTTCGGCCTGCACGCCGGTGTAGCCATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCCGTAAGGAAAGCGAgcggcg                                                                                                                                                                                                             >  2:305372/1‑251 (MQ=255)
                                   ttCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTg                                                                                                                                                                                                         <  1:49905/246‑1 (MQ=255)
                                    ccccTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGGCATGAAAGCGAGCGGCGTTGATAccc                                                                                                                                                                                                   <  2:114642/250‑1 (MQ=255)
                                                           aTCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCAt                                                                                                                                                                            >  1:294486/1‑251 (MQ=255)
                                                            tCCCCATTTCCAGACGATAAGGATGCCTTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGAAAATGCTGTTGGGCACAATCGAAATCGTCTTCGATGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCCTGCAGTCTGGCTGCCAGTCCTGCAAGGGGGCGGCGTTTGAACCCGCCAGtttttggcttatgagcat                                                                                                                                                                            >  2:436236/1‑235 (MQ=255)
                                                                cattgcaaGACGATAAGGCGGCATAAAATCGAGCGGGCGCAGTACCCCATACAAGCCGGAAAGTATTCGCAAATGCGGTTGGGCAAAAGCGAAATCGTCTGCGCGTAAGTTTTCGGCCTGCAAGCCGGTGTATACATCACCTTTAAACGCCATAATCGCCTGGCGGGCATTCTCCGTCGTGAAATCTGGCGGCCAGTCAGGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCg                                                                                                                                                                       <  1:436236/244‑1 (MQ=255)
                                                                 aTTTCCAGACGATAAGAATGCATTAAATCGAGCGGGCGCAGTACTCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGAAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTACGTAGTGTAATCTGGCTGCCAGTCATTAAAGCGCGCGGCATTGATACCCGCCAGTTTGACGCTGATGCGCATccgatt                                                                                                                                                                      >  2:44535/1‑246 (MQ=255)
                                                                         gcaaTAAGGCAGAATTAAAGCGAGCGGGCGCAGTACCCCAACCAGGCCGGAAGCCTTGCGAAAATGCTGTTGGGAAAAATCGAAATCTTCTTCGCTGAGGGTTCCGGCCGGAAAGCCGGTGTAGACATCACGTTTAACCGCCAGAATCGCCTGGCGGGCATTCCCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGGGAGGGGCGTTAATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCt                                                                                                                                                             <  2:145433/248‑1 (MQ=255)
                                                                                     cATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCAACTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGAAGCGCATCAGCGTGCTCATCTGCGGCGGCGTc                                                                                                                                                  >  1:99905/1‑250 (MQ=255)
                                                                                          aaTCGAACGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACACCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGGCATGAAAGAGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAACCTGCGGAGGCGTCAGTTTc                                                                                                                                            >  1:291768/1‑251 (MQ=255)
                                                                                              gagtgggcgcAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCTGCGttt                                                                                                                                                                                                          <  2:636121/183‑2 (MQ=255)
                                                                                              gAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTg                                                                                                                                                                                                         >  1:636121/1‑187 (MQ=255)
                                                                                              aaGCGGACGCAGTACGCCATACAAAACGGAAAGCATTCGCAAATGCTGTTGGGCACAATCGAAATCGTCTTCGCTGAAGGTTTAGGCCTGCAAGCCGGTGTAGACATCACCTTTACACGCCAGAATCGCCTGGCGGGCCTTCTCCGGCGTGAAATCTGGCGGCACGTCCAGAAAGCTCGCGGATTTGATCCCCGCAAGTTTGTGTCTTGTGCGATAAGCATGGCTAATCTGCGGAGGGGTCAGTTTccacg                                                                                                                                        >  2:501103/2‑248 (MQ=255)
                                                                                                 ggggcgcAGTCCGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCgcgctt                                                                                                                                     <  1:45848/250‑3 (MQ=255)
                                                                                                              cctcaacccccctcttccattctgCAAATGTTGTTGGGCAAAACCAAAACCGTCTTCGCTGAAGGTTTCGGCTGCCAAGCCGTTGTAGACATCACCTTAAACCGCCAGAATCGCCTGGCGGGAATTCCCCGGCGTAAAATTTGGCTGCCAGTCATAAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAAct                                                                                                                         <  2:386061/228‑1 (MQ=255)
                                                                                                               cATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGAAAAATCGAAATCGTCTCCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgc                                                                                                                       <  1:659388/251‑1 (MQ=255)
                                                                                                                aTACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgct                                                                                                                      >  1:153083/1‑251 (MQ=255)
                                                                                                                        cGGAAAGCATTCGCAAATGATGTTGGGCAAAATCGAAATCGTCTTAGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCCGTCATGAAGGAGAACGGCGTTGATACCCGCCAATTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAAGTTCCGCGCTTCATGGATAAACGTCTGGGAATTg                                                                                                              >  2:676507/1‑251 (MQ=255)
                                                                                                                           aaaGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTGTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCt                                                                                                           <  1:166533/251‑1 (MQ=255)
                                                                                                                            aaGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCt                                                                                                           >  1:615531/1‑250 (MQ=255)
                                                                                                                             aGCATTCGCAACTGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCAGTGTAGACATCACCTTTAAACGCAAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGGTTACGCGCTTCATGGATCAACTGCTGGGAATTGTCTa                                                                                                          >  1:372901/1‑250 (MQ=255)
                                                                                                                                aTTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAATCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAAc                                                                                                        >  1:631213/1‑249 (MQ=255)
                                                                                                                                    agCAAATGCTGTTAGGCAAAATGGAAATAGTCTTCGCTGAAGGTCTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGACAAAATCGCATGGCGGGCATTCTACGGCGTGACATATGGCTGCAAGTCATGAAAGAGAGGGACGTTGATACCCACAATATTGACGCTGATGCGGAACAACGTGCATATATGAGGAGGAGTCGGGTTACTCGCTTCATGGCACAACGGCTGGGAATGTTCTAACAGCTc                                                                                                   >  2:652653/2‑251 (MQ=255)
                                                                                                                                     gCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAAct                                                                                                                         >  1:386061/1‑228 (MQ=255)
                                                                                                                                     gCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTcc                                                                                                  >  1:396470/1‑251 (MQ=255)
                                                                                                                                        aaTGCTGTTGGGAAAAATCGAAATCGTCTCCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACACCCCCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGttt                                                                                                                                                                                           <  1:374004/160‑1 (MQ=255)
                                                                                                                                          tGCTGTTGGGCAAAATCGAACTCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGACAGCGAGCGGCGTTGATACCCGCCAGTTTTTCGCTGATGCGCATCCGAGTGCTACACTGCGGAGGCGTAAGTTTCAGCGCTTCCTGGATCAACTGATGGGAATTGTTTACCAGCTCCGGccg                                                                                             >  2:634000/1‑249 (MQ=255)
                                                                                                                                                    cAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAgcgc                                                                                   >  1:466424/1‑251 (MQ=255)
                                                                                                                                                     aaaaTCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGcgcg                                                                                  >  1:494217/1‑251 (MQ=255)
                                                                                                                                                         ccgaaatcgTCTCAGCTGAAGGTTCCGGCCTGAAAGCCGGTGTAGACATCACTTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCgtggt                                                                              <  1:664721/250‑1 (MQ=255)
                                                                                                                                                            aaatcgTCGTCACTGAAGGTTTCGGACTGCAAGCCGGTGTAGACATCACCTTTAACCGCCAGAATCGCCTGGCGGGGATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGGGCGGCGTTGATACCCGCCAGTTTGGCGCTGATGCGCATCAGCGTGCTAATATGCGGAGGCGTGAGTTTCCGAGCTTCAAGGATCAAATGAGGGGCATTATCTGACAGCTCGGGCAGCGTAAAGCGCGtggtggg                                                                           >  2:333223/1‑250 (MQ=255)
                                                                                                                                                                     tCGCGGAAGGTTTCGGCCTGCAAGCCGGTGTAGCCCTCACCTTTAAACGCCAGAATCGCCTGGCGGCCATTCTCCGGCGTGAAATCGGGCTGCAAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACgg                                                                     <  1:373925/248‑1 (MQ=255)
                                                                                                                                                                       gCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCttt                                                                >  1:434204/1‑251 (MQ=255)
                                                                                                                                                                           atggTTGCGGCCTGCAAGCCGTCGGAGAAACCCCCTTAAAACGCCAGAATTGCGTGCCGGGCACTCGCCGGCGATAAGTATGGCTGCATGTCATGAAAGCGAGTGTGGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGGGGTGGTCAACGTGCTTTGGTa                                                            <  2:631213/249‑1 (MQ=255)
                                                                                                                                                                               tttCGGCCTGCAGGCCGGTGGAGACATCACCTTTAAACGCAAGAACCGCCTGGCGGGCATGCTCCGGCGTGAAATCGGGCTGCCGGTCATGAAAGCGAGGGGGGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATGGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATc                                                         <  1:305372/250‑1 (MQ=255)
                                                                                                                                                                                   cggctggaaaGCCGGTGTAGACATCACCTTAAAACGCCATAAGCCTCTGGCGGGCATTCTCCGGTGTGAAATCTGGCTGCCCGTCATGAAAGCGGGCGGGGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGc                                                     <  1:126040/243‑1 (MQ=255)
                                                                                                                                                                                              ccGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAgg                                          >  1:109737/1‑251 (MQ=255)
                                                                                                                                                                                                                                ccgggCGGGCATTCTCTGGCGGTAAATCTGGCTGCCAGTCAGGAAAGCGAGCGGTGTTGATACCCGCCAGTTTGGCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAgg        <  1:634000/248‑1 (MQ=255)
                                                                                                                                                                                                                                   ggCGGGCATTCGCCGGCGTGAAGTCGGCTCGCCAGTCATGAAAGCGGGGGGCGTTTACACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCCGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGa       <  1:263979/249‑1 (MQ=255)
                                                                                                                                                                                                                                      gggCATTCCCCGGCGGAAAATCTGGCTGCCAGTTATAAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAAttt  <  1:518608/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
TGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTT  >  NC_000913/6003‑6481

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