Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,379,266:1 (C)7→9 100% coding (732/1128 nt) zapE ← divisome ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,379,2591.C86.0% 261.4 / 36.5 86coding (739/1128 nt)zapEdivisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (5/7);  new base C (35/39);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,379,2592.C54.7% 157.9 / 155.4 86coding (739/1128 nt)zapEdivisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (13/26);  new base C (27/20);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 3.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTG  >  NC_000913/3379161‑3379348
                                                                                                      ||                                                                                                  
 cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGATGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑c                                                                                               <  1:2222764/100‑1 (MQ=255)
 cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑c                                                                                               <  2:652465/100‑1 (MQ=255)
 cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑c                                                                                               >  2:2106202/1‑100 (MQ=255)
  acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTCATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑cc                                                                                              >  1:2404053/1‑100 (MQ=255)
  acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑cc                                                                                              >  1:1467386/1‑100 (MQ=255)
   caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑ccc                                                                                             >  1:1512789/1‑100 (MQ=255)
    aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑ccca                                                                                            <  2:765530/100‑2 (MQ=255)
       gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑ccc‑‑cc‑‑cc                                                                                     >  1:284425/1‑99 (MQ=255)
         agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑ccc‑‑cc‑‑ccgc                                                                                   <  1:1350388/100‑4 (MQ=255)
         agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTCGCCAATGGCCGATGGTTGATTTCTAACGGCGGTGAATTTTCTCGTT‑‑T‑‑‑C‑‑GCC‑‑CA‑‑C‑‑Gcc                                                                                 <  2:161322/100‑1 (MQ=255)
          gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑ccc‑‑cc‑‑cccgc                                                                                  >  2:1916795/1‑96 (MQ=255)
          gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCAATGTTGCTAATGGCCGATGGTTAATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑cac‑‑cc‑‑ccc                                                                                    >  1:2404219/1‑96 (MQ=255)
           taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑ccc‑‑cc‑‑cccgcc                                                                                 <  1:1855111/100‑6 (MQ=255)
            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTACCGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑ccc‑‑cc‑‑cccgcca                                                                                >  2:492714/1‑94 (MQ=255)
             aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAgc                                                                              <  2:2019138/100‑1 (MQ=255)
             aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑ccc‑‑cc‑‑cccgccag                                                                               >  2:1220254/1‑93 (MQ=255)
             aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CACCC‑‑CC‑‑C‑‑GCCAgc                                                                              <  2:1707432/100‑1 (MQ=255)
               agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAgcgc                                                                            >  2:2367959/1‑100 (MQ=255)
               agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑ccc‑‑cc‑‑cccgccagcg                                                                             <  2:1591206/100‑10 (MQ=255)
               agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑ccc‑‑cc‑‑cccgccagcg                                                                             <  1:2327922/100‑10 (MQ=255)
                gagaCCGCCAGCGTCTGTTTCTCGACGCCCATTGTTGCTAATGTCCGATTGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑ccc‑‑cc‑‑cccgccagcg                                                                             <  2:1595794/99‑10 (MQ=255)
                gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAgcgc                                                                            >  2:438816/1‑100 (MQ=255)
                gagaACGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c‑‑ccc‑‑cc‑‑cccgccagcg                                                                             <  2:707470/99‑10 (MQ=255)
                  gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCa                                                                          >  2:1194426/1‑100 (MQ=255)
                  gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCa                                                                          <  2:2254395/100‑1 (MQ=255)
                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGCTTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAAc                                                                        >  1:587618/1‑100 (MQ=255)
                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAAc                                                                        >  2:225945/1‑100 (MQ=255)
                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAAcc                                                                       <  2:2176935/100‑1 (MQ=255)
                      gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCa                                                                      >  1:1796397/1‑100 (MQ=255)
                      gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGACCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCa                                                                      >  1:1381255/1‑100 (MQ=255)
                       ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGGGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGGCaa                                                                         >  1:2175426/1‑95 (MQ=255)
                       ccAGAGTCTGGTTCTCGACGCCCCTTCTTGTTAATGGCCGATGGTTGATTTCTACCGTCGGTGAATTTTCTCGTT‑‑TCCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAAcc                                                                       <  1:2198094/99‑1 (MQ=255)
                        cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAg                                                                   <  2:211454/100‑1 (MQ=255)
                          gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGt                                                                  >  1:1853249/1‑100 (MQ=255)
                            gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGtt                                                                 <  1:567149/100‑1 (MQ=255)
                             tCTGGTTCTCGCCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGtt                                                                 <  2:266289/97‑1 (MQ=255)
                             tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGACTTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGtt                                                                 <  1:582132/98‑1 (MQ=255)
                              cTGTTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTAt                                                              <  1:1811627/99‑1 (MQ=255)
                              cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTAt                                                              <  2:2117379/100‑1 (MQ=255)
                              cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTAt                                                              <  1:1496770/100‑1 (MQ=255)
                               tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATc                                                             >  2:758785/1‑100 (MQ=255)
                               tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATc                                                             >  1:1964096/1‑100 (MQ=255)
                               tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATcc                                                            <  2:391940/100‑1 (MQ=255)
                               tGGTTCTCGACGCCCATTGTTGCTAAGGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATc                                                             <  2:1967811/100‑1 (MQ=255)
                               tGGTTCTCAACGCCCCTTGTTGCTAATGGCCAATGGTAGATTTCTAACGTCGGTGAATTTGCGCGAG‑‑T‑‑‑C‑‑GCG‑‑CCTGC‑‑GCCAGCGCCAGCCGTAGTTTAt                                                              >  1:1506691/1‑99 (MQ=255)
                                ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTAAATTTTCTCGTT‑‑T‑‑‑C‑‑CCCCACC‑‑C‑‑ACCAGCGGAAACCATTGTatat                                                              >  1:1276322/1‑98 (MQ=255)
                                 gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCa                                                           >  2:520526/1‑100 (MQ=255)
                                   tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCAttt                                                        >  2:1140562/1‑100 (MQ=255)
                                      cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTgc                                                      <  2:1073865/100‑1 (MQ=255)
                                         cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGcc                                                   >  1:1996027/1‑100 (MQ=255)
                                         cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                  >  2:2427690/1‑100 (MQ=255)
                                           cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCCACCATAGTTTATCCATTTGCGCCCgg                                                >  1:428224/1‑100 (MQ=255)
                                            ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                >  2:62787/1‑100 (MQ=255)
                                            ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                               <  2:1258643/100‑1 (MQ=255)
                                              aTTGTTTCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                              <  2:2568359/100‑1 (MQ=255)
                                              aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                              <  2:1796406/100‑1 (MQ=255)
                                                tgttgCTACTGCCCGATGGTTGATTTCTAACGTTGGTGAATTTTCTAGTT‑‑T‑‑‑CCTCCC‑‑CC‑‑C‑‑GCCATCGCCAACCCTAGTTTATCCATTTGCGCCCGGGtt                                             <  1:569758/99‑1 (MQ=255)
                                                tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATGGTTTATCCATTTGCGCCCGGGTTTc                                           >  1:894469/1‑100 (MQ=255)
                                                     cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                       <  2:2180457/100‑1 (MQ=255)
                                                     cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                       <  1:732944/100‑1 (MQ=255)
                                                      tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                      >  1:208372/1‑100 (MQ=255)
                                                          ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                  >  2:1029734/1‑100 (MQ=255)
                                                           gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                 >  2:1375362/1‑100 (MQ=255)
                                                              gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                             <  2:811635/100‑1 (MQ=255)
                                                               gctcttccgatCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                           <  2:1174594/90‑1 (MQ=255)
                                                                tGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                          >  1:1004983/1‑100 (MQ=255)
                                                                  gTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                          >  1:409214/1‑100 (MQ=255)
                                                                  gTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAAACATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                          >  2:1942134/1‑100 (MQ=255)
                                                                  gTTGATTTCTAACGGCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                          <  1:1215323/99‑1 (MQ=255)
                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                        >  2:2539019/1‑100 (MQ=255)
                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGctttcgcc                       <  2:1779908/100‑8 (MQ=255)
                                                                    tGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                        >  2:834300/1‑100 (MQ=255)
                                                                     gATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGctttcgcc                       <  2:1385342/100‑8 (MQ=255)
                                                                     gATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                       <  2:1197218/100‑1 (MQ=255)
                                                                      aTTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                     <  2:2326580/100‑1 (MQ=255)
                                                                       tttCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGctttcgccact                    <  2:895889/100‑11 (MQ=255)
                                                                        ttCTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                    <  2:526393/100‑1 (MQ=255)
                                                                        ttCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                   <  2:2071925/100‑1 (MQ=255)
                                                                        ttCTAACGTCGGTGAATTTTCTCGGTGCC‑‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTGCATCGTGAAGTGGCGAAAGCcac                     <  1:49839/99‑1 (MQ=255)
                                                                          cTAACGTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                  >  2:287309/1‑100 (MQ=255)
                                                                          cTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                 <  1:984042/100‑1 (MQ=255)
                                                                               gTCGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCGGCGCCAACCTTATTTTATCCATTTGCTCCCGGGTTTCATCGTGAGGTGGCGAAAGCCACAGATgc               >  1:88848/1‑98 (MQ=255)
                                                                               gTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCCACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc            <  1:395152/100‑1 (MQ=255)
                                                                               gTCGGTGAATTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc            <  2:1802688/100‑1 (MQ=255)
                                                                                 cGGTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCt           <  1:1527543/100‑1 (MQ=255)
                                                                                 cGGTCAATTTTCTCGCTT‑‑C‑‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTCGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc            <  1:1612231/98‑1 (MQ=255)
                                                                                   gTGAATTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga         >  2:924231/1‑100 (MQ=255)
                                                                                    tGAATTTTCTCGTT‑‑T‑‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg      >  2:162492/1‑100 (MQ=255)
                                                                                      tcttttCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg      <  2:2016294/98‑1 (MQ=255)
                                                                                      aaTTTTCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg      <  2:27050/100‑1 (MQ=255)
                                                                                      aaTTTTCTCGTT‑‑GC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg      <  1:1420972/99‑1 (MQ=255)
                                                                                       aTTTTCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCATAGATGCGCCtgagtgag    >  2:663046/1‑100 (MQ=255)
                                                                                        ttttCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag    >  1:1521742/1‑100 (MQ=255)
                                                                                        ttttCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt   <  1:2367944/100‑1 (MQ=255)
                                                                                         tttCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCCTTTGCGCCCGGGTTTCATCGTGAAGTGGAGAAAGCCACAGATGCGCCTgagtgagt   >  2:640502/1‑100 (MQ=255)
                                                                                         tttCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag    >  1:2016282/1‑99 (MQ=255)
                                                                                         tttCTCGTT‑‑TC‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgggt   <  1:1773136/99‑1 (MQ=255)
                                                                                          ttCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg  >  2:2570379/1‑100 (MQ=255)
                                                                                          ttCTCGTT‑‑TCCC‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg  <  2:1670251/100‑1 (MQ=255)
                                                                                           tctcGTT‑‑T‑‑‑C‑‑CCC‑‑CC‑‑CCAGCCAGCGCCAACCATGATTTAACCAATTCGGCCCGGGTTTCCTCGTGAACTTGCGAAAGCAACAGATGCACCagaggaaga   >  1:1092411/1‑93 (MQ=255)
                                                                                                      ||                                                                                                  
GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑C‑‑CCC‑‑CC‑‑C‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTG  >  NC_000913/3379161‑3379348

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