Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,978,570 (G)6→7 coding (71/1386 nt) yifK → putative amino acid/amine transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,978,5641.G97.6% 1184.1 / 16.3 374coding (65/1386 nt)yifKputative amino acid/amine transport protein
Reads supporting (aligned to +/- strand):  ref base . (3/6);  new base G (138/227);  total (141/233)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 455 total aligned reads displayed.

AAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCAT  >  NC_000913/3978472‑3978651
                                                                                             |                                                                                          
aaaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                                                             <  2:4432299/93‑1 (MQ=255)
  aGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gg                                                                                        >  1:3999068/1‑93 (MQ=255)
  aGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gg                                                                                        <  2:2703204/93‑1 (MQ=255)
  aGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gg                                                                                        >  1:2490545/1‑93 (MQ=255)
  aGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gg                                                                                        >  1:8530891/1‑93 (MQ=255)
   gCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGCACTCATCCCCCT‑‑ggg                                                                                       >  1:9810036/1‑93 (MQ=255)
   gCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggg                                                                                       <  2:2575076/93‑1 (MQ=255)
    cATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggg                                                                                     <  2:1506390/93‑1 (MQ=255)
     aTGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggggg                                                                                    <  2:9537448/93‑1 (MQ=255)
     aTGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggggg                                                                                    <  2:7668389/93‑1 (MQ=255)
      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggggg                                                                                   >  2:4643551/1‑93 (MQ=255)
      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggggg                                                                                   >  1:949975/1‑93 (MQ=255)
        cAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggggggc                                                                                 <  1:1756168/93‑3 (MQ=255)
            cacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                                                              <  1:2167111/93‑1 (MQ=255)
             acaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                                                             <  2:4761691/93‑1 (MQ=255)
              caacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                                                                            <  2:6383314/93‑1 (MQ=255)
                acacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                                                                          <  2:3850773/93‑1 (MQ=255)
                 cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                                                                         >  1:9311584/1‑93 (MQ=255)
                 cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                                                                         <  1:2800057/93‑1 (MQ=255)
                  acGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                                                                        >  2:9402393/1‑93 (MQ=255)
                   cGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                                                                       <  2:2021449/93‑1 (MQ=255)
                    gAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                                                                      >  1:10038702/1‑93 (MQ=255)
                         tttATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                                                                 <  2:4304185/93‑1 (MQ=255)
                           tATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                                                               <  1:223919/93‑1 (MQ=255)
                            aTGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                                                              <  2:350169/93‑1 (MQ=255)
                              tgCAGATAACAAACCAGAGGTACAGCGTGGGCTGGAAACTCGACATATCGAACTCATCCCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                                                            <  2:1668754/92‑1 (MQ=255)
                              ggCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                                                            <  1:9156056/93‑1 (MQ=255)
                              ggCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                                                            <  2:152978/93‑1 (MQ=255)
                              ggCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                                                            <  2:3032062/93‑1 (MQ=255)
                              ggCACAAAACAAACCAGAGCTACAGCGGGGGCTGGAAGCTCAACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                                                            <  2:3932702/93‑1 (MQ=255)
                               gCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGGCACCATTGGCGTCGGCCTGTTTa                                                            <  2:4182893/93‑1 (MQ=255)
                                cAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg                                                          >  1:7678747/1‑93 (MQ=255)
                                 aGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg                                                         >  1:6878750/1‑93 (MQ=255)
                                 aGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg                                                         <  2:9006851/93‑1 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                        <  1:6228151/93‑1 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                        >  1:2630355/1‑93 (MQ=255)
                                   aTAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                       <  2:3169948/93‑1 (MQ=255)
                                     aaCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc                                                     <  2:5996025/93‑1 (MQ=255)
                                     aaCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGAGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc                                                     >  1:9683788/1‑93 (MQ=255)
                                      aCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc                                                    >  2:2897747/1‑93 (MQ=255)
                                        aaaCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc                                                  <  1:368889/93‑1 (MQ=255)
                                           ccAGAGTTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg                                               <  2:5206511/93‑1 (MQ=255)
                                            cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt                                              <  1:3498716/93‑1 (MQ=255)
                                             agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa                                             <  1:8414706/93‑1 (MQ=255)
                                             agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa                                             >  1:3079377/1‑93 (MQ=255)
                                             agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa                                             >  2:3258359/1‑93 (MQ=255)
                                             agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa                                             <  1:6692647/93‑1 (MQ=255)
                                              gagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACGATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAc                                            >  2:9138607/1‑93 (MQ=255)
                                                gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAccc                                          <  2:9829441/93‑1 (MQ=255)
                                                gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAccc                                          <  1:4255537/93‑1 (MQ=255)
                                                gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAccc                                          <  1:4441572/93‑1 (MQ=255)
                                                 cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCt                                         <  1:9834193/93‑1 (MQ=255)
                                                 cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCt                                         <  2:10002024/93‑1 (MQ=255)
                                                   aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCACTACCCTGa                                       >  2:7564265/1‑93 (MQ=255)
                                                      gCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAAt                                    >  2:7485810/1‑93 (MQ=255)
                                                      gCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAAt                                    >  1:5850981/1‑93 (MQ=255)
                                                      gCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAAt                                    >  1:7274130/1‑93 (MQ=255)
                                                       cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGTGCCGCCAGTACCCTGAAATg                                   <  1:5994702/93‑1 (MQ=255)
                                                       cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATg                                   <  2:4338762/93‑1 (MQ=255)
                                                       cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATg                                   >  1:1878811/1‑93 (MQ=255)
                                                        gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgg                                  >  1:4543458/1‑93 (MQ=255)
                                                         tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATggg                                 <  1:9552829/93‑1 (MQ=255)
                                                         tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATggg                                 <  1:5498694/93‑1 (MQ=255)
                                                           ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggcc                               <  1:4161539/93‑1 (MQ=255)
                                                            gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccg                              <  1:6378501/93‑1 (MQ=255)
                                                            gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccg                              <  2:8888028/93‑1 (MQ=255)
                                                             cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgg                             <  1:8802138/93‑1 (MQ=255)
                                                             cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgg                             <  1:213114/93‑1 (MQ=255)
                                                             cTGGAAGCTCGACATAACGAACTCATCGCCCT‑‑GAGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgg                             <  1:8985369/93‑1 (MQ=255)
                                                              tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccggg                            <  1:703387/93‑1 (MQ=255)
                                                              tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccggg                            <  1:3607013/93‑1 (MQ=255)
                                                                  aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCAt                        <  1:7340331/93‑1 (MQ=255)
                                                                   gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATc                       <  1:1208740/93‑1 (MQ=255)
                                                                     tCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCg                     <  2:6676335/93‑1 (MQ=255)
                                                                       gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTa                   <  2:1068401/93‑1 (MQ=255)
                                                                        aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAt                  <  1:781339/93‑1 (MQ=255)
                                                                        aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAt                  >  2:8054089/1‑93 (MQ=255)
                                                                        aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAt                  <  1:4795074/93‑1 (MQ=255)
                                                                           tctCGAACTCATCGCCGTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGGACCCTGAAATGGGCGGGGCCATCCGTAttgt               <  2:4920469/91‑1 (MQ=255)
                                                                           tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt               <  1:1285749/93‑1 (MQ=255)
                                                                           tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt               >  2:1816162/1‑93 (MQ=255)
                                                                            atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgtt              <  2:8209017/93‑1 (MQ=255)
                                                                            atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgtt              <  2:2404762/93‑1 (MQ=255)
                                                                             tCGAACTCATCGCCCTGGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgtt              >  1:7995774/1‑93 (MQ=255)
                                                                             tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgttg             >  2:7377828/1‑93 (MQ=255)
                                                                              cGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTgg            <  1:5243540/93‑1 (MQ=255)
                                                                                aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGcc          <  1:3468078/93‑1 (MQ=255)
                                                                                aaCTCATCGCCCTGGGGGGGCACCATTGGCATCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGcc          >  2:5460435/1‑93 (MQ=255)
                                                                                 aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCt         <  2:9303928/93‑1 (MQ=255)
                                                                                  cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCta        <  1:1511525/93‑1 (MQ=255)
                                                                                     ccgccccG‑‑GGGGGGCCCCATTGGCTTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatc    <  2:8608097/89‑1 (MQ=255)
                                                                                     aTCGCCCTGGGGGGGCACCATGGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTatat     <  2:716556/93‑1 (MQ=255)
                                                                                      tCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatc    <  1:1659708/93‑1 (MQ=255)
                                                                                      tCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatc    >  1:1066556/1‑93 (MQ=255)
                                                                                       cGCCCTGGGGGGGCCCCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatca   <  2:3627385/93‑1 (MQ=255)
                                                                                       cGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatca   <  1:3193369/93‑1 (MQ=255)
                                                                                       cGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatca   <  2:213058/93‑1 (MQ=255)
                                                                                        gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatcat  >  2:7808157/1‑93 (MQ=255)
                                                                                        gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatcat  <  2:6085332/93‑1 (MQ=255)
                                                                                        gCCCTGGGGGGGCACCATTGGCGACGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatcat  <  2:6891632/93‑1 (MQ=255)
                                                                                             |                                                                                          
AAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCAT  >  NC_000913/3978472‑3978651

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