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% 575.3 / 1.8 186coding (65/1386 nt)yifKputative amino acid/amine transport protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (89/94);  total (91/95)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

Only 100 of 238 total aligned reads displayed.

AAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATC  >  NC_000913/3978472‑3978652
                                                                                              |                                                                                         
aaaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                                                            <  1:4155885/93‑1 (MQ=255)
aaaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                                                            <  1:2211168/93‑1 (MQ=255)
aaaGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                                                            <  1:2060131/93‑1 (MQ=255)
  aGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgg                                                                                         <  1:1668010/93‑1 (MQ=255)
   gCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggg                                                                                        >  2:289475/1‑93 (MQ=255)
   gCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggg                                                                                        >  2:2647688/1‑93 (MQ=255)
   gCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggg                                                                                        <  1:1645202/93‑1 (MQ=255)
    cATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                                                                      <  2:1588361/93‑1 (MQ=255)
    cATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                                                                      <  2:2340968/93‑1 (MQ=255)
      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                                                                    >  1:1195080/1‑93 (MQ=255)
      tGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                                                                    >  1:4110568/1‑93 (MQ=255)
       gCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggggg                                                                                   <  2:1470367/93‑2 (MQ=255)
        cAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggggc                                                                                  >  2:136746/1‑91 (MQ=255)
         aaaCACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggggca                                                                                 >  1:1698614/1‑90 (MQ=255)
          aacacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggggcac                                                                                >  2:3695779/1‑89 (MQ=255)
          aacacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggggcac                                                                                >  2:1098209/1‑89 (MQ=255)
            cacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                                                               <  2:1732894/93‑1 (MQ=255)
            cacaacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                                                               <  2:3630075/93‑1 (MQ=255)
               aacgcGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                                                            >  2:2566183/1‑93 (MQ=255)
               aacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                                                            <  1:2798597/93‑1 (MQ=255)
               aacacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                                                            >  1:2852416/1‑93 (MQ=255)
                acacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                                                                           <  2:4044503/93‑1 (MQ=255)
                acacGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                                                                           >  1:635171/1‑93 (MQ=255)
                  acGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                                                                         >  2:1258782/1‑93 (MQ=255)
                   cGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                                                                        <  1:484528/93‑1 (MQ=255)
                   cGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                                                                        >  2:4269181/1‑93 (MQ=255)
                      ggATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                                                                     >  2:3953257/1‑93 (MQ=255)
                      ggATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                                                                     >  2:900510/1‑93 (MQ=255)
                      ggATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                                                                     >  1:3156592/1‑93 (MQ=255)
                         tttATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                                                                  >  1:3342122/1‑93 (MQ=255)
                             tGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                                                              >  1:2435014/1‑93 (MQ=255)
                              ggCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                                                             <  2:1902196/93‑1 (MQ=255)
                               gCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                                                            <  1:1447009/93‑1 (MQ=255)
                                cAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg                                                           <  2:440626/93‑1 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                         <  1:848396/93‑1 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                         <  1:2651702/93‑1 (MQ=255)
                                  gATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                                                         <  1:2625986/93‑1 (MQ=255)
                                   aTAACAACCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATAGCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                                                        <  2:2620511/93‑1 (MQ=255)
                                    tAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg                                                       <  2:2777782/93‑1 (MQ=255)
                                    tAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg                                                       >  1:4273528/1‑93 (MQ=255)
                                    aaacaacCCAGAGTTACAGCGGGGCCTGAAACTCCGAAATATCGACCTAACCCCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg                                                       <  2:3874970/92‑1 (MQ=255)
                                     aaCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc                                                      <  2:21054/93‑1 (MQ=255)
                                      aCAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc                                                     >  1:3265773/1‑93 (MQ=255)
                                        aaaCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc                                                   >  1:1487770/1‑93 (MQ=255)
                                         aaCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc                                                  >  2:4558434/1‑93 (MQ=255)
                                         aaCCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc                                                  >  1:686993/1‑93 (MQ=255)
                                            cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt                                               >  1:3620688/1‑93 (MQ=255)
                                             agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa                                              >  2:3144097/1‑93 (MQ=255)
                                               agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAcc                                            >  2:2390102/1‑93 (MQ=255)
                                                 cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCt                                          <  1:1215082/93‑1 (MQ=255)
                                                 cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCt                                          <  2:470504/93‑1 (MQ=255)
                                                 cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCt                                          <  2:1195212/93‑1 (MQ=255)
                                                  tACAGCGTGGGCTGGAAGCTCGACATATTGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGGCGCCAGTACCCTg                                         >  2:4088964/1‑93 (MQ=255)
                                                  tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTg                                         >  1:559268/1‑93 (MQ=255)
                                                     aGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGaaa                                      >  2:3214762/1‑93 (MQ=255)
                                                      gCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAAt                                     >  1:1384956/1‑93 (MQ=255)
                                                      gCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAAt                                     <  2:2668574/93‑1 (MQ=255)
                                                      gCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGAGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAAt                                     >  1:2509744/1‑93 (MQ=255)
                                                       ggTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATg                                    <  2:516081/92‑1 (MQ=255)
                                                          gggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggc                                 >  1:954572/1‑93 (MQ=255)
                                                           ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggcc                                <  1:1058744/93‑1 (MQ=255)
                                                             cTGGAAGCTCGACATATCGAACTCATCGCCCGGGGGGGGCACCATTGGCGTCGGCCTGTTTAGGGGGGCCGCCAGTACCCTGAAATgggccgg                              <  2:2079873/93‑1 (MQ=255)
                                                             cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgg                              <  2:3479104/93‑1 (MQ=255)
                                                             cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgg                              <  2:1043685/93‑1 (MQ=255)
                                                              tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccggg                             >  1:1299303/1‑93 (MQ=255)
                                                               ggAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                            >  1:3612579/1‑93 (MQ=255)
                                                               ggAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                            >  1:1361984/1‑93 (MQ=255)
                                                               ggAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggc                            >  1:1147626/1‑93 (MQ=255)
                                                                gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATgggccgggcc                           >  1:406382/1‑93 (MQ=255)
                                                                 aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCa                          <  2:152466/93‑1 (MQ=255)
                                                                  aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCTCCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCAt                         <  1:2302011/93‑1 (MQ=255)
                                                                  aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCAt                         <  1:873328/93‑1 (MQ=255)
                                                                  aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCAt                         <  2:1551927/93‑1 (MQ=255)
                                                                  aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCAt                         >  1:3042966/1‑93 (MQ=255)
                                                                  aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGCGCCAt                         >  1:3711092/1‑93 (MQ=255)
                                                                    cTCGCCATATCGAACTCATCGCCCTGGGGGGGCATCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATcc                       <  1:4318575/93‑1 (MQ=255)
                                                                    cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATcc                       <  1:962288/93‑1 (MQ=255)
                                                                     tCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCg                      >  1:1832506/1‑93 (MQ=255)
                                                                     tCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCg                      >  1:2169266/1‑93 (MQ=255)
                                                                       gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTa                    >  1:1103495/1‑93 (MQ=255)
                                                                       gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTa                    <  1:2897567/93‑1 (MQ=255)
                                                                         cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAtt                  <  2:3400704/93‑1 (MQ=255)
                                                                          atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttg                 >  2:4440542/1‑93 (MQ=255)
                                                                           tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgt                <  1:4007305/93‑1 (MQ=255)
                                                                           tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGGACCCTGAAATGGGCCGGGCCATCCGTCttgt                >  1:493910/1‑93 (MQ=255)
                                                                            atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgtt               <  2:2688860/93‑1 (MQ=255)
                                                                             tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgttg              <  1:4402440/93‑1 (MQ=255)
                                                                             tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTAttgttg              <  1:1500202/93‑1 (MQ=255)
                                                                              cGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTgg             <  1:3952478/93‑1 (MQ=255)
                                                                               gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGc            >  2:3534165/1‑93 (MQ=255)
                                                                                aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGcc           >  2:2402891/1‑93 (MQ=255)
                                                                                aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGcc           <  1:2694529/93‑1 (MQ=255)
                                                                                  cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCta         <  2:3039758/93‑1 (MQ=255)
                                                                                    cATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCtata       <  1:926465/93‑1 (MQ=255)
                                                                                    cATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGCCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCtata       <  2:167901/93‑1 (MQ=255)
                                                                                      tCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatc     >  1:300436/1‑93 (MQ=255)
                                                                                       cGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGCGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatca    <  2:1024423/93‑1 (MQ=255)
                                                                                        gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatcat   >  2:1146393/1‑93 (MQ=255)
                                                                                        gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatcat   <  2:3890257/93‑1 (MQ=255)
                                                                                        gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATatcat   >  2:724015/1‑93 (MQ=255)
                                                                                              |                                                                                         
AAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATC  >  NC_000913/3978472‑3978652

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: