Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 488,917 A→G E1053G (GAG→GGG)  aefA → putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913488,9170AG98.8% 527.3 / ‑2.5 165E1053G (GAG→GGG) aefAputative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (86/77);  total (87/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 168 total aligned reads displayed.

TGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACG  >  NC_000913/488827‑489006
                                                                                          |                                                                                         
tGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGc                                                                                         >  2:3634944/1‑93 (MQ=255)
tGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGc                                                                                         >  2:3503244/1‑93 (MQ=255)
 gAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCt                                                                                        >  2:188129/1‑93 (MQ=255)
  aaGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTg                                                                                       >  1:1059081/1‑93 (MQ=255)
   aGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGc                                                                                      >  1:294349/1‑93 (MQ=255)
   aGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGc                                                                                      <  2:4703527/93‑1 (MQ=255)
   aGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGCGCTGc                                                                                      >  1:3441930/1‑93 (MQ=255)
    ggcggcGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCg                                                                                     <  2:637456/93‑1 (MQ=255)
     gcggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGt                                                                                    <  1:419407/93‑1 (MQ=255)
          gACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTa                                                                               >  1:707695/1‑93 (MQ=255)
          gACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTa                                                                               <  1:2304975/93‑1 (MQ=255)
           aCTGACCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAt                                                                              <  1:2638422/93‑1 (MQ=255)
             tGAGCACCCAAGGGTGATGCACGAATCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgt                                                                            >  2:3396129/1‑93 (MQ=255)
             tGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgt                                                                            >  2:986309/1‑93 (MQ=255)
             tGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgt                                                                            >  2:3147749/1‑93 (MQ=255)
              gAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgtg                                                                           >  2:1156010/1‑93 (MQ=255)
                  aCCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTg                                                                       <  2:1799112/93‑1 (MQ=255)
                  aCCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTg                                                                       <  1:1137445/93‑1 (MQ=255)
                    ccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGaa                                                                     >  1:1206992/1‑93 (MQ=255)
                    ccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGaa                                                                     <  1:357537/93‑1 (MQ=255)
                    ccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGaa                                                                     <  1:2902060/93‑1 (MQ=255)
                    ccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGaa                                                                     >  2:1827146/1‑93 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAAc                                                                    >  1:4494432/1‑93 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAAc                                                                    <  2:1400077/93‑1 (MQ=255)
                     cAAGGGTGATACACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAAc                                                                    <  2:820031/93‑1 (MQ=255)
                      aaGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACt                                                                   <  2:1713774/93‑1 (MQ=255)
                      aaGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACt                                                                   >  1:2164026/1‑93 (MQ=255)
                      aaGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACt                                                                   <  2:3920432/93‑1 (MQ=255)
                       aGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTg                                                                  >  1:2259860/1‑93 (MQ=255)
                        gggTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGc                                                                 >  1:1677848/1‑93 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCGTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGt                                                               <  2:609057/93‑1 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGt                                                               <  2:289995/93‑1 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGt                                                               <  1:3045209/93‑1 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGt                                                               >  2:2005754/1‑93 (MQ=255)
                            gATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAATTGCGTGa                                                             >  2:1750774/1‑93 (MQ=255)
                              tGCACTAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGAcc                                                           >  1:4694941/1‑93 (MQ=255)
                              tGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGAcc                                                           <  1:2113660/93‑1 (MQ=255)
                              tGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGAcc                                                           >  2:3040397/1‑93 (MQ=255)
                                 aCGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTa                                                        <  2:134583/93‑1 (MQ=255)
                                 aCGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTa                                                        <  2:1232246/93‑1 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAg                                                       <  2:1894944/93‑1 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAg                                                       <  2:262002/93‑1 (MQ=255)
                                   gAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGt                                                      <  1:1384303/93‑1 (MQ=255)
                                   gAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGt                                                      >  2:2797864/1‑93 (MQ=255)
                                      ccAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGt                                                   <  2:1864485/93‑1 (MQ=255)
                                      ccAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGt                                                   >  2:3928110/1‑93 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTa                                                  <  2:1558622/93‑1 (MQ=255)
                                        aaggCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTAc                                                 <  1:3338629/90‑1 (MQ=255)
                                        aaTGTCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTAc                                                 <  1:600873/93‑1 (MQ=255)
                                           gCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGt                                              >  1:2169249/1‑93 (MQ=255)
                                           gCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGt                                              >  1:1921875/1‑93 (MQ=255)
                                           gCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGt                                              <  1:1192038/93‑1 (MQ=255)
                                             cGGAAGTCTTCTTTACGGCATATGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCg                                            <  2:3062316/93‑1 (MQ=255)
                                              ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                           <  1:2710521/93‑1 (MQ=255)
                                               gAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGAt                                          >  2:2445710/1‑93 (MQ=255)
                                               gAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGAt                                          >  2:4338519/1‑93 (MQ=255)
                                                gaGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATg                                         <  1:3816121/92‑1 (MQ=255)
                                                aaGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATg                                         <  1:333380/93‑1 (MQ=255)
                                                  gTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAg                                       >  2:2961138/1‑93 (MQ=255)
                                                  gTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAg                                       >  1:3994355/1‑93 (MQ=255)
                                                     ttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                    >  2:3019599/1‑93 (MQ=255)
                                                      tcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGa                                   >  1:746754/1‑93 (MQ=255)
                                                      tcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGa                                   <  1:4685008/93‑1 (MQ=255)
                                                      tcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGa                                   >  1:4342295/1‑93 (MQ=255)
                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 >  1:993652/1‑93 (MQ=255)
                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 <  2:3875146/93‑1 (MQ=255)
                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 >  1:2384782/1‑93 (MQ=255)
                                                         ttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAcc                                >  2:2776745/1‑93 (MQ=255)
                                                         tgaCGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAcc                                <  1:378963/91‑1 (MQ=255)
                                                          tACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCg                               <  2:673767/93‑1 (MQ=255)
                                                            cGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTa                             <  2:1155421/93‑1 (MQ=255)
                                                            cGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTa                             >  1:305440/1‑93 (MQ=255)
                                                             ggCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTAc                            <  2:2226320/93‑1 (MQ=255)
                                                               gaTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTa                          <  1:4082060/92‑1 (MQ=255)
                                                               cATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTa                          >  1:3562137/1‑93 (MQ=255)
                                                                aTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTAt                         >  1:1084948/1‑93 (MQ=255)
                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        >  1:1328731/1‑93 (MQ=255)
                                                                   tGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGa                      <  2:3996337/93‑1 (MQ=255)
                                                                    ggTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGAt                     <  1:2854781/93‑1 (MQ=255)
                                                                    ggTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGAt                     <  1:4006184/93‑1 (MQ=255)
                                                                      tGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCa                   <  2:768997/93‑1 (MQ=255)
                                                                      tGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCa                   <  2:3756835/93‑1 (MQ=255)
                                                                      tGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCa                   <  2:2114659/93‑1 (MQ=255)
                                                                          aacACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtg               <  1:2117906/91‑1 (MQ=255)
                                                                          aGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtg               >  1:3683261/1‑93 (MQ=255)
                                                                          aGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtg               >  2:1716762/1‑93 (MQ=255)
                                                                           gCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtgt              >  2:4532965/1‑93 (MQ=255)
                                                                           gCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtgt              >  1:2996727/1‑93 (MQ=255)
                                                                             aCGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGc            >  2:351915/1‑93 (MQ=255)
                                                                               gTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCg          >  2:594768/1‑93 (MQ=255)
                                                                               gTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCg          >  2:372402/1‑93 (MQ=255)
                                                                               gTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCg          >  2:3282235/1‑93 (MQ=255)
                                                                                ttGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGt         >  2:3114339/1‑93 (MQ=255)
                                                                                ttGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGt         <  1:2671479/93‑1 (MQ=255)
                                                                                   gATCATGGGCTGCGTCTGTATGTGCGTGATCTGCGGGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaa      <  1:1285067/93‑1 (MQ=255)
                                                                                    aTCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaa     >  1:3955260/1‑93 (MQ=255)
                                                                                    aTCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaa     >  1:69593/1‑93 (MQ=255)
                                                                                     tCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa    <  1:2447017/93‑1 (MQ=255)
                                                                                     tCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa    <  1:3276742/93‑1 (MQ=255)
                                                                                       aTGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACg  >  2:2361470/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACG  >  NC_000913/488827‑489006

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