Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,457,566 C→T A100A (GCG→GCA)  yjgA ← ribosome‑associated UPF0307 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,5660CT98.4% 232.6 / ‑4.1 63A100A (GCG→GCA) yjgAribosome‑associated UPF0307 family protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (32/30);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

ACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTT  >  NC_000913/4457419‑4457713
                                                                                                                                                   |                                                                                                                                                   
aCGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAgtgc                                                                                                                                                   <  1:555286/150‑1 (MQ=255)
    aGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGAc                                                                                                                                              >  1:329852/1‑151 (MQ=255)
    aGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGAc                                                                                                                                              >  1:671851/1‑151 (MQ=255)
        gCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGaa                                                                                                                                           <  1:853010/150‑1 (MQ=255)
        gCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAAt                                                                                                                                          <  1:674459/151‑1 (MQ=255)
           gACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATa                                                                                                                                         <  1:784465/149‑1 (MQ=255)
            aCGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGc                                                                                                                                      <  1:164027/151‑1 (MQ=255)
            aCGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGc                                                                                                                                      <  1:91494/151‑1 (MQ=255)
             cGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGct                                                                                                                                     >  1:628807/1‑151 (MQ=255)
              gATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCTTCATCACCTTGATCGCTCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGTATACGc                                                                                                                                      >  1:370011/1‑149 (MQ=255)
              gATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGctc                                                                                                                                    <  2:128724/151‑1 (MQ=255)
               aTCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTTCCTGACGAATAGGctct                                                                                                                                   >  2:486941/1‑151 (MQ=255)
                  ggCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACg                                                                                                                                <  2:887471/151‑1 (MQ=255)
                  ggCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACg                                                                                                                                <  2:59811/151‑1 (MQ=255)
                    cATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTc                                                                                                                              >  2:334822/1‑151 (MQ=255)
                        gggCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCgcgc                                                                                                                          >  1:531723/1‑151 (MQ=255)
                            cacaGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGgc                                                                                                                      <  1:626321/151‑1 (MQ=255)
                            cacaGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGgc                                                                                                                      >  2:232775/1‑151 (MQ=255)
                               aGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGgcgc                                                                                                                    >  2:139308/1‑150 (MQ=255)
                                  tttAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCgcagc                                                                                                                 >  1:49945/1‑150 (MQ=255)
                                   ttAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCAt                                                                                                               >  1:763316/1‑151 (MQ=255)
                                    tAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCAt                                                                                                               <  2:269019/150‑1 (MQ=255)
                                       tACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCtt                                                                                                            >  2:799330/1‑150 (MQ=255)
                                             aGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATg                                                                                                     >  2:692026/1‑151 (MQ=255)
                                                gATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGc                                                                                                  >  1:136544/1‑151 (MQ=255)
                                                gATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGc                                                                                                  >  2:5359/1‑151 (MQ=255)
                                                   ggCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGc                                                                                               >  1:76060/1‑151 (MQ=255)
                                                   ggCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTATACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGc                                                                                               <  1:618890/151‑1 (MQ=255)
                                                     cGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCaa                                                                                             <  1:136513/151‑1 (MQ=255)
                                                      gTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAAc                                                                                            <  1:727345/151‑1 (MQ=255)
                                                         atcaCCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCCTCTTACCGATGAGCTGCAACTgg                                                                                         <  1:386639/151‑1 (MQ=255)
                                                            aCCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTggcgg                                                                                      >  2:797723/1‑151 (MQ=255)
                                                                tGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTggcggcgg                                                                                   <  1:741314/150‑1 (MQ=255)
                                                                 gatcgatcAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTggcggcggc                                                                                  <  2:857489/150‑1 (MQ=255)
                                                                 gatcgatcAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGgcggcggcg                                                                                 >  1:504629/1‑151 (MQ=255)
                                                                        cAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTcc                                                                          >  1:89122/1‑151 (MQ=255)
                                                                            cGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCt                                                                      >  1:702915/1‑151 (MQ=255)
                                                                            cGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCt                                                                      >  2:1720/1‑151 (MQ=255)
                                                                             gATCGCGCAGGATTTCCAGGTTATGGAACAGCACCACCTGCTGGTTGCGACGGGTTTTCAGCTTATCCAGTGCCTGACGAATAGCCTCTACGTCGCGCTGGCGCACCATCTTACCGACGAGCTGCAACTGGCGGCGGCGACCTTCCATc                                                                       >  2:630666/1‑149 (MQ=255)
                                                                                cgcgCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTaa                                                                   <  2:201973/150‑1 (MQ=255)
                                                                                  cgcAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATAc                                                                <  2:173070/151‑1 (MQ=255)
                                                                                  cgcAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATAc                                                                >  2:41770/1‑151 (MQ=255)
                                                                                    cAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGc                                                              >  2:696989/1‑151 (MQ=255)
                                                                                         tttCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTggg                                                          <  1:329118/150‑1 (MQ=255)
                                                                                          ttCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCg                                                        >  1:684005/1‑151 (MQ=255)
                                                                                              aGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCt                                                    <  2:531026/151‑1 (MQ=255)
                                                                                                   aTGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAAt                                                >  1:828236/1‑150 (MQ=255)
                                                                                                      gAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGcc                                            <  2:490109/151‑1 (MQ=255)
                                                                                                       aaCAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCcg                                           >  1:712280/1‑151 (MQ=255)
                                                                                                          aGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCcgcgc                                        >  1:538251/1‑151 (MQ=255)
                                                                                                                     ctgGTTGTGACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGc                              >  2:443945/1‑150 (MQ=255)
                                                                                                                          tgtgACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTa                         <  1:795824/150‑1 (MQ=255)
                                                                                                                          tgtgACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTa                         <  1:692075/150‑1 (MQ=255)
                                                                                                                           gtgACGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAAc                       >  2:16960/1‑151 (MQ=255)
                                                                                                                              aCGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACggg                    >  1:830989/1‑151 (MQ=255)
                                                                                                                               cGGTTTTTCAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGa                   <  2:400193/151‑1 (MQ=255)
                                                                                                                                 gTTTTTCAGCTTATCCATTTCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATATCCGCGCGCAGATCCGCATCTAACTGGAt                  <  1:471321/150‑1 (MQ=255)
                                                                                                                                       cAGCTTATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATcc           <  2:750463/151‑1 (MQ=255)
                                                                                                                                           ttATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAgcg       >  2:400415/1‑151 (MQ=255)
                                                                                                                                           ttATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGCGGACATTCCAGCTTGATACGCTGGGCGAGCTCAATATCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAgcg       >  1:616768/1‑151 (MQ=255)
                                                                                                                                            tATCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAgc        <  2:763266/149‑1 (MQ=255)
                                                                                                                                             aTCCAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGcgcg     <  2:477082/151‑1 (MQ=255)
                                                                                                                                               ccAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGtt   <  2:322412/151‑1 (MQ=255)
                                                                                                                                                cAGTGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGttt  <  2:628888/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCTTCCATCTTAATACGCTGGGCGAGCTCAATAGCCGCGCGCAGATCCGCATCTAACGGGATCTTATCCAGCGCGTTT  >  NC_000913/4457419‑4457713

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