Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 701,381 (C)5→4 coding (214/1221 nt) nagC ← N‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913701,3770C.98.4% 291.2 / ‑3.5 62coding (218/1221 nt)nagCN‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (35/26);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCG  >  NC_000913/701231‑701522
                                                                                                                                                  |                                                                                                                                                 
gggTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGccccg                                                                                                                                               >  2:13827/1‑150 (MQ=255)
  gTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGccccggt                                                                                                                                             <  1:261801/151‑4 (MQ=255)
    aCGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGccccggtg                                                                                                                                            <  1:35503/150‑5 (MQ=255)
      gTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTgg                                                                                                                                          >  2:601882/1‑149 (MQ=255)
          cGGCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCct                                                                                                                                    >  2:980307/1‑151 (MQ=255)
              aGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCctgct                                                                                                                                 >  2:787245/1‑150 (MQ=255)
                 gggTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATc                                                                                                                             <  1:722036/151‑1 (MQ=255)
                   gTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCa                                                                                                                            >  2:169558/1‑150 (MQ=255)
                    tAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCaa                                                                                                                           >  2:31158/1‑150 (MQ=255)
                      aTGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAActt                                                                                                                        <  2:227625/151‑1 (MQ=255)
                      aTGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAActt                                                                                                                        >  1:391823/1‑151 (MQ=255)
                         ttcttcTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAActtctt                                                                                                                     >  2:993085/1‑151 (MQ=255)
                         ttcttcTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAActtctt                                                                                                                     >  1:935618/1‑151 (MQ=255)
                             tctGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCttt                                                                                                                    <  2:1060265/148‑1 (MQ=255)
                             tctGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCttt                                                                                                                    <  1:491665/148‑1 (MQ=255)
                             tctGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTg                                                                                                                   <  1:791917/149‑1 (MQ=255)
                              ctGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATc                                                                                                                >  2:97997/1‑151 (MQ=255)
                                gCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCa                                                                                                               <  1:52856/151‑1 (MQ=255)
                                gCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGAt                                                                                                                 <  2:1059343/148‑1 (MQ=255)
                                 ccAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGc                                                                                                             <  1:4746/151‑1 (MQ=255)
                                     cACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCcgc                                                                                                         >  2:479765/1‑151 (MQ=255)
                                       cTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCcgcg                                                                                                        >  1:925145/1‑150 (MQ=255)
                                        tttGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGtt                                                                                                      >  1:472380/1‑151 (MQ=255)
                                        tttGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGtt                                                                                                      >  2:530136/1‑151 (MQ=255)
                                           gctgctGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCg                                                                                                    <  2:939911/150‑1 (MQ=255)
                                           gctgctGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCg                                                                                                    >  2:887937/1‑150 (MQ=255)
                                             tgctgAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGAt                                                                                                  >  2:696237/1‑150 (MQ=255)
                                              gctgAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGAt                                                                                                  <  2:45650/149‑1 (MQ=255)
                                              gctgAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGAt                                                                                                  <  1:973014/149‑1 (MQ=255)
                                                tgAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGc                                                                                              >  1:59975/1‑151 (MQ=255)
                                                      cAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGAcgc                                                                                        >  2:943814/1‑151 (MQ=255)
                                                      cAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGAcg                                                                                         <  2:895401/150‑1 (MQ=255)
                                                         aCAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGt                                                                                      >  1:366231/1‑150 (MQ=255)
                                                           agagTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAAt                                                                                   >  1:190349/1‑151 (MQ=255)
                                                             agTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAAttt                                                                                 >  2:124151/1‑151 (MQ=255)
                                                                gATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTg                                                                               <  2:774726/150‑1 (MQ=255)
                                                                 aTGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGt                                                                             <  1:192184/151‑1 (MQ=255)
                                                                      ggCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGc                                                                        >  1:798841/1‑151 (MQ=255)
                                                                        cGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTg                                                                      >  1:721693/1‑151 (MQ=255)
                                                                         gTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTg                                                                      <  2:387164/150‑1 (MQ=255)
                                                                               acgacCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCggggg                                                               <  2:441688/151‑1 (MQ=255)
                                                                                  acCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCaa                                                            <  2:951965/151‑1 (MQ=255)
                                                                                        ccGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGc                                                      >  1:565128/1‑151 (MQ=255)
                                                                                           tACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGAGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCttt                                                   >  2:506468/1‑151 (MQ=255)
                                                                                            aCGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTg                                                  >  2:1007157/1‑151 (MQ=255)
                                                                                              gCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGc                                                      >  1:789453/1‑145 (MQ=255)
                                                                                              gCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCt                                                >  2:559395/1‑151 (MQ=255)
                                                                                               ccGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTc                                               >  2:557003/1‑151 (MQ=255)
                                                                                                 gATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCgg                                             >  1:761950/1‑151 (MQ=255)
                                                                                                  aTTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGc                                            >  2:36588/1‑151 (MQ=255)
                                                                                                    tGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCa                                           <  1:549585/150‑1 (MQ=255)
                                                                                                       gTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCt                                       >  1:945858/1‑151 (MQ=255)
                                                                                                       gTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCt                                       >  2:38270/1‑151 (MQ=255)
                                                                                                                 cGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCga                             <  2:798973/151‑1 (MQ=255)
                                                                                                                    gTTTCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCcgcg                              >  2:126057/1‑147 (MQ=255)
                                                                                                                      ttCGGTGACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGAtt                         <  2:266242/150‑1 (MQ=255)
                                                                                                                            gACGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGt                  <  1:914616/151‑1 (MQ=255)
                                                                                                                             aCGATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTa                 <  1:887788/151‑1 (MQ=255)
                                                                                                                               gATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACt               >  2:37315/1‑151 (MQ=255)
                                                                                                                               gATGGAGATAGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACt               <  1:549680/151‑1 (MQ=255)
                                                                                                                                        aGCGCGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATc       >  1:1053277/1‑150 (MQ=255)
                                                                                                                                         gcgcGGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAg     <  1:43238/151‑1 (MQ=255)
                                                                                                                                          cgcgGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAgg    <  2:826930/151‑1 (MQ=255)
                                                                                                                                          cgcgGCGG‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAgg    <  1:175332/151‑1 (MQ=255)
                                                                                                                                            cggcgg‑CCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCg  <  2:294043/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCG  >  NC_000913/701231‑701522

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