Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 267,653 +T coding (532/864 nt) ykfA ← putative GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913267,6491.T97.4% 824.0 / 10.5 233coding (536/864 nt)ykfAputative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (5/1);  new base T (106/121);  total (111/122)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 266 total aligned reads displayed.

TCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAA  >  NC_000913/267557‑267736
                                                                                              |                                                                                       
 ccGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAt                                                                                        >  2:4439731/1‑93 (MQ=255)
  cGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAtt                                                                                       <  2:7506630/93‑1 (MQ=255)
   gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAttt                                                                                      <  1:6424606/93‑1 (MQ=255)
    cTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAtttt                                                                                     <  2:3422859/93‑1 (MQ=255)
     tCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAttttt                                                                                    <  1:1565967/93‑2 (MQ=255)
     tCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAttttt                                                                                    <  2:7213807/93‑2 (MQ=255)
     tCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAttttt                                                                                    >  2:1087080/1‑92 (MQ=255)
     tCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAttttt                                                                                    <  1:4437725/93‑2 (MQ=255)
     tCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAttttt                                                                                    <  1:3790672/93‑2 (MQ=255)
     tCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAttttt                                                                                    <  2:5752022/93‑2 (MQ=255)
      cGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTCGAACAGCTCATGGAGCAGGCAGAtttttg                                                                                   >  1:1130751/1‑91 (MQ=255)
        gCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAtttttgcg                                                                                 >  2:6540006/1‑89 (MQ=255)
        gCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAtttttgcg                                                                                 >  2:4071833/1‑89 (MQ=255)
           aTCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCt                                                                               >  1:2800769/1‑93 (MQ=255)
            tCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTg                                                                              <  1:5983017/93‑1 (MQ=255)
            tCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTg                                                                              <  2:6375346/93‑1 (MQ=255)
             cACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGa                                                                             <  1:4200654/93‑1 (MQ=255)
              aCCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGat                                                                            >  1:3130696/1‑93 (MQ=255)
              aCCTTAAGCCCCCATTGGAAACGGACCGACAAAGCACACCCCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGat                                                                            <  1:2148723/93‑1 (MQ=18)
               ccTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGata                                                                           <  2:1201324/93‑1 (MQ=255)
               ccTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGata                                                                           <  2:1379615/93‑1 (MQ=255)
               ccTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGata                                                                           <  2:3337142/93‑1 (MQ=255)
                 ttAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAtt                                                                         <  1:2983919/93‑1 (MQ=255)
                 ttAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAtt                                                                         >  1:2895568/1‑93 (MQ=255)
                  tAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttc                                                                        <  2:5487819/93‑1 (MQ=255)
                  tAAGCCCCCATTATAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttc                                                                        <  1:68682/93‑1 (MQ=255)
                   aaGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttct                                                                       <  2:3401646/93‑1 (MQ=255)
                   aaGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttct                                                                       <  2:3640718/93‑1 (MQ=255)
                   aaGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttct                                                                       <  2:439875/93‑1 (MQ=255)
                     gCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgc                                                                     <  1:2122161/93‑1 (MQ=255)
                     gCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgc                                                                     >  1:5093318/1‑93 (MQ=255)
                      cccccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGGAGGCAGATTTTTGCGGCTGATAttctgct                                                                    >  1:7323011/1‑93 (MQ=255)
                      cccccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgct                                                                    >  1:3653800/1‑93 (MQ=255)
                        cccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttc                                                                  <  2:6952852/93‑1 (MQ=255)
                         ccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttct                                                                 <  2:2261993/93‑1 (MQ=255)
                         ccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttct                                                                 <  1:1349973/93‑1 (MQ=255)
                          cATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttctg                                                                <  1:895618/93‑1 (MQ=255)
                           aTTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttctgc                                                               <  2:7591763/93‑1 (MQ=255)
                           aTTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttctgc                                                               <  1:4736851/93‑1 (MQ=255)
                           aTTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttctgc                                                               >  2:4024230/1‑93 (MQ=255)
                            ttGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTgcg                                                              >  2:7668998/1‑93 (MQ=255)
                            ttGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTgcg                                                              >  2:4625056/1‑93 (MQ=255)
                             tGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTgcgc                                                             >  2:3228163/1‑93 (MQ=255)
                              gTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGCTTTTTGCGGCTGATATTCTGCTTCTgcgcg                                                            >  1:2180944/1‑92 (MQ=18)
                                 gACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTg                                                         >  1:5621741/1‑93 (MQ=255)
                                 gACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTg                                                         >  2:4682191/1‑93 (MQ=255)
                                  aCGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTgg                                                        <  1:2590179/93‑1 (MQ=255)
                                    ggACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGAc                                                      <  1:6061496/93‑1 (MQ=255)
                                     gACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACa                                                     <  1:4343569/93‑1 (MQ=255)
                                      aCCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGTAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACaa                                                    <  2:4516450/93‑1 (MQ=255)
                                       ccGACACAGCACACGCCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAAc                                                   <  2:816290/93‑1 (MQ=255)
                                       ccGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAAc                                                   <  1:1543748/93‑1 (MQ=255)
                                        cGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGCGGACAACt                                                  >  2:2003096/1‑93 (MQ=255)
                                         gACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTg                                                 >  1:6452397/1‑93 (MQ=255)
                                         gACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCACCGTGGCCAACTg                                                 >  2:3470885/1‑93 (MQ=18)
                                          acacAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGt                                                >  2:4894873/1‑93 (MQ=255)
                                          acacAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGt                                                <  2:1564187/93‑1 (MQ=255)
                                            acaGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTcc                                              <  1:541639/93‑1 (MQ=255)
                                            acaGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTcc                                              >  2:7294658/1‑93 (MQ=255)
                                              aGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTcccc                                            <  1:2715781/93‑1 (MQ=255)
                                              aGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTcccc                                            <  1:2407557/93‑1 (MQ=255)
                                               gCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTccccc                                           <  1:7156272/93‑1 (MQ=255)
                                                  cacacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCAcc                                        >  2:2229568/1‑93 (MQ=255)
                                                   ccacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCg                                       <  2:1703857/92‑1 (MQ=255)
                                                   acacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCg                                       <  2:5535350/93‑1 (MQ=255)
                                                    cacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGc                                      <  2:1000141/93‑1 (MQ=255)
                                                    cacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGc                                      >  2:4294614/1‑93 (MQ=255)
                                                    cacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGc                                      <  2:3766736/93‑1 (MQ=255)
                                                    cacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGc                                      >  2:5179232/1‑93 (MQ=255)
                                                      ccGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTg                                    >  1:2551337/1‑93 (MQ=255)
                                                      ccGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTg                                    <  2:2875752/93‑1 (MQ=255)
                                                      ccGAATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGAACAACTGTCCCCCACCGCTg                                    <  2:295598/93‑1 (MQ=25)
                                                         gATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTa                                 <  2:6608597/93‑1 (MQ=255)
                                                         gATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTa                                 <  2:2405749/93‑1 (MQ=255)
                                                          aTGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAg                                <  2:3100522/93‑1 (MQ=255)
                                                             gACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGtt                             <  2:5990118/93‑1 (MQ=255)
                                                               cgcgCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCg                           <  1:3229668/90‑1 (MQ=255)
                                                               cGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCg                           >  2:416082/1‑93 (MQ=255)
                                                                gggCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCgg                          >  2:4185917/1‑93 (MQ=255)
                                                                 ggCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGc                         <  1:2360088/93‑1 (MQ=255)
                                                                  gCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGcc                        >  1:4610983/1‑93 (MQ=255)
                                                                   cTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCt                       <  1:7644068/93‑1 (MQ=255)
                                                                     ggAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTa                     <  1:7743824/93‑1 (MQ=255)
                                                                       aaCAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATc                   <  1:5840711/93‑1 (MQ=255)
                                                                       aaCAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATc                   <  1:5262304/93‑1 (MQ=255)
                                                                       aaCAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGGAGGTTCGGCCTTATc                   >  1:484920/1‑93 (MQ=255)
                                                                        aCAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATcc                  >  1:3263173/1‑93 (MQ=255)
                                                                           gttCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGAc               <  2:3182841/91‑1 (MQ=255)
                                                                           gCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGAc               >  1:3395234/1‑93 (MQ=255)
                                                                           gCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGAc               >  1:6967230/1‑93 (MQ=255)
                                                                            cTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACt              <  2:2728482/93‑1 (MQ=255)
                                                                            cTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACt              >  1:1795946/1‑93 (MQ=255)
                                                                             tCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCTGCCTTATCCGACTg             >  1:40216/1‑93 (MQ=255)
                                                                              cATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTgg            >  1:38585/1‑93 (MQ=255)
                                                                               aTGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCGCCGCTGGTAGGTTCGGCCTTATCCGACTGGc           <  2:5638195/93‑1 (MQ=255)
                                                                               aTGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGCAACTGGTCCCCAACGCCGGGAGGGTCGGGCTTATCCGGCTgggt          >  1:489972/1‑91 (MQ=16)
                                                                                tGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCt          >  1:2344588/1‑93 (MQ=255)
                                                                                 gcagcaGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTg         >  1:6206351/1‑93 (MQ=255)
                                                                                    gcaggcagATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATa      >  1:7456391/1‑93 (MQ=255)
                                                                                       ggcagATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTGGGTTCGGCCTTATCCGACTGGCTGATaaca   >  1:591452/1‑93 (MQ=255)
                                                                                              |                                                                                       
TCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAA  >  NC_000913/267557‑267736

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