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.3% 466.6 / NA 149E1053G (GAG→GGG) aefAputative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (80/68);  minor base T (0/1);  total (80/69)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

TGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACG  >  NC_000913/488827‑489006
                                                                                          |                                                                                         
tGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGc                                                                                         <  1:2926714/93‑1 (MQ=255)
    ggcggcGACTGAGCACCCAAGGGAGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCg                                                                                     >  1:1968765/1‑93 (MQ=255)
      cggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTc                                                                                   <  2:522554/93‑1 (MQ=255)
      cggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTc                                                                                   >  1:2241499/1‑93 (MQ=255)
       ggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCt                                                                                  >  2:4390502/1‑93 (MQ=255)
        gcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTg                                                                                 >  1:2077972/1‑93 (MQ=255)
          gACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTa                                                                               >  1:1855085/1‑93 (MQ=255)
          gACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTa                                                                               >  2:3811695/1‑93 (MQ=255)
            atGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtg                                                                             <  1:4617728/92‑1 (MQ=255)
             tGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgt                                                                            >  2:4385400/1‑93 (MQ=255)
             tGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgt                                                                            <  2:2743015/93‑1 (MQ=255)
              gAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgtg                                                                           >  1:3325536/1‑93 (MQ=255)
                gCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCg                                                                         >  2:2804939/1‑93 (MQ=255)
                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAAc                                                                    <  2:2856918/93‑1 (MQ=255)
                        gggTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGc                                                                 >  2:225737/1‑93 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGt                                                               >  1:1546181/1‑93 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGt                                                               <  1:2198742/93‑1 (MQ=255)
                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGt                                                               >  2:4319220/1‑93 (MQ=255)
                            gATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGa                                                             >  1:4602523/1‑93 (MQ=255)
                            gATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGa                                                             >  2:2381426/1‑93 (MQ=255)
                             aTGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGAc                                                            <  1:2027703/93‑1 (MQ=255)
                              tGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGAcc                                                           <  1:1675307/93‑1 (MQ=255)
                                cACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGt                                                         <  2:3021686/93‑1 (MQ=255)
                                cACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGt                                                         >  1:323328/1‑93 (MQ=255)
                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAg                                                       <  2:3662862/93‑1 (MQ=255)
                                   gAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGt                                                      <  1:3685358/93‑1 (MQ=255)
                                    aaCCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTc                                                     >  1:2942547/1‑93 (MQ=255)
                                     aCCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCg                                                    <  2:2121496/93‑1 (MQ=255)
                                       cAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTa                                                  >  1:1469821/1‑93 (MQ=255)
                                        aaTGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTAc                                                 <  1:554812/93‑1 (MQ=255)
                                        aaTGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTAc                                                 >  1:2589434/1‑93 (MQ=255)
                                             cGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCg                                            <  1:1343546/93‑1 (MQ=255)
                                              ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                           <  1:4299959/93‑1 (MQ=255)
                                              ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                           <  2:2360432/93‑1 (MQ=255)
                                               gAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGAt                                          >  2:2022727/1‑93 (MQ=255)
                                                aaGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATg                                         <  2:3489940/93‑1 (MQ=255)
                                                 aggcttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGa                                        >  1:2327189/4‑93 (MQ=255)
                                                 aGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGGCCGTAGTCGTACTGTCGATGa                                        <  2:4100481/93‑1 (MQ=255)
                                                 aGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGa                                        >  2:1472327/1‑93 (MQ=255)
                                                 aGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACAGTAGTCGTACTGTCGAtaa                                        >  1:3082573/1‑91 (MQ=255)
                                                  gTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGCGAACTGCGTGACCGTAGTCGTACTGTCGATCAg                                       >  1:3853930/1‑93 (MQ=255)
                                                   tcttctTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGc                                      >  1:2340639/1‑93 (MQ=255)
                                                   tcttctTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGc                                      >  1:321108/1‑93 (MQ=255)
                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCt                                     <  2:60819/93‑1 (MQ=255)
                                                    cttctaTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCt                                     <  2:3274113/93‑1 (MQ=255)
                                                     ttcttTACGGCATTTGGTGCCATCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                    <  2:620240/93‑1 (MQ=255)
                                                     ttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                    >  1:1361181/1‑93 (MQ=255)
                                                     ttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                    >  1:3504793/1‑93 (MQ=255)
                                                      tcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCgatgagatga                                   >  1:1738227/1‑93 (MQ=255)
                                                       cttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGaa                                  <  2:3315598/93‑1 (MQ=255)
                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 >  1:1089293/1‑93 (MQ=255)
                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 >  2:3136500/1‑93 (MQ=255)
                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 <  2:3710669/93‑1 (MQ=255)
                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 >  2:31997/1‑93 (MQ=255)
                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 >  2:1918873/1‑93 (MQ=255)
                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACAGTAGTCGTACTGTCGATGAGCTGAAc                                 >  2:1249650/1‑93 (MQ=255)
                                                            cGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTa                             >  2:1420989/1‑93 (MQ=255)
                                                             ggCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTAc                            >  1:3030203/1‑93 (MQ=255)
                                                             ggCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTAc                            <  1:3600926/93‑1 (MQ=255)
                                                             ggCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTAc                            <  1:1086259/93‑1 (MQ=255)
                                                               gaTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTa                          <  1:3279584/92‑1 (MQ=255)
                                                               cATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTa                          <  2:2628418/93‑1 (MQ=255)
                                                                aTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTAt                         >  2:4621615/1‑93 (MQ=255)
                                                                aTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTAt                         >  1:3871843/1‑93 (MQ=255)
                                                                aTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTAt                         >  1:4491040/1‑93 (MQ=255)
                                                                aTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTAt                         <  1:4388198/93‑1 (MQ=255)
                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        <  1:1774571/93‑1 (MQ=255)
                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        <  1:3584794/93‑1 (MQ=255)
                                                                   tGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGa                      <  1:19235/93‑1 (MQ=255)
                                                                   tGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGa                      <  2:2442811/93‑1 (MQ=255)
                                                                   tGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGa                      <  2:1284568/93‑1 (MQ=255)
                                                                     gTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc                    <  1:4183433/93‑1 (MQ=255)
                                                                     gTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc                    <  2:4644230/93‑1 (MQ=255)
                                                                     gTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc                    >  1:1551697/1‑93 (MQ=255)
                                                                     gTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc                    <  2:1795820/93‑1 (MQ=255)
                                                                      tGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCa                   >  1:348381/1‑93 (MQ=255)
                                                                      tGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCa                   <  2:1372114/93‑1 (MQ=255)
                                                                        ccAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc                 <  2:1870150/93‑1 (MQ=255)
                                                                        ccAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc                 <  2:2971498/93‑1 (MQ=255)
                                                                        ccAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc                 <  2:4470526/93‑1 (MQ=255)
                                                                        ccAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc                 <  1:2414154/93‑1 (MQ=255)
                                                                        ccAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc                 <  1:2008842/93‑1 (MQ=255)
                                                                        ccAGCACGTTGGATCATGGGCTGCGGCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc                 <  2:1698396/93‑1 (MQ=255)
                                                                          aGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtg               >  2:142584/1‑93 (MQ=255)
                                                                          aGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtg               >  2:4275695/1‑93 (MQ=255)
                                                                           gCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtgt              <  2:590674/93‑1 (MQ=255)
                                                                           gCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtgt              >  1:2425481/1‑93 (MQ=255)
                                                                               gTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCg          >  2:2538956/1‑93 (MQ=255)
                                                                               gTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCg          <  2:4408446/93‑1 (MQ=255)
                                                                                ttGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGt         <  2:4206372/93‑1 (MQ=255)
                                                                                ttGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGt         >  1:3418828/1‑93 (MQ=255)
                                                                                  ggATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGa       >  1:627247/1‑93 (MQ=255)
                                                                                   gATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaa      >  2:161717/1‑93 (MQ=255)
                                                                                   gATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaa      >  1:538987/1‑93 (MQ=255)
                                                                                   gATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaa      <  1:250730/93‑1 (MQ=255)
                                                                                    aTCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaa     >  1:3763227/1‑93 (MQ=255)
                                                                                     tCATGGGCTGCGTCTGTGTGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa    >  2:3441417/1‑93 (MQ=255)
                                                                                     tCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCGGTGaaaa    >  1:1373195/1‑93 (MQ=255)
                                                                                      cATGTGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAAc   <  2:716777/93‑1 (MQ=255)
                                                                                       aTGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACg  >  2:1171082/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACG  >  NC_000913/488827‑489006

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: