Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,336,411 2 bp→AT coding (83‑84/1260 nt) murA ← UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,336,4110GA94.6% 128.2 / ‑0.4 37coding (84/1260 nt)murAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/22);  minor base . (0/2);  total (13/24)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,336,4120AT94.6% 122.4 / ‑3.7 37I28N (ATC→AAC) murAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base T (13/22);  total (13/24)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGAT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCA  >  NC_000913/3336168‑3336654
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                    
aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCa                                                                                                                                                                                                                                              <  2:119575/251‑1 (MQ=255)
    tGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCTCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCagag                                                                                                                                                                                                                                           <  1:187924/250‑1 (MQ=255)
              cAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTa                                                                                                                                                                                                                                <  1:105788/251‑1 (MQ=255)
              cAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTa                                                                                                                                                                                                                                <  1:418882/251‑1 (MQ=255)
                           cAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGTTTTTAACCAGATCGTACTGTGCGCAGAATACATTAACGGCGCGGGCATCAATATGCACAGAACCAGTGCGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGGCTTTCAGTTGCGGGACGTTCTGGATCTCTACCGTTCCTTCCGCCAGTATTGCGGCAAAAAG‑AT‑AAGGCAGAGCAGCATTTTTAGCGCCGTAAAttt                                                                                                                                                                                                                   <  2:284859/251‑2 (MQ=255)
                                gCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACt                                                                                                                                                                                                              <  1:161025/251‑1 (MQ=255)
                                 ccccAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGCCGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACtt                                                                                                                                                                                                             >  1:195539/1‑251 (MQ=255)
                                                    ccACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACTTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCTCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTgt                                                                                                                                                                                          <  2:235564/250‑2 (MQ=255)
                                                           tgggtttATCCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGTTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATGTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTg                                                                                                                                                                                    <  1:429231/246‑1 (MQ=255)
                                                                 tAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCAGTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCGTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGGTGGCCCCTGAACACg                                                                                                                                                                              <  2:44626/250‑1 (MQ=255)
                                                                                  gcgcAGAATCCATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCGACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATGTTTAGCGCCGGAAATTGTGACGTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGt                                                                                                                                                             <  2:223036/250‑1 (MQ=255)
                                                                                     cAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGt                                                                                                                                                         >  1:7638/1‑251 (MQ=255)
                                                                                      aGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGtt                                                                                                                                                        >  1:421651/1‑251 (MQ=255)
                                                                                                                  atGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCg                                                                                                                            >  1:357163/1‑251 (MQ=255)
                                                                                                                      acaGAACCATTACTTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCTACGTCTTTCAGTTTCGGGACGTGCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGTCTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGa                                                                                                                        <  2:172687/251‑1 (MQ=255)
                                                                                                                              aTTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGccc                                                                                                                 >  1:191381/1‑250 (MQ=255)
                                                                                                                               ttACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCAt                                                                                                               >  1:329945/1‑251 (MQ=255)
                                                                                                                                      cTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCATGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCAATAGCTCa                                                                                                         >  1:363227/1‑250 (MQ=255)
                                                                                                                                        aCTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCTAATCGCCCATAGCTCaaaa                                                                                                      <  2:279548/251‑1 (MQ=255)
                                                                                                                                                   cAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGtt                                                                                           <  2:191381/251‑1 (MQ=255)
                                                                                                                                                     tcTGGCTTAGCAGCTCTATTTATGTATTGTCTTTGTTCAGTTTCTGGACGTTTTGGATCTCTACCGGTTCTGCCGCCAGTAGTGCGTCAAAAAG‑ATGA‑GTGAGAGCAACATTGTTAGCGTCGGAAATTTTGACTTCGGCCTTGTTCTGCGTTGGCCCCTTAACACGAAATTGATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCTCTTCCGGCGTATCGCCCTGTGTTCAAAAGCCGTTCAGttt                                                                                          <  2:428726/249‑1 (MQ=255)
                                                                                                                                                         gCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGAt                                                                                     >  1:147104/1‑251 (MQ=255)
                                                                                                                                                            tAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATcg                                                                                   <  1:166152/250‑1 (MQ=255)
                                                                                                                                                                     cATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCGCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCAc                                                                          <  2:174563/250‑1 (MQ=255)
                                                                                                                                                                       ttgctGTATCGACGTTTTTCAGTTTCGAGCCGTTCTGGCTCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCAAATCCGCTACCGGCGAATCGCCCATCGCTCACAAGCCGTTAAGTTTGAGATCGCGCGCCCACTcc                                                                       >  2:405202/5‑251 (MQ=255)
                                                                                                                                                                            gTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAg                                                                   <  2:215975/250‑1 (MQ=255)
                                                                                                                                                                                 gACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGtat                                                             >  1:80487/1‑251 (MQ=255)
                                                                                                                                                                                           aGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGAt                                                   <  1:112191/251‑1 (MQ=255)
                                                                                                                                                                                            gTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATc                                                  <  1:185343/251‑1 (MQ=255)
                                                                                                                                                                                                        tCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAACTTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATg                                       >  2:255800/1‑250 (MQ=255)
                                                                                                                                                                                                           ggATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGCCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATg                                   <  1:101892/251‑1 (MQ=255)
                                                                                                                                                                                                              tctcTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCAt                                                                                                                                                                  >  2:184102/1‑121 (MQ=255)
                                                                                                                                                                                                              tctcTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCAt                                                                                                                                                                  <  1:184102/121‑1 (MQ=255)
                                                                                                                                                                                                                          catcCGCCAGTAGTGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAat                     <  1:136473/248‑1 (MQ=255)
                                                                                                                                                                                                                              cGCCAGAAGTGCGGCAAAAAGATT‑A‑TGCAGAGCAGCATTTTTAGCGCCGGAACTTGTGGCTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGCACTTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCCGCGAATCGCCCCTAGCTCAAAACCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATAtt                >  2:204195/1‑251 (MQ=255)
                                                                                                                                                                                                                                  agtagtGCGGCAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCAt            <  1:438949/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cAAAAAGATT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCa  <  1:444194/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                    
AACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGAT‑A‑GGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCA  >  NC_000913/3336168‑3336654

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