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.4% 400.0 / ‑0.9 124coding (65/1386 nt)yifKputative amino acid/amine transport protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (51/71);  total (51/73)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 160 total aligned reads displayed.

AAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCAT  >  NC_000913/3978472‑3978651
                                                                                              |                                                                                          
 aaGCGTGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑g                                                                                          >  2:4543692/1‑93 (MQ=255)
 aaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑g                                                                                          <  1:3025301/93‑1 (MQ=255)
 aaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑g                                                                                          <  2:4185198/93‑1 (MQ=255)
  aGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gg                                                                                         <  2:2280411/93‑1 (MQ=255)
  aGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gg                                                                                         <  1:684845/93‑1 (MQ=255)
   gCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggg                                                                                        <  2:3327818/93‑1 (MQ=255)
   gCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggg                                                                                        <  2:436175/93‑1 (MQ=255)
   gCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggg                                                                                        >  1:970553/1‑93 (MQ=255)
   gCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggg                                                                                        <  1:1653893/93‑1 (MQ=255)
     aTGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggggg                                                                                    >  1:1971331/1‑93 (MQ=255)
      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggggg                                                                                   <  1:3996537/93‑1 (MQ=255)
      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggggg                                                                                   >  1:331716/1‑93 (MQ=255)
       gCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggggggg                                                                                  >  2:1985476/1‑92 (MQ=255)
       gCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggggggg                                                                                  >  2:216266/1‑92 (MQ=255)
       gCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggggggg                                                                                  >  2:4569979/1‑92 (MQ=255)
         aaaCACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggggggca                                                                                >  2:2112199/1‑90 (MQ=255)
         aaaCACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggggggca                                                                                <  2:3534009/93‑4 (MQ=255)
          aacacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggggggcac                                                                               >  2:4245032/1‑89 (MQ=255)
          aacacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggggggcac                                                                               >  1:4244990/1‑89 (MQ=255)
            cacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                                                              <  1:1965211/93‑1 (MQ=255)
              ccacacGAGGGTTAATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCAAACTCATCGCCCTGGGGGGGCACCAtt                                                                            <  2:1154544/91‑1 (MQ=255)
              caacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                                                                            <  1:1064135/93‑1 (MQ=255)
                acacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                                                                          <  1:2712430/93‑1 (MQ=255)
                acacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                                                                          <  2:1932883/93‑1 (MQ=255)
                    gAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                                                                      >  2:1759980/1‑93 (MQ=255)
                    gAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                                                                      <  1:3571023/93‑1 (MQ=255)
                    gAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                                                                      <  1:2306240/93‑1 (MQ=255)
                     aGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                                                                     >  1:3272005/1‑93 (MQ=255)
                      ggATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                                                                    <  1:2470962/93‑1 (MQ=255)
                       gATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                                                                   >  2:4165052/1‑93 (MQ=255)
                         tttATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                                                                 >  1:535408/1‑93 (MQ=255)
                          ttATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTg                                                                <  2:213183/93‑1 (MQ=255)
                           tATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                                                               <  2:3277305/93‑1 (MQ=255)
                             tGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                                                             >  2:2273255/1‑93 (MQ=255)
                             tGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                                                             <  2:4354289/93‑1 (MQ=255)
                              ggCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                                                            <  2:3476747/93‑1 (MQ=255)
                              ggCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                                                            <  2:4262503/93‑1 (MQ=255)
                              ggCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                                                            <  1:1147137/93‑1 (MQ=255)
                               gCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGGACCATTGGCGTCGGCCTGTTTAt                                                           >  2:3377866/1‑93 (MQ=255)
                               gCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                                                           <  1:1607419/93‑1 (MQ=255)
                                cAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                                                           <  2:1311776/93‑1 (MQ=255)
                                cAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg                                                          >  1:1483753/1‑93 (MQ=255)
                                cAGACAACAACCCAGAGCTGCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAtc                                                          <  1:3153795/93‑2 (MQ=255)
                                 aGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg                                                         <  2:3363165/93‑1 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                        <  2:2594527/93‑1 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                        >  2:4684646/1‑93 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCATGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                        >  1:2974556/1‑93 (MQ=255)
                                   aTAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                       <  1:2895800/93‑1 (MQ=255)
                                     aaCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc                                                     <  1:731772/93‑1 (MQ=255)
                                     aaCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc                                                     >  1:156838/1‑93 (MQ=255)
                                         aaCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc                                                 >  1:1219566/1‑93 (MQ=255)
                                         aaCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc                                                 >  1:2380559/1‑93 (MQ=255)
                                         aaCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc                                                 <  1:3351936/93‑1 (MQ=255)
                                          aCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa                                                <  2:2184394/93‑1 (MQ=255)
                                           ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg                                               >  2:3525766/1‑93 (MQ=255)
                                           ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg                                               <  2:1189972/93‑1 (MQ=255)
                                            cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt                                              >  2:39165/1‑93 (MQ=255)
                                             agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGtc                                             >  1:3467738/1‑92 (MQ=255)
                                               agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAcc                                           <  2:609493/93‑1 (MQ=255)
                                                gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAccc                                          >  2:490638/1‑93 (MQ=255)
                                                gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAccc                                          <  2:2803486/93‑1 (MQ=255)
                                                gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAccc                                          <  1:2439267/93‑1 (MQ=255)
                                                 cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCt                                         >  1:2786513/1‑93 (MQ=255)
                                                  tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTg                                        >  1:2606713/1‑93 (MQ=255)
                                                     aGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGaaa                                     <  1:1802302/93‑1 (MQ=255)
                                                     aGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGaaa                                     <  2:813824/93‑1 (MQ=255)
                                                      gCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTTTGGGGGCCGCCAGTACCCTGCAAt                                    >  2:3692110/1‑93 (MQ=255)
                                                      gCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAAt                                    >  1:541099/1‑93 (MQ=255)
                                                          gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggc                                <  1:537381/93‑1 (MQ=255)
                                                          gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggc                                <  2:3712826/93‑1 (MQ=255)
                                                            gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccg                              <  1:2080679/93‑1 (MQ=255)
                                                            gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccg                              <  1:1056724/93‑1 (MQ=255)
                                                             cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgg                             <  1:3562497/93‑1 (MQ=255)
                                                                gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGGACCATTGGCGTCGGCCTGTTTATGGGGGCAGCCAGTACCCTGAAATgggcagggcc                          >  1:895353/1‑93 (MQ=255)
                                                                gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggcc                          <  1:3459047/93‑1 (MQ=255)
                                                                    cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATcc                      <  2:2767404/93‑1 (MQ=255)
                                                                     tCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCg                     <  1:675714/93‑1 (MQ=255)
                                                                       gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTa                   >  1:3598984/1‑93 (MQ=255)
                                                                        aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAt                  <  1:861608/93‑1 (MQ=255)
                                                                        aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAt                  <  1:4457476/93‑1 (MQ=255)
                                                                         cATATCGATCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAtt                 <  1:4642254/93‑1 (MQ=255)
                                                                         cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAtt                 <  2:4550939/93‑1 (MQ=255)
                                                                          atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGGCATCCGTAttg                >  1:4148077/1‑93 (MQ=255)
                                                                          atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttg                <  1:133168/93‑1 (MQ=255)
                                                                          atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttg                >  2:1285968/1‑93 (MQ=255)
                                                                           tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt               >  1:3833681/1‑93 (MQ=255)
                                                                           tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt               >  2:4404369/1‑93 (MQ=255)
                                                                           aatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt               <  1:998689/92‑1 (MQ=255)
                                                                            atCGACCTCATCGCCCTGGGGGGGCACCATAGGCGTCGACCTGTTTATGGGTGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgtt              <  1:3078216/93‑1 (MQ=255)
                                                                            atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgtt              <  1:3975602/93‑1 (MQ=255)
                                                                            atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgtt              <  1:4160700/93‑1 (MQ=255)
                                                                             tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgttg             >  2:2627376/1‑93 (MQ=255)
                                                                                 aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCt         <  1:624180/93‑1 (MQ=255)
                                                                                  cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCta        >  2:1253490/1‑93 (MQ=255)
                                                                                   tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCtat       >  2:2405224/1‑93 (MQ=255)
                                                                                   tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCtat       >  2:3223517/1‑93 (MQ=255)
                                                                                    cATCGCCCGG‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCtata      <  1:3606334/93‑1 (MQ=255)
                                                                                    cATCGCCCT‑‑GGGGAGGCACCATTGGCGTCGGCCTGTTTATGGGGCCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCtata      <  2:1925594/93‑1 (MQ=255)
                                                                                      tCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatc    <  1:757348/93‑1 (MQ=255)
                                                                                       cGCCCT‑‑GGGTGGGCCCCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatca   <  2:869415/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 < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: