Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,608,707 C→T 15.8% D106N (GAT→AAT)  yneE ← bestrophin family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,608,7070CT15.8% 364.6 / 51.3 152D106N (GAT→AAT) yneEbestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (69/59);  new base T (16/8);  total (85/67)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 155 total aligned reads displayed.

GCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATT  >  NC_000913/1608610‑1608806
                                                                                                 |                                                                                                   
gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTAGACTTGCCGAATCCg                                                                                                   >  1:3000208/1‑100 (MQ=255)
 ccAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCgg                                                                                                  >  2:2523298/1‑100 (MQ=255)
 ccAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCgg                                                                                                  <  2:1458442/100‑1 (MQ=255)
 ccAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCgg                                                                                                  >  2:3506268/1‑100 (MQ=255)
  cAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGc                                                                                                 <  2:3108133/100‑1 (MQ=255)
    gCACTTCCGCGTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGGGAAAGCGATTTGCAGCCGGGCTAACTGCCTTACATTTGCCGAATCCAGCa                                                                                                >  1:3814885/1‑99 (MQ=255)
    gCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCaa                                                                                               >  2:1308145/1‑100 (MQ=255)
    gCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCaa                                                                                               <  2:2445274/100‑1 (MQ=255)
      aCTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATg                                                                                             <  1:3523310/100‑1 (MQ=255)
       cTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGt                                                                                            >  2:1268838/1‑100 (MQ=255)
         tCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCTGCaagtgta                                                                                          <  2:2191098/100‑6 (MQ=255)
         tCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCg                                                                                          <  1:3506244/100‑1 (MQ=255)
         tCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGACAAACTCCCTTACACTTGCCGGATTCGGCaagtgta                                                                                          <  2:1292505/100‑6 (MQ=255)
          ccGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgtaa                                                                                         <  1:2423236/100‑7 (MQ=255)
           cGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACCCTTGCCGGATTCGGCaagtgtaag                                                                                        <  1:296209/100‑8 (MQ=255)
           cGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtt                                                                                        >  1:2536330/1‑100 (MQ=255)
            gCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGatt                                                                                       >  1:109637/1‑100 (MQ=255)
             ccTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtttt                                                                                      <  1:844031/100‑1 (MQ=255)
              cTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGttttt                                                                                     <  2:1830879/100‑1 (MQ=255)
               tgtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGttttt                                                                                     >  1:1177140/1‑99 (MQ=255)
                gtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAc                                                                                   <  1:500781/100‑1 (MQ=255)
                gtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAc                                                                                   <  2:3576941/100‑1 (MQ=255)
                 tgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAcc                                                                                  <  2:1776114/100‑1 (MQ=255)
                    cTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcg                                                                               >  2:2156055/1‑100 (MQ=255)
                       tttGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCgcgc                                                                            <  2:1666098/100‑1 (MQ=255)
                       tttGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCgcgc                                                                            <  1:1982050/100‑1 (MQ=255)
                        ttGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCa                                                                           <  1:1006044/100‑1 (MQ=255)
                         tGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAg                                                                          >  2:2876761/1‑100 (MQ=255)
                         tGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAg                                                                          >  2:3768917/1‑100 (MQ=255)
                          gCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGt                                                                         >  1:2542256/1‑100 (MQ=255)
                          gCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGt                                                                         <  2:3264661/100‑1 (MQ=255)
                           cGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTa                                                                        >  2:382427/1‑100 (MQ=255)
                            gTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTaa                                                                       >  2:1348327/1‑100 (MQ=255)
                              aaTGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACg                                                                     <  2:4414902/100‑1 (MQ=255)
                               aTGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACg                                                                     <  1:2464667/99‑1 (MQ=255)
                               aTGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGa                                                                    >  2:482769/1‑100 (MQ=255)
                               aTGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGa                                                                    >  2:780045/1‑100 (MQ=255)
                                tGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGAc                                                                   <  2:2536343/100‑1 (MQ=255)
                                 gTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGAcc                                                                  >  1:650876/1‑100 (MQ=255)
                                  tCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCg                                                                 <  1:2463840/100‑1 (MQ=255)
                                      gCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAg                                                             <  2:2430548/100‑1 (MQ=255)
                                        gTAAACAGTGGGCGAAAGCGATTTGCAGCCGGTCAAACTACCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGc                                                           <  1:3091209/100‑1 (MQ=255)
                                        gTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgg                                                            <  1:1031990/99‑1 (MQ=255)
                                        gTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGc                                                           >  1:3295134/1‑100 (MQ=255)
                                        gTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGc                                                           <  1:3900976/100‑1 (MQ=255)
                                             cAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatca                                                      <  2:4154494/100‑1 (MQ=255)
                                             cAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatca                                                      <  2:1965572/100‑1 (MQ=255)
                                             cAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACACCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAGCTCGCGCAGTAACGACCGTGAGGCAatca                                                      >  1:2651776/1‑100 (MQ=255)
                                             cAGTGGGCGAAAGCGATTTGCAGCCAGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAat                                                        >  1:2680996/1‑98 (MQ=255)
                                               gTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatcatc                                                    <  1:1873676/100‑1 (MQ=255)
                                               gTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatcatc                                                    >  2:2593463/1‑100 (MQ=255)
                                               gTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatcatc                                                    >  2:439356/1‑100 (MQ=255)
                                                tGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAAtcatca                                                   <  1:3055127/100‑1 (MQ=255)
                                                   gCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACt                                                >  2:4291093/1‑100 (MQ=255)
                                                   gCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACt                                                >  2:3179114/1‑100 (MQ=255)
                                                       aaGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGcc                                             >  1:1763673/1‑99 (MQ=255)
                                                         gCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGccccc                                          >  2:3568448/1‑100 (MQ=255)
                                                             tttGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAg                                      <  2:2498927/100‑1 (MQ=255)
                                                               tGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtt                                    >  2:2815675/1‑100 (MQ=255)
                                                               tGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtt                                    >  2:2916025/1‑100 (MQ=255)
                                                               tGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtt                                    <  2:295613/100‑1 (MQ=255)
                                                               tGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtt                                    >  1:2418181/1‑100 (MQ=255)
                                                                gCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGttt                                   >  2:2699193/1‑100 (MQ=255)
                                                                     cGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcgcg                              >  2:2603999/1‑100 (MQ=255)
                                                                     cGGGCAAACTCCCGTACACTTGCCGAAGCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcgc                               <  1:820244/99‑1 (MQ=255)
                                                                       ggCAAACTCCCTTACACTTGCCGAATCCGGCGATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCACCCAAAGTTTTCgcgc                             >  1:4258061/1‑99 (MQ=255)
                                                                        gCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCtt                           >  2:2163024/1‑100 (MQ=255)
                                                                        gCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCtt                           <  1:2163017/100‑1 (MQ=255)
                                                                        gCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCtt                           >  2:3159550/1‑100 (MQ=255)
                                                                             cTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACg                      <  1:2228316/100‑1 (MQ=255)
                                                                               cccTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTa                    <  1:2523285/100‑1 (MQ=255)
                                                                               cccTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTa                    >  2:298186/1‑100 (MQ=255)
                                                                               aaaacgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTa                    >  1:3231066/10‑100 (MQ=255)
                                                                                 aacgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGCAAc                  >  2:781129/8‑100 (MQ=255)
                                                                                  ttACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACg                 >  1:3868430/1‑100 (MQ=255)
                                                                                  acgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACg                 >  2:366772/7‑100 (MQ=255)
                                                                                   cgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACgg                >  2:4199571/6‑100 (MQ=255)
                                                                                   cgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACgg                >  2:1035584/6‑100 (MQ=255)
                                                                                   cgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACgg                >  2:712183/6‑100 (MQ=255)
                                                                                   cgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTAAGGCAATCATCAACTGCCCACAAAGTTTTCGCGCTTCAACGTAGCgg                >  2:641123/6‑100 (MQ=255)
                                                                                    gacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACggg               >  2:3413203/5‑100 (MQ=255)
                                                                                    gacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACggg               >  1:575455/5‑100 (MQ=255)
                                                                                    acacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACggg               <  1:4307729/100‑1 (MQ=255)
                                                                                    acacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACggg               >  2:2495559/1‑100 (MQ=255)
                                                                                     cacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGc              >  2:1865031/1‑100 (MQ=255)
                                                                                     acattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTAAGGCAATCATCAACTACCCCCAAAGTTTTCGCGCTTCAACGTAACGGGc              >  1:1861127/4‑100 (MQ=255)
                                                                                       cTTTCCGAATCCTACAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAATGGGCg             <  1:4194447/99‑1 (MQ=255)
                                                                                       cTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGt            >  1:3622683/1‑100 (MQ=255)
                                                                                       cTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGt            >  2:3600262/1‑100 (MQ=255)
                                                                                       cTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGt            >  2:2653434/1‑100 (MQ=255)
                                                                                       cTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCACAAAGTTTTCGGGCTTCAACGTAACGGGCg             >  1:763640/1‑99 (MQ=255)
                                                                                          gCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc         >  2:2732098/1‑100 (MQ=255)
                                                                                          gCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc         >  2:2473382/1‑100 (MQ=255)
                                                                                           ccGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAccc        >  1:1639976/1‑100 (MQ=255)
                                                                                            cGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCg       <  1:455414/100‑1 (MQ=255)
                                                                                            cGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCg       >  2:855670/1‑100 (MQ=255)
                                                                                               aTCCGGGAATGTCGTTTTTACCTCGCGCAGTAACGTCCGTGAGGCAATCATCAACTGCCCCGAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCTGGCa    <  1:1432503/100‑1 (MQ=255)
                                                                                               aTCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCa    <  1:1432400/100‑1 (MQ=255)
                                                                                                 ccGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCatt  >  2:4164679/1‑100 (MQ=255)
                                                                                                 ccGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTACCCCCAAAGTTTTCGCGCTTCAACGTAACGAGCCTACCCGGCatt  >  1:356890/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATT  >  NC_000913/1608610‑1608806

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