Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,727,394 T→C noncoding (71/76 nt) gltW ← tRNA‑Sec

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,727,3940TC93.1% 389.4 / 18.8 145noncoding (71/76 nt)gltWtRNA‑Sec
Reads supporting (aligned to +/- strand):  ref base T (8/2);  new base C (84/51);  total (92/53)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

Only 100 of 153 total aligned reads displayed.

TTTTAAAGAGCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCC  >  NC_000913/2727304‑2727448
                                                                                          |                                                      
 tttAAAGAGCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGcccc                                                     <  2:3997601/93‑1 (MQ=255)
  ttAAAGAGCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCcccc                                                    <  1:6129856/93‑1 (MQ=255)
  ttAAAGAGCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCcccc                                                    >  2:1158282/1‑93 (MQ=255)
  ttAAAGAGCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCcccc                                                    >  1:2145713/1‑93 (MQ=255)
  taaaaGAGCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCcccc                                                    <  2:3603158/91‑1 (MQ=255)
   tAAAGAGCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCt                                                   >  1:3372654/1‑93 (MQ=255)
   tAAAGAGCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGAAAGTGGCGCCCCCt                                                   >  1:989273/1‑93 (MQ=255)
    aaaGAGCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTa                                                  <  1:703397/93‑1 (MQ=255)
    aaaGAGCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTa                                                  <  2:3286658/93‑1 (MQ=255)
      agagCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTTACAAACCAGCAAGTGGCGTCCCCTAgg                                                <  1:3522840/93‑1 (MQ=255)
      agagCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAgg                                                <  1:3564544/93‑1 (MQ=255)
      agagCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAgg                                                <  2:4051759/93‑1 (MQ=255)
      agagCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGCCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAgg                                                <  1:3575192/93‑1 (MQ=255)
       gagCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAggg                                               >  1:2038618/1‑93 (MQ=255)
       gagCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAggg                                               >  2:3217404/1‑93 (MQ=255)
       gagCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAggg                                               >  1:2924099/1‑93 (MQ=255)
       gagCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAggg                                               <  1:5094704/93‑1 (MQ=255)
        agCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAgggg                                              <  2:3724506/93‑1 (MQ=255)
         gCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGa                                             <  1:4310163/93‑1 (MQ=255)
         gCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGa                                             <  2:2928067/93‑1 (MQ=255)
         gCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGa                                             >  2:5930951/1‑93 (MQ=255)
         gCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGa                                             >  1:1880866/1‑93 (MQ=255)
          cAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGAt                                            <  2:2208796/93‑1 (MQ=255)
           aaaTATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGAtt                                           >  1:2685565/1‑93 (MQ=255)
           aaaTATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGAtt                                           <  2:1035188/93‑1 (MQ=255)
            aaTATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTc                                          >  1:6139982/1‑93 (MQ=255)
            aaTATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTc                                          >  2:5475481/1‑93 (MQ=255)
            aaTATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCACCCTAGGGGATTc                                          >  1:5606799/1‑93 (MQ=255)
             atatatCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCg                                         >  2:5460622/1‑93 (MQ=255)
              tatatCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGa                                        >  1:3451208/1‑93 (MQ=255)
              tatatCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGa                                        >  1:1122324/1‑93 (MQ=255)
              tatatCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGa                                        >  2:2482057/1‑93 (MQ=255)
               atatCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGaa                                       <  1:3220810/93‑1 (MQ=255)
               atatCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGaa                                       >  1:5881051/1‑93 (MQ=255)
                tatCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAAc                                      >  2:2697376/1‑93 (MQ=255)
                tatCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAAc                                      <  2:1812869/93‑1 (MQ=255)
                tatCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAAc                                      <  2:2150414/93‑1 (MQ=255)
                tatCAAACAAGACTTAACAGTCTGATTTGAGATATTGAGGATGGCGACTTTCACTCACAAATCAGAAAGTGGCGCCCCCTAGGGGCTTCGAAc                                      >  1:3797722/1‑93 (MQ=18)
                 atCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAAcc                                     >  2:5050119/1‑93 (MQ=255)
                   cAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAAcccc                                   >  2:4175207/1‑93 (MQ=255)
                   cAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAAcccc                                   <  2:4073827/93‑1 (MQ=255)
                   cAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAAcccc                                   >  2:5751928/1‑93 (MQ=255)
                      aCAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGt                                >  2:3996076/1‑93 (MQ=255)
                      aCAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGt                                >  2:1293543/1‑93 (MQ=255)
                       cAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGtt                               >  2:5620206/1‑93 (MQ=255)
                       cAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGtt                               >  2:3276891/1‑93 (MQ=255)
                       cAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGtt                               >  2:950034/1‑93 (MQ=255)
                        aaGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCAAACCCCTGttc                              >  1:4466447/1‑92 (MQ=255)
                        aaGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTa                              >  1:1066264/1‑93 (MQ=255)
                        aaGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTa                              >  2:4809497/1‑93 (MQ=255)
                        aaGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTa                              >  1:2399856/1‑93 (MQ=255)
                        aaGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTa                              >  1:461251/1‑93 (MQ=255)
                         aGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTAc                             >  1:179855/1‑93 (MQ=255)
                            cTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTAccgc                          >  1:4725143/1‑93 (MQ=255)
                             ttAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTAccgcc                         >  2:1913836/1‑93 (MQ=255)
                               aaCAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGt                       >  2:5603773/1‑93 (MQ=255)
                                aCAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTg                      >  2:928936/1‑93 (MQ=255)
                                aCAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACACCTGTTACCGCCGTg                      >  2:882239/1‑93 (MQ=255)
                                 cAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGa                     >  1:1068798/1‑93 (MQ=255)
                                 cAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGa                     <  2:450730/93‑1 (MQ=255)
                                  aGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGaa                    >  1:5336904/1‑93 (MQ=255)
                                  aGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGaa                    <  1:1548463/93‑1 (MQ=255)
                                   gTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGaaa                   <  1:5097312/93‑1 (MQ=255)
                                   gTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGaaa                   <  2:2390825/93‑1 (MQ=255)
                                   gTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGaaa                   >  2:767969/1‑93 (MQ=255)
                                   gTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGaaa                   >  1:3099829/1‑93 (MQ=255)
                                   gTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGaaa                   <  2:1334218/93‑1 (MQ=255)
                                    tCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAg                  <  2:1518559/93‑1 (MQ=37)
                                    tCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAg                  >  1:4505584/1‑93 (MQ=37)
                                     cTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAgg                 <  1:2571018/93‑1 (MQ=25)
                                     cTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAgg                 <  1:3347429/93‑1 (MQ=25)
                                      tGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGATCCCCTGTTACCGCCGTGAAAggg                <  1:5694238/93‑1 (MQ=255)
                                      tGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAggg                >  1:3676643/1‑93 (MQ=25)
                                      tGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAggg                >  2:1993605/1‑93 (MQ=25)
                                       gTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGc               <  2:2197427/93‑1 (MQ=18)
                                       gTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGc               >  1:5167324/1‑93 (MQ=21)
                                        ttttGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCg              >  2:4519120/1‑93 (MQ=21)
                                         tttGAGATATTGATGTCTGCGATTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGTTTCGAACTCCTGGTTCTGCTGTGAAAGGGCgg             >  1:5467057/1‑93 (MQ=11)
                                         tttGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCgg             <  2:3383231/93‑1 (MQ=18)
                                         tttGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCgg             >  2:3244135/1‑93 (MQ=21)
                                          ttGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGgt            >  1:1461851/1‑93 (MQ=21)
                                           tGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGgtg           <  2:3006908/93‑1 (MQ=33)
                                           tGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGAGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGgtg           >  2:2275361/1‑93 (MQ=255)
                                            gagaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGgtgt          <  2:2299444/93‑1 (MQ=18)
                                             agaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTc         >  1:3246418/1‑93 (MQ=21)
                                              gaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTcc        >  2:1453931/1‑93 (MQ=18)
                                              gaTAGTGAGGTCGGCGACTTTCACTCACAACCCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTcc        <  1:2531560/93‑1 (MQ=16)
                                               atatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCt       <  2:4429626/93‑1 (MQ=18)
                                                tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTg      >  2:3499684/1‑93 (MQ=18)
                                                tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTg      >  1:675868/1‑93 (MQ=18)
                                                tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTg      <  1:5860556/93‑1 (MQ=18)
                                                tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTATCCTg      >  2:4943705/1‑93 (MQ=16)
                                                  ttGAGGTCGGCGACTTTCCCTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTCCCGCCGTGACAGGGCGGTGTCCTggg    >  1:1869908/1‑93 (MQ=11)
                                                  ttGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTggg    >  2:5660157/1‑93 (MQ=18)
                                                  ttGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTggg    >  2:670802/1‑93 (MQ=18)
                                                  ttGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTggg    >  2:730943/1‑93 (MQ=18)
                                                  ttGAGGTCGACGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTggg    >  1:1181267/1‑93 (MQ=16)
                                                   tGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGc   <  1:5742560/93‑1 (MQ=17)
                                                    gAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGcc  >  1:1779879/1‑93 (MQ=14)
                                                    gAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGCCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGcc  <  1:3016131/93‑1 (MQ=14)
                                                                                          |                                                      
TTTTAAAGAGCAAATATATCAAACAAGACTTAACAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCC  >  NC_000913/2727304‑2727448

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