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.T98.6% 518.4 / ‑0.8 146coding (536/864 nt)ykfAputative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base T (76/68);  total (76/70)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

TCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAA  >  NC_000913/267557‑267736
                                                                                             |                                                                                          
tcCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGa                                                                                             <  1:236383/93‑1 (MQ=255)
 ccGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAt                                                                                           <  2:3272694/93‑1 (MQ=255)
  cGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAtt                                                                                          <  1:4198453/93‑1 (MQ=255)
  cGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAtt                                                                                          <  1:1834845/93‑1 (MQ=255)
   gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAttt                                                                                        >  1:3530842/1‑93 (MQ=255)
   gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAttt                                                                                        >  2:683399/1‑93 (MQ=255)
    cTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAtttt                                                                                      <  2:4299686/93‑1 (MQ=255)
     tCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAttttt                                                                                     <  1:571525/93‑2 (MQ=255)
     tCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAttttt                                                                                     <  2:2976826/93‑2 (MQ=255)
     tCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAttttt                                                                                     >  1:1049872/1‑92 (MQ=255)
       ggCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAtttttgc                                                                                   <  2:2898964/93‑4 (MQ=255)
       ggCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAtttttgc                                                                                   <  1:4022985/93‑4 (MQ=255)
       ggCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGAtttttgc                                                                                   >  2:3335029/1‑90 (MQ=255)
            tCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTg                                                                              >  2:1967729/1‑93 (MQ=255)
             cACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGa                                                                             <  1:1426270/93‑1 (MQ=255)
               ccTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGata                                                                           >  2:3014211/1‑93 (MQ=255)
               ccTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGata                                                                           <  1:3895657/93‑1 (MQ=255)
                cTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGatat                                                                          <  2:3051090/93‑1 (MQ=255)
                cTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGatat                                                                          <  1:395921/93‑1 (MQ=255)
                cTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGatat                                                                          <  2:2574479/93‑1 (MQ=255)
                cTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGatat                                                                          <  2:2395657/93‑1 (MQ=255)
                 ttAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAtt                                                                         >  2:3594657/1‑93 (MQ=255)
                 ttAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAtt                                                                         >  1:4162996/1‑93 (MQ=255)
                 ttAAGCACCGATTGTAAACGCAGCAAAACAGCACACACACGCTCGACGGGGTAGAACAGCTCATGCAGAAGCAAGATTGTTTCGTCTGAGAtt                                                                         >  2:1791529/1‑93 (MQ=255)
                  tAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttc                                                                        <  2:3830524/93‑1 (MQ=255)
                  tAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttc                                                                        <  1:375366/93‑1 (MQ=255)
                  tAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGCCTGATAtta                                                                        >  1:3462845/1‑92 (MQ=18)
                   aaGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttct                                                                       <  1:856390/93‑1 (MQ=255)
                     gCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgc                                                                     <  1:4639018/93‑1 (MQ=255)
                      cccccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgct                                                                    <  2:3991978/93‑1 (MQ=255)
                       ccccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgctt                                                                   >  1:2910594/1‑93 (MQ=255)
                        cccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttc                                                                  <  2:14931/93‑1 (MQ=255)
                        cccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttc                                                                  <  2:55582/93‑1 (MQ=255)
                        cccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttc                                                                  <  2:851701/93‑1 (MQ=255)
                         ccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttttgcttct                                                                 >  2:4452522/1‑93 (MQ=255)
                         ccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttct                                                                 >  1:1423044/1‑93 (MQ=255)
                          cATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttctg                                                                >  1:2909976/1‑93 (MQ=255)
                          cATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttctg                                                                >  2:346177/1‑93 (MQ=255)
                           aTTGTAGACGGACCGACACATCACCCACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATAttctgcttctgc                                                               <  2:590106/93‑1 (MQ=18)
                            ttGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTgcg                                                              >  2:1341898/1‑93 (MQ=255)
                            ttGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTgcg                                                              >  1:3701948/1‑93 (MQ=255)
                             tGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTgcgc                                                             >  2:2774117/1‑93 (MQ=255)
                              gTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGcc                                                            >  2:4513627/1‑93 (MQ=255)
                              gTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGcc                                                            >  2:2871908/1‑93 (MQ=255)
                                aGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGt                                                          <  2:750329/93‑1 (MQ=255)
                                aGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGAAGATTTTTGCGGCTGATATTCTGCTTTTGCGCCgg                                                          >  2:3125300/1‑92 (MQ=18)
                                 gACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTg                                                         >  1:1469193/1‑93 (MQ=255)
                                 gACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTg                                                         <  1:1762345/93‑1 (MQ=255)
                                 gACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTg                                                         <  1:1978927/93‑1 (MQ=255)
                                   cGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGa                                                       <  1:947789/93‑1 (MQ=255)
                                   cGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGa                                                       >  2:2068726/1‑93 (MQ=255)
                                     gACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACa                                                     <  2:3203952/93‑1 (MQ=255)
                                     gACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACa                                                     >  2:1276089/1‑93 (MQ=255)
                                      aCCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACaa                                                    <  1:1217701/93‑1 (MQ=255)
                                      aCCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACaa                                                    <  1:3986782/93‑1 (MQ=255)
                                         gACACAGCACACACCGGATGGACGGGCTGGGACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAAGTg                                                 >  1:3970253/1‑93 (MQ=18)
                                         gACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTg                                                 <  2:1820910/93‑1 (MQ=255)
                                         gACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTg                                                 >  1:18544/1‑93 (MQ=255)
                                          acacAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGt                                                >  2:2221772/1‑93 (MQ=255)
                                           cacaGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTc                                               >  2:2108022/1‑93 (MQ=255)
                                           cacaGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTc                                               >  2:3974433/1‑93 (MQ=255)
                                            acaGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTcc                                              >  1:417040/1‑93 (MQ=255)
                                               gCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTccccc                                           >  1:3808889/1‑93 (MQ=255)
                                               gCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTccccc                                           >  2:2656695/1‑93 (MQ=255)
                                               gCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGAAATTCTGCTTCTGCGCCGTGGACAACTGTccccc                                           >  2:1246440/1‑93 (MQ=255)
                                                   acacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCg                                       >  2:2481043/1‑93 (MQ=255)
                                                    cacCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGc                                      >  1:4536108/1‑93 (MQ=255)
                                                       cGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTgg                                   >  2:1929374/1‑93 (MQ=255)
                                                        ggATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGt                                  >  1:1855439/1‑93 (MQ=255)
                                                        ggATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGt                                  >  1:831853/1‑93 (MQ=255)
                                                         gATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTa                                 >  2:2699433/1‑93 (MQ=255)
                                                         gATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTa                                 <  2:3402661/93‑1 (MQ=255)
                                                         gATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTa                                 <  2:1783011/93‑1 (MQ=255)
                                                         gATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTa                                 <  1:3249581/93‑1 (MQ=255)
                                                         gATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTa                                 >  1:30596/1‑93 (MQ=255)
                                                           tGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAgg                               >  1:1157566/1‑93 (MQ=255)
                                                            ggACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGt                              >  2:1176083/1‑93 (MQ=255)
                                                               cGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCg                           <  2:2313413/93‑1 (MQ=255)
                                                                   cTGGAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCt                       <  1:2634571/93‑1 (MQ=255)
                                                                     ggAACAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCCCCGTCGACAACTGTCCCCCACCGCTGTTAGGTTCGGCCTTa                     <  2:1364324/93‑1 (MQ=37)
                                                                       aaCAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATc                   <  2:2752761/93‑1 (MQ=255)
                                                                         cAGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCg                 <  2:4105243/93‑1 (MQ=255)
                                                                          aGCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGa                <  2:4004906/93‑1 (MQ=255)
                                                                           gCTCATTCAGCAGCCAGATTTGTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGAc               <  2:4400637/93‑1 (MQ=255)
                                                                           gCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCTCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGAc               <  1:155526/93‑1 (MQ=255)
                                                                           gCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGAc               >  2:1576911/1‑93 (MQ=255)
                                                                           gCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGAc               <  1:891436/93‑1 (MQ=255)
                                                                           gCTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGAc               <  1:1909471/93‑1 (MQ=255)
                                                                            cTCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACt              <  1:4175622/93‑1 (MQ=255)
                                                                             tCATGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTg             <  1:2240684/93‑1 (MQ=255)
                                                                               aTGCAGCAGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGc           >  2:2581882/1‑93 (MQ=255)
                                                                                 gcagcaGGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTg         <  1:1505091/93‑1 (MQ=255)
                                                                                  cagcagGCAGATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGa        <  1:2456417/93‑1 (MQ=255)
                                                                                  cagcagGCAGATTTTTGCGCCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGa        <  2:3259886/93‑1 (MQ=255)
                                                                                     caggcagATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaa     <  2:64654/93‑1 (MQ=255)
                                                                                     caggcagATTTTGGCGTCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaa     <  2:3857482/93‑1 (MQ=255)
                                                                                      aggcagATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaac    >  1:4529306/1‑93 (MQ=255)
                                                                                       ggcagATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaaca   <  1:4447805/93‑1 (MQ=255)
                                                                                        gcagATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaacaa  >  2:110154/1‑93 (MQ=255)
                                                                                        gcagATTTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATaacaa  <  1:130011/93‑1 (MQ=255)
                                                                                             |                                                                                          
TCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTCGGCCTTATCCGACTGGCTGATAACAA  >  NC_000913/267557‑267736

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