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,9170AG100.0% 893.6 / NA 281E1053G (GAG→GGG) aefAputative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (118/163);  total (118/163)

Only 100 of 285 total aligned reads displayed.

TGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACG  >  NC_000913/488827‑489006
                                                                                          |                                                                                         
tGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGc                                                                                         >  1:2521926/1‑93 (MQ=255)
tGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGc                                                                                         >  1:1991458/1‑93 (MQ=255)
 gAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCt                                                                                        >  2:6706417/1‑93 (MQ=255)
      cggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTc                                                                                   <  1:6234382/93‑1 (MQ=255)
        gcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTg                                                                                 >  2:3003512/1‑93 (MQ=255)
          gACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTa                                                                               <  1:5750479/93‑1 (MQ=255)
           aCTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAt                                                                              <  2:2811537/93‑1 (MQ=255)
           aCTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAt                                                                              <  1:4287023/93‑1 (MQ=255)
           aCTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAt                                                                              <  1:2076344/93‑1 (MQ=255)
            cTGAGCACCCAAGTGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTTCCAGCACGTTGGATCATGGGCTGCGTCTGTAtg                                                                             <  2:3482295/93‑1 (MQ=255)
            cTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtg                                                                             <  2:2008089/93‑1 (MQ=255)
            cTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtg                                                                             <  1:7171948/93‑1 (MQ=255)
            cTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtg                                                                             <  1:4005555/93‑1 (MQ=255)
             tgcgCACCCAAGGGTGATGCACGAACCAATGCCGGATGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgt                                                                            <  2:7024209/90‑1 (MQ=255)
             tGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgt                                                                            >  1:5081576/1‑93 (MQ=255)
             tGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgt                                                                            <  2:2318763/93‑1 (MQ=255)
              gAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgtg                                                                           <  2:5289741/93‑1 (MQ=255)
               aGCCCCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGc                                                                          <  2:3082558/93‑1 (MQ=255)
                 cACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGt                                                                        >  1:3201398/1‑93 (MQ=255)
                   cccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGa                                                                      >  2:2751564/1‑93 (MQ=255)
                    ccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGaa                                                                     >  2:5208569/1‑93 (MQ=255)
                    ccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGaa                                                                     >  2:2030818/1‑93 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAAc                                                                    <  2:1410454/93‑1 (MQ=255)
                      aaGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTGCGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACt                                                                   >  2:3941592/1‑93 (MQ=255)
                      aaGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACt                                                                   <  2:7148735/93‑1 (MQ=255)
                       aGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTg                                                                  >  1:1685014/1‑93 (MQ=255)
                       aGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTg                                                                  >  2:7362262/1‑93 (MQ=255)
                         ggTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCg                                                                <  1:6559799/93‑1 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGt                                                               <  2:6094339/93‑1 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGt                                                               <  1:5292088/93‑1 (MQ=255)
                             aTGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGAc                                                            >  1:7281263/1‑93 (MQ=255)
                               gCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCg                                                          <  1:2168954/93‑1 (MQ=255)
                                 aCGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTa                                                        <  2:5892472/93‑1 (MQ=255)
                                 aCGAACCAATGCCGGAAGTCTTCTGTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTa                                                        <  1:3718825/93‑1 (MQ=255)
                                      ccAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGt                                                   <  1:9887/93‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTa                                                  >  2:3595832/1‑93 (MQ=255)
                                        aaTGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTAc                                                 <  1:101169/93‑1 (MQ=255)
                                         aTGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACt                                                <  2:7465873/93‑1 (MQ=255)
                                         aTGCCGGAAGTCTTCTTTACGGAATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACt                                                <  1:4152586/93‑1 (MQ=255)
                                             cGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCg                                            >  2:6622921/1‑93 (MQ=255)
                                              ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAGCTGCGTGACCGTAGTCGTACTGTCGa                                           >  1:547788/1‑93 (MQ=255)
                                              ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                           <  2:3766282/93‑1 (MQ=255)
                                              ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                           <  1:1051615/93‑1 (MQ=255)
                                               gAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGAt                                          <  1:7385688/93‑1 (MQ=255)
                                                aaGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATg                                         <  2:3532015/93‑1 (MQ=255)
                                                aaGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATg                                         <  2:815394/93‑1 (MQ=255)
                                                 aGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGa                                        <  1:2734708/93‑1 (MQ=255)
                                                   tcttctTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGc                                      >  2:1644487/1‑93 (MQ=255)
                                                   tcttctTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGc                                      <  2:6337343/93‑1 (MQ=255)
                                                   tcttctTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGc                                      <  1:6396817/93‑1 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGAGCGTAGTCGTACTGTCGATGAGCt                                     >  2:4319422/1‑93 (MQ=255)
                                                    cttattTACGGCATTTGGTGCCAGCACGTTGGATCATTGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCt                                     <  2:6569260/93‑1 (MQ=255)
                                                     ttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                    >  1:248858/1‑93 (MQ=255)
                                                     ttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                    <  2:5727247/93‑1 (MQ=255)
                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 >  1:4804363/1‑93 (MQ=255)
                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 <  2:3881019/93‑1 (MQ=255)
                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 <  1:4025250/93‑1 (MQ=255)
                                                          tACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCg                               <  1:6895376/93‑1 (MQ=255)
                                                            cGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTa                             <  2:3212849/93‑1 (MQ=255)
                                                             ggCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTAc                            <  2:4986008/93‑1 (MQ=255)
                                                             ggCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTAc                            <  1:1863983/93‑1 (MQ=255)
                                                               cATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTa                          <  2:4718836/93‑1 (MQ=255)
                                                               cATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTa                          >  2:7516767/1‑93 (MQ=255)
                                                               cATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTa                          <  1:1347715/93‑1 (MQ=255)
                                                                aTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTAt                         >  1:3270886/1‑93 (MQ=255)
                                                                aTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTAt                         <  2:7893227/93‑1 (MQ=255)
                                                                aTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTAt                         >  1:2873449/1‑93 (MQ=255)
                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        <  1:5821005/93‑1 (MQ=255)
                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        >  1:1543149/1‑93 (MQ=255)
                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        <  2:2980682/93‑1 (MQ=255)
                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        >  1:6212157/1‑93 (MQ=255)
                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        >  2:2689869/1‑93 (MQ=255)
                                                                  ttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCg                       >  2:7734769/1‑93 (MQ=255)
                                                                   tGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGa                      <  1:7791753/93‑1 (MQ=255)
                                                                    ggTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGAt                     <  2:5519820/93‑1 (MQ=255)
                                                                     gTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc                    >  2:3378766/1‑93 (MQ=255)
                                                                      tGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCa                   <  2:372641/93‑1 (MQ=255)
                                                                      tGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCa                   <  1:7883863/93‑1 (MQ=255)
                                                                        ccAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc                 >  2:4537685/1‑93 (MQ=255)
                                                                         cAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCt                <  1:3535259/93‑1 (MQ=255)
                                                                          aGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtg               >  1:4551288/1‑93 (MQ=255)
                                                                           gCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtgt              >  2:1184717/1‑93 (MQ=255)
                                                                           gCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtgt              >  1:925908/1‑93 (MQ=255)
                                                                             ccGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGc            <  2:2468044/92‑1 (MQ=255)
                                                                              cGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGcc           <  2:1079902/93‑1 (MQ=255)
                                                                               gTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCg          >  1:1616100/1‑93 (MQ=255)
                                                                               gTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCg          >  2:4364899/1‑93 (MQ=255)
                                                                                 tGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTg        <  2:4869115/93‑1 (MQ=255)
                                                                                  ggATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGa       <  1:1845818/93‑1 (MQ=255)
                                                                                  ggATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGa       <  1:7215762/93‑1 (MQ=255)
                                                                                  ggATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGa       <  2:6517092/93‑1 (MQ=255)
                                                                                    aTCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaa     >  1:3319104/1‑93 (MQ=255)
                                                                                     tCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGGCGTGaaca    >  2:599302/1‑91 (MQ=255)
                                                                                     tCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa    <  2:1806982/93‑1 (MQ=255)
                                                                                     tCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa    >  2:4200398/1‑93 (MQ=255)
                                                                                     tCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa    <  1:4889336/93‑1 (MQ=255)
                                                                                       aTGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACg  >  2:1505438/1‑93 (MQ=255)
                                                                                       aTGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACg  >  2:5795820/1‑93 (MQ=255)
                                                                                       aTGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACg  <  1:6063948/93‑1 (MQ=255)
                                                                                       aTGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACg  <  1:4422946/93‑1 (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: