Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 95,418 T→G S257A (TCG→GCG)  murF → UDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091395,4180TG99.5% 764.5 / ‑5.1 216S257A (TCG→GCG) murFUDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (106/109);  total (106/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 216 total aligned reads displayed.

TGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATGCGCTGGCAG  >  NC_000913/95323‑95514
                                                                                               |                                                                                                
tGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCa                                                                                              >  1:4103468/1‑100 (MQ=255)
 gAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCAc                                                                                             <  1:1775257/100‑1 (MQ=255)
 gAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCAc                                                                                             <  1:2846025/100‑1 (MQ=255)
  aaCTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACg                                                                                            >  1:4584927/1‑100 (MQ=255)
    cTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCGATGTGACCGCGCACGGt                                                                                          >  1:4713690/1‑100 (MQ=255)
       gCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACg                                                                                       <  1:4603122/100‑1 (MQ=255)
        cAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACgg                                                                                      <  1:3254501/100‑1 (MQ=255)
        cAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACgg                                                                                      <  2:1044141/100‑1 (MQ=255)
         agagCGTAATTGGCGCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGa                                                                                     <  1:5154750/100‑1 (MQ=255)
           agCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAAt                                                                                   <  1:3276952/100‑1 (MQ=255)
            gCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAAtt                                                                                  >  1:3798099/1‑100 (MQ=255)
             cGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCAGCCAATGTCGTCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAAttt                                                                                 <  2:5246829/100‑1 (MQ=255)
             cGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATACATGTGACCGCGCACGGTACGGAAttt                                                                                 >  2:3375734/1‑100 (MQ=255)
              gTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTa                                                                                >  1:3339625/1‑100 (MQ=255)
               tAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTAc                                                                               >  2:5153639/1‑100 (MQ=255)
               tAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTAc                                                                               <  1:5418673/100‑1 (MQ=255)
               tAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAAAGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTAc                                                                               >  1:4003781/1‑100 (MQ=255)
                aaTTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTAcc                                                                              <  2:2393197/100‑1 (MQ=255)
                aaTTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGACAACAGCGATTTCACCGTCACCAATATCAATGTGACCGCGCACGGTAAGGAATTTaga                                                                              >  2:4710002/1‑98 (MQ=255)
                 aTTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATTTGACCGCGCACGGTACGGAATTTAccc                                                                             <  1:3924788/100‑1 (MQ=255)
                  ttGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCt                                                                            <  2:4535214/100‑1 (MQ=255)
                  ttGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCt                                                                            <  2:3496693/100‑1 (MQ=255)
                   tGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTa                                                                           <  2:4027797/100‑1 (MQ=255)
                   tGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTa                                                                           >  2:1531794/1‑100 (MQ=255)
                    ggCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGGCCGCGCACGGTACGGAATTTACCCTAc                                                                          >  2:2643296/1‑100 (MQ=255)
                    ggCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTAc                                                                          >  2:1000484/1‑100 (MQ=255)
                    gccTCACGAAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTAc                                                                          <  1:1649344/98‑1 (MQ=255)
                       tCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACaaa                                                                       >  1:4661634/1‑100 (MQ=255)
                        cACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAAc                                                                      <  2:3763241/100‑1 (MQ=255)
                          cGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAAccc                                                                    >  2:3161271/1‑100 (MQ=255)
                          cGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAAccc                                                                    >  2:2916008/1‑100 (MQ=255)
                           gCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAAcccc                                                                   <  1:2756099/100‑1 (MQ=255)
                           gCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAAcccc                                                                   >  1:2444241/1‑100 (MQ=255)
                            cAAAGTGTGGCGTTTCTCGCCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGAACGGTACGGAATTTACCCTACAAACCCCa                                                                  <  2:4437129/100‑1 (MQ=255)
                              aaGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAAc                                                                <  1:3727432/100‑1 (MQ=255)
                               aGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAAcc                                                               <  1:5468622/100‑1 (MQ=255)
                               aGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAAcc                                                               <  1:1413878/100‑1 (MQ=255)
                                 tgtgGCGTTTCTCACCCAGTGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCgg                                                             <  1:3018971/100‑1 (MQ=255)
                                 tgtgGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCgg                                                             >  1:3611509/1‑100 (MQ=255)
                                 tgtgGCGTTGCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCgg                                                             <  2:2887176/100‑1 (MQ=255)
                                  gtgGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGt                                                            >  2:303758/1‑100 (MQ=255)
                                   tgGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTa                                                           <  1:5373994/100‑1 (MQ=255)
                                    ggCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAg                                                          >  1:2923527/1‑100 (MQ=255)
                                    ggCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAg                                                          >  1:1823885/1‑100 (MQ=255)
                                    ggCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAg                                                          >  2:3737184/1‑100 (MQ=255)
                                    ggCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGCTAg                                                          >  2:1311483/1‑100 (MQ=255)
                                    ggCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATCTACCCTACAAACCCCAACCGGTAg                                                          >  2:1613958/1‑100 (MQ=255)
                                     gCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGc                                                         >  1:1585433/1‑100 (MQ=255)
                                     gCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGc                                                         <  2:4890294/100‑1 (MQ=255)
                                      cGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCg                                                        >  1:3841488/1‑100 (MQ=255)
                                      cGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCg                                                        >  1:912531/1‑100 (MQ=255)
                                      cGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCg                                                        <  2:1440021/100‑1 (MQ=255)
                                       gTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGt                                                       >  2:2731079/1‑100 (MQ=255)
                                         ttCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCg                                                     <  1:2579359/100‑1 (MQ=255)
                                          tctcACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGa                                                    >  2:4362100/1‑100 (MQ=255)
                                           ctcACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCTACCGGTAGCGTCGAt                                                   >  2:1934143/1‑100 (MQ=255)
                                           ctcACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGAt                                                   <  2:495391/100‑1 (MQ=255)
                                           ctcACCCAATGCCGCCAACAGCGATTTCACCGACACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCGACAAGCCCCAACCGGTAGCGTCGAt                                                   <  2:3653220/100‑1 (MQ=255)
                                               cccAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTc                                               <  1:2061572/100‑1 (MQ=255)
                                               cccAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTc                                               >  2:2536893/1‑100 (MQ=255)
                                                ccAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTct                                              <  1:1911305/100‑1 (MQ=255)
                                                ccAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAACGTCGATGTTct                                              <  1:308974/100‑1 (MQ=255)
                                                 cAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTctg                                             <  2:2961475/100‑1 (MQ=255)
                                                  aaTGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTctgc                                            <  1:3858484/100‑1 (MQ=255)
                                                    ggccgccAACTGCGATGCCACCGGCACCCATATCCATGTGACCGCGCACGGTACGGAAGTCACCCTACAAACCCCAACCGGTAGCGTCGATGTTcggctg                                          <  2:4653161/99‑1 (MQ=255)
                                                     gccgccAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCtgctgc                                         >  1:4765380/1‑100 (MQ=255)
                                                     gccgccAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCtgctgc                                         >  1:4128038/1‑100 (MQ=255)
                                                     gccgccAACAGCAATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGAAGCGTCGATGTTCtgctgc                                         <  2:2483733/100‑1 (MQ=255)
                                                      ccgccAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGcc                                        >  1:3121087/1‑100 (MQ=255)
                                                       cgccAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCg                                       <  1:643141/100‑1 (MQ=255)
                                                         ccAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTCCCGtt                                     >  2:705250/1‑100 (MQ=255)
                                                           naCAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGc                                   >  1:341714/2‑100 (MQ=255)
                                                             cAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCg                                 >  2:559273/1‑100 (MQ=255)
                                                             cAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCg                                 >  2:5010221/1‑100 (MQ=255)
                                                              aGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCgg                                <  2:4213251/100‑1 (MQ=255)
                                                               gCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCggg                               >  2:5330803/1‑100 (MQ=255)
                                                               gCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCggg                               <  1:1150503/100‑1 (MQ=255)
                                                               gCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCggg                               <  1:5396748/100‑1 (MQ=255)
                                                               gCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCggg                               <  2:3884718/100‑1 (MQ=255)
                                                                cGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCgggg                              <  1:4523272/100‑1 (MQ=255)
                                                                  aTTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCg                            >  1:2721314/1‑100 (MQ=255)
                                                                  aTTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCg                            <  2:3960859/100‑1 (MQ=255)
                                                                   tttCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGt                           >  1:176241/1‑100 (MQ=255)
                                                                          gCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAata                    >  1:4342664/1‑100 (MQ=255)
                                                                             accaATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTg                 >  1:3195344/1‑100 (MQ=255)
                                                                             accaATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTg                 <  1:2257300/100‑1 (MQ=255)
                                                                             accaATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTg                 <  1:532201/100‑1 (MQ=255)
                                                                              ccaATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCACCCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGc                <  2:5414112/100‑1 (MQ=255)
                                                                               caATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCg               <  1:1451540/100‑1 (MQ=255)
                                                                                 atatCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGaa             <  1:2316574/100‑1 (MQ=255)
                                                                                   atCCATGTGACCGCGCACGGTACGGAATTGACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATg           <  1:5088098/100‑1 (MQ=255)
                                                                                    tacATGTGACCGCGCAAGGTACGGAATTTAACCTCCAAACCCCAAAGGGTAGCGTCGGTGTTCTGCTGCCGTTGCCGGGGCGTCACAAAATTGCACATgc          >  2:118432/3‑100 (MQ=255)
                                                                                    atcttgtgACCGCGCACGGTACGGAATTTACCCTGCAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATgc          <  2:5345528/96‑1 (MQ=255)
                                                                                      cATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATgcgc        >  1:5518353/1‑100 (MQ=255)
                                                                                      cATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATgcgc        >  2:2949504/1‑100 (MQ=255)
                                                                                       aTGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATGCGCt       <  1:716838/100‑1 (MQ=255)
                                                                                        tgtgACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATGCGCTg      <  1:4013634/100‑1 (MQ=255)
                                                                                         gtgACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATGCGCTgg     <  1:2800618/100‑1 (MQ=255)
                                                                                           gACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATGCGCTGGCa   <  2:2450478/100‑1 (MQ=255)
                                                                                           gACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATGCGCTGGCa   >  2:4059007/1‑100 (MQ=255)
                                                                                               |                                                                                                
TGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATGCGCTGGCAG  >  NC_000913/95323‑95514

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