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.G98.2% 708.6 / 4.2 221coding (65/1386 nt)yifKputative amino acid/amine transport protein
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base G (117/100);  total (120/101)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

AAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCAT  >  NC_000913/3978472‑3978651
                                                                                             |                                                                                          
aaaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                                                             <  1:4400294/93‑1 (MQ=255)
   gCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggg                                                                                       <  2:3500054/93‑1 (MQ=255)
     aTGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggg                                                                                    >  1:2107971/1‑93 (MQ=255)
     aTGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggg                                                                                    <  1:3037397/93‑1 (MQ=255)
      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                                                                   >  2:4843818/1‑93 (MQ=255)
      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                                                                   <  1:231967/93‑1 (MQ=255)
        cAAACACAACACGAGGATTTATGGCAGATAACAACCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggggc                                                                                 <  2:6263199/93‑3 (MQ=255)
          aacacaacacGAGGATTTATGGCAGATAACAAACCCGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggggcac                                                                               >  2:4500047/1‑89 (MQ=255)
          aacacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggggcac                                                                               >  2:5236114/1‑89 (MQ=255)
            cacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                                                              <  1:5826030/93‑1 (MQ=255)
             acaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                                                             <  1:2355320/93‑1 (MQ=255)
             acaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                                                             <  2:5561130/93‑1 (MQ=255)
             acaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                                                             <  1:2662547/93‑1 (MQ=255)
             acaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                                                             >  2:5443517/1‑93 (MQ=255)
               aacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                                                           >  2:2324387/1‑93 (MQ=255)
               aacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                                                           >  1:6372445/1‑93 (MQ=255)
               aacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                                                           >  1:6548403/1‑93 (MQ=255)
               aacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATAGCCCTGGGGGGGGACCATTg                                                                           >  1:1375568/1‑93 (MQ=255)
                acacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                                                                          >  1:1707366/1‑93 (MQ=255)
                 cacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                                                                         >  2:4189744/1‑93 (MQ=255)
                  acGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                                                                        >  1:5208160/1‑93 (MQ=255)
                   cGAGGATTTATGCCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCAACGCCCTGGGGGGGCACCATTGGCGt                                                                       <  1:3179057/93‑1 (MQ=255)
                    gAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                                                                      >  1:545828/1‑93 (MQ=255)
                     aGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                                                                     <  2:1021502/93‑1 (MQ=255)
                      ggATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                                                                    >  2:6467978/1‑93 (MQ=255)
                      ggATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                                                                    >  2:1589948/1‑93 (MQ=255)
                      ggATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                                                                    >  1:4320677/1‑93 (MQ=255)
                       gATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                                                                   >  2:3019076/1‑93 (MQ=255)
                         tttATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                                                                 >  2:2968015/1‑93 (MQ=255)
                         tttATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                                                                 <  1:2561890/93‑1 (MQ=255)
                         tttATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                                                                 <  2:2076253/93‑1 (MQ=255)
                             tGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                                                             >  1:1144089/1‑93 (MQ=255)
                              agCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                                                            <  2:3664712/92‑1 (MQ=255)
                               gCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                                                           >  2:4891042/1‑93 (MQ=255)
                                cAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg                                                          <  2:1833425/93‑1 (MQ=255)
                                cAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg                                                          <  1:3251728/93‑1 (MQ=255)
                                cAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg                                                          <  1:154668/93‑1 (MQ=255)
                                 aGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg                                                         >  2:4059960/1‑93 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                        <  2:925448/93‑1 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                        <  1:5915791/93‑1 (MQ=255)
                                   aTAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                       <  1:1882049/93‑1 (MQ=255)
                                   aTAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCCGCCTGTTTATgggg                                                       <  1:4063705/93‑1 (MQ=255)
                                    tAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg                                                      >  2:2010326/1‑93 (MQ=255)
                                     aaCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc                                                     <  1:116579/93‑1 (MQ=255)
                                      aCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc                                                    >  1:4783583/1‑93 (MQ=255)
                                         aaCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc                                                 >  1:1256882/1‑93 (MQ=255)
                                         aaCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc                                                 >  2:2261720/1‑93 (MQ=255)
                                         aaCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgca                                                 >  2:2108157/1‑92 (MQ=255)
                                          aCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                >  2:1147875/1‑93 (MQ=255)
                                          aCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                >  1:5008102/1‑93 (MQ=255)
                                           ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGTGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg                                               >  1:5284048/1‑93 (MQ=255)
                                            cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt                                              >  2:4229211/1‑93 (MQ=255)
                                             agagCTACAGCGTGGGCTGTAAGCTCGACATATCGAACTCATCTCCCTGTGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGGCGCCAGtc                                             >  2:2502771/1‑92 (MQ=255)
                                             agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa                                             >  2:1310682/1‑93 (MQ=255)
                                              gagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAc                                            >  2:3525640/1‑93 (MQ=255)
                                              gagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAc                                            >  1:335465/1‑93 (MQ=255)
                                                gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTTGGCCTGTTTATGGGGGCCGCCAGTAccc                                          >  1:6717951/1‑93 (MQ=255)
                                                gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAccc                                          >  2:148851/1‑93 (MQ=255)
                                                gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAccc                                          <  2:3864930/93‑1 (MQ=255)
                                                gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAccc                                          <  2:3420369/93‑1 (MQ=255)
                                                 cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCt                                         <  1:5756527/93‑1 (MQ=255)
                                                  tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTg                                        <  1:4759346/93‑1 (MQ=255)
                                                   aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGa                                       >  2:272040/1‑93 (MQ=255)
                                                     aGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGaaa                                     <  2:1542696/93‑1 (MQ=255)
                                                      gCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAAt                                    <  1:1833131/93‑1 (MQ=255)
                                                       ggTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATg                                   <  2:1000389/92‑1 (MQ=255)
                                                       cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATg                                   <  2:6397965/93‑1 (MQ=255)
                                                        gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgg                                  <  2:5366112/93‑1 (MQ=255)
                                                        gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgg                                  <  1:3737689/93‑1 (MQ=255)
                                                         tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATggg                                 <  2:4878307/93‑1 (MQ=255)
                                                          gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggc                                <  1:4691793/93‑1 (MQ=255)
                                                             cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgg                             >  2:4993254/1‑93 (MQ=255)
                                                             cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgg                             >  1:5966936/1‑93 (MQ=255)
                                                              tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccggg                            <  2:694864/93‑1 (MQ=255)
                                                                gAAGCTCGACTTAGCGACCTCAGCGCCCTGGGGGGGGACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggcc                          <  1:2620414/93‑1 (MQ=255)
                                                                 aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCa                         >  2:6008513/1‑93 (MQ=255)
                                                                   gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTCCCCTGAAATGGGCCGGGCCATc                       <  2:3372116/93‑1 (MQ=255)
                                                                      cGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGt                    >  1:5707287/1‑93 (MQ=255)
                                                                      cGACATATCGAACTCATCGCCCTGGGTGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGt                    >  2:774259/1‑93 (MQ=255)
                                                                       gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTa                   <  1:4216156/93‑1 (MQ=255)
                                                                        aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAt                  >  1:609561/1‑93 (MQ=255)
                                                                         cATATCGACCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAtt                 <  2:6580528/93‑1 (MQ=255)
                                                                         cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAtt                 <  1:706284/93‑1 (MQ=255)
                                                                          atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGACGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttg                >  2:360232/1‑93 (MQ=255)
                                                                           tatCGAACTCATCGCCCTGGGGGGGCACCGTTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGACGGGCCATCCGTAttgt               <  2:415925/93‑1 (MQ=255)
                                                                           tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt               <  1:2146126/93‑1 (MQ=255)
                                                                            atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgtt              <  2:4752233/93‑1 (MQ=255)
                                                                            atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgtt              >  2:1224008/1‑93 (MQ=255)
                                                                             tCGACCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgttg             <  2:5844021/93‑1 (MQ=255)
                                                                             tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgttg             >  1:87304/1‑93 (MQ=255)
                                                                               gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGc           >  1:3556834/1‑93 (MQ=255)
                                                                                 aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCt         <  2:4574090/93‑1 (MQ=255)
                                                                                 aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATGGTTGGCCt         <  1:6533874/93‑1 (MQ=255)
                                                                                   tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCtat       >  1:3778907/1‑93 (MQ=255)
                                                                                   tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCtat       >  1:4432165/1‑93 (MQ=255)
                                                                                    cATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCtata      <  1:4088961/93‑1 (MQ=255)
                                                                                     aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTTTTGTTGGCCTatat     >  1:5142125/1‑93 (MQ=255)
                                                                                        gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatcat  <  1:2918677/93‑1 (MQ=255)
                                                                                        gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatcat  <  1:4885273/93‑1 (MQ=255)
                                                                                        gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatcat  >  2:1874057/1‑93 (MQ=255)
                                                                                             |                                                                                          
AAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCAT  >  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: