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,9170AG99.6% 858.5 / NA 262E1053G (GAG→GGG) aefAputative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (143/118);  minor base C (0/1);  total (143/119)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

TGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACG  >  NC_000913/488827‑489006
                                                                                          |                                                                                         
tGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGc                                                                                         >  2:1661571/1‑93 (MQ=255)
tGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGc                                                                                         <  1:4831459/93‑1 (MQ=255)
tGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGc                                                                                         >  2:2924230/1‑93 (MQ=255)
tGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGc                                                                                         >  1:5868814/1‑93 (MQ=255)
 gAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCt                                                                                        >  1:2215493/1‑93 (MQ=255)
 gAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCt                                                                                        <  2:2045175/93‑1 (MQ=255)
  aaGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTg                                                                                       >  1:6016842/1‑93 (MQ=255)
  aaGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTg                                                                                       >  2:6205359/1‑93 (MQ=255)
   aGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGc                                                                                      >  1:6590697/1‑93 (MQ=255)
   aGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGc                                                                                      >  2:603381/1‑93 (MQ=255)
    ggcggcGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCg                                                                                     >  2:4006798/1‑93 (MQ=255)
       cgcgACTGAGCACCCAAGGATGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCt                                                                                  <  1:3696743/92‑1 (MQ=255)
        gcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTg                                                                                 >  1:2891775/1‑93 (MQ=255)
        gcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTg                                                                                 >  1:3221326/1‑93 (MQ=255)
         cgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGt                                                                                >  2:2805645/1‑93 (MQ=255)
           aCTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAt                                                                              >  2:5876732/1‑93 (MQ=255)
           aCTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAt                                                                              >  1:541594/1‑93 (MQ=255)
              gAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgtg                                                                           <  2:2154606/93‑1 (MQ=255)
              gAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgtg                                                                           >  1:4383440/1‑93 (MQ=255)
                 cACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGt                                                                        <  1:3637870/93‑1 (MQ=255)
                 cACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGt                                                                        <  1:6663112/93‑1 (MQ=255)
                  aCCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTg                                                                       >  1:3356951/1‑93 (MQ=255)
                  aCCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTg                                                                       <  1:6182745/93‑1 (MQ=255)
                  aCCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTg                                                                       >  1:1545545/1‑93 (MQ=255)
                   cccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGa                                                                      <  1:4447197/93‑1 (MQ=255)
                    ccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGaa                                                                     >  2:4091332/1‑93 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAAc                                                                    <  1:3497508/93‑1 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAAc                                                                    <  2:569593/93‑1 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAAc                                                                    <  2:502109/93‑1 (MQ=255)
                      aaGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACt                                                                   <  2:2362306/93‑1 (MQ=255)
                      aaGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACt                                                                   <  2:4578537/93‑1 (MQ=255)
                       aGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTg                                                                  <  2:2403142/93‑1 (MQ=255)
                       aGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTg                                                                  >  2:1050819/1‑93 (MQ=255)
                         ggTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCg                                                                >  1:6237396/1‑93 (MQ=255)
                         ggTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCg                                                                <  1:61742/93‑1 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGt                                                               <  2:5133910/93‑1 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGt                                                               <  2:3113298/93‑1 (MQ=255)
                             aTGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGAc                                                            >  1:2347134/1‑93 (MQ=255)
                             aTGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGAc                                                            >  1:2784416/1‑93 (MQ=255)
                                   gAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGt                                                      >  2:4262878/1‑93 (MQ=255)
                                      ccAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGt                                                   <  1:387090/93‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTa                                                  <  2:3699426/93‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTa                                                  <  1:1194522/93‑1 (MQ=255)
                                         aTGCCGGAAGTCTTCTTTACGTCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACt                                                <  1:2289150/93‑1 (MQ=255)
                                         aTGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACt                                                <  2:6131963/93‑1 (MQ=255)
                                            ccGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTc                                             <  2:2679788/93‑1 (MQ=255)
                                             cGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCg                                            <  1:2696009/93‑1 (MQ=255)
                                             cGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCg                                            >  1:78293/1‑93 (MQ=255)
                                              ggAAGTCTTCTTTACGGTATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                           <  1:4508394/93‑1 (MQ=255)
                                              ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                           <  1:5953481/93‑1 (MQ=255)
                                              ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                           >  1:4805008/1‑93 (MQ=255)
                                              ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                           >  1:5288278/1‑93 (MQ=255)
                                               gAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGAt                                          >  1:607764/1‑93 (MQ=255)
                                               gAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGAt                                          >  1:4125208/1‑93 (MQ=255)
                                               gAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCAGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGAt                                          <  1:5431085/93‑1 (MQ=255)
                                                aaGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATg                                         <  1:2971201/93‑1 (MQ=255)
                                                aaGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATg                                         <  2:4416540/93‑1 (MQ=255)
                                                 aGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGa                                        <  2:760594/93‑1 (MQ=255)
                                                 aGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGa                                        >  2:6631089/1‑93 (MQ=255)
                                                 aGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGa                                        >  1:3619085/1‑93 (MQ=255)
                                                   tcttctTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGc                                      <  2:1186361/93‑1 (MQ=255)
                                                   tcttctTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGc                                      <  2:4726897/93‑1 (MQ=255)
                                                     ttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                    <  1:1883831/93‑1 (MQ=255)
                                                     ttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                    >  2:4932873/1‑93 (MQ=255)
                                                     ttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                    <  1:4215796/93‑1 (MQ=255)
                                                      tcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGa                                   <  1:1323850/93‑1 (MQ=255)
                                                          tACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCg                               >  2:5164158/1‑93 (MQ=255)
                                                           aCGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGt                              >  2:1394946/1‑93 (MQ=255)
                                                           aCGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGt                              <  2:1875012/93‑1 (MQ=255)
                                                            cGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTa                             >  1:1130588/1‑93 (MQ=255)
                                                            cGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTa                             >  2:3230747/1‑93 (MQ=255)
                                                             ggCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTAc                            >  1:3079498/1‑93 (MQ=255)
                                                              gCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACt                           <  2:5441300/93‑1 (MQ=255)
                                                               cATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTa                          >  1:1459300/1‑93 (MQ=255)
                                                                aTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTAt                         >  2:495983/1‑93 (MQ=255)
                                                                aTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTAt                         <  1:378527/93‑1 (MQ=255)
                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        >  2:251076/1‑93 (MQ=255)
                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        >  1:6393674/1‑93 (MQ=255)
                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        >  1:5616368/1‑93 (MQ=255)
                                                                    ggTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGAt                     >  1:6726254/1‑93 (MQ=255)
                                                                    ggTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGAt                     <  1:5095046/93‑1 (MQ=255)
                                                                     gTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc                    <  1:1239423/93‑1 (MQ=255)
                                                                     gTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc                    <  2:1201996/93‑1 (MQ=255)
                                                                     gTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc                    <  1:2689464/93‑1 (MQ=255)
                                                                      tGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCa                   <  2:5487675/93‑1 (MQ=255)
                                                                        ccAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc                 <  2:3937783/93‑1 (MQ=255)
                                                                        ccAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc                 <  2:6006216/93‑1 (MQ=255)
                                                                          aGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtg               >  1:4546869/1‑93 (MQ=255)
                                                                          aGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtg               <  1:878972/93‑1 (MQ=255)
                                                                           gCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATAAGCTGAACCGTACTATCGATCAGCtgt              >  1:5307463/1‑93 (MQ=255)
                                                                            cACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtgtg             <  1:1661571/93‑1 (MQ=255)
                                                                             aCGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGc            >  2:3330077/1‑93 (MQ=255)
                                                                                ttGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGt         <  2:5353125/93‑1 (MQ=255)
                                                                                 tGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTg        >  1:1858610/1‑93 (MQ=255)
                                                                                  ggATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGa       >  1:156157/1‑93 (MQ=255)
                                                                                   gATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaa      >  2:635706/1‑93 (MQ=255)
                                                                                     tCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa    <  1:6509060/93‑1 (MQ=255)
                                                                                     tCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa    <  2:4439729/93‑1 (MQ=255)
                                                                                       aTGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACg  <  1:5120125/93‑1 (MQ=255)
                                                                                       aTGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACg  >  1:4979160/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACG  >  NC_000913/488827‑489006

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