Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,608,704 A→G 19.0% S107P (TCG→CCG)  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,7040AG19.0% 243.8 / 54.5 121S107P (TCG→CCG) yneEbestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (54/44);  new base G (14/9);  total (68/53)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

Only 100 of 125 total aligned reads displayed.

TGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGC  >  NC_000913/1608607‑1608803
                                                                                                 |                                                                                                   
tGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAAt                                                                                                   <  1:2316534/100‑1 (MQ=255)
  aGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGaa                                                                                                    >  1:895579/1‑97 (MQ=255)
  aGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATcc                                                                                                 >  2:1719772/1‑100 (MQ=255)
   gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTATACTTGCCGAATcgg                                                                                                >  1:1761816/1‑98 (MQ=255)
   gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCg                                                                                                >  2:1992843/1‑100 (MQ=255)
    ccAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCgg                                                                                               >  1:1046802/1‑100 (MQ=255)
    ccAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCgg                                                                                               <  2:2587944/100‑1 (MQ=255)
    ccAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCgg                                                                                               >  2:2500787/1‑100 (MQ=255)
     cAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGc                                                                                              <  2:562370/100‑1 (MQ=255)
       gCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCaa                                                                                            <  1:1530287/100‑1 (MQ=255)
        cACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGcgag                                                                                           <  2:1168021/100‑4 (MQ=255)
        cACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaag                                                                                           <  2:1732022/100‑2 (MQ=255)
         aCTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATg                                                                                          <  1:1799672/100‑1 (MQ=255)
          cTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtg                                                                                         <  1:5792/100‑4 (MQ=255)
          cTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATg                                                                                          >  2:2145278/1‑99 (MQ=255)
           ttCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgt                                                                                        <  2:1710974/100‑5 (MQ=255)
            tCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCg                                                                                       >  1:1878062/1‑100 (MQ=255)
               gCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgtaagg                                                                                    <  1:139698/100‑9 (MQ=255)
               gCCTGGGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgtaagg                                                                                    <  1:1053250/100‑9 (MQ=255)
                ccTGTGGCTGTTTGCTTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGCCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGttt                                                                                    <  1:425258/99‑1 (MQ=255)
                ccTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtttt                                                                                   >  1:2115310/1‑100 (MQ=255)
                 cTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAGACTCCCTTACACTTGCCGAATCCGGCAATGTCGttttt                                                                                  <  1:1595642/100‑1 (MQ=255)
                 cTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGttttt                                                                                  <  1:2115804/100‑1 (MQ=255)
                  tgtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtttt                                                                                   <  2:318170/98‑1 (MQ=255)
                   gtgGCTGTGTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAc                                                                                <  1:2446312/100‑1 (MQ=255)
                    tgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAcc                                                                               <  2:70987/100‑1 (MQ=255)
                        tGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcgc                                                                           >  2:513500/1‑100 (MQ=255)
                         gTTTGCGTGATGGCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGACGGGCAAAATCGCGTACACTAACAGAATACGGCAATGTCGTTTTTa                                                                                 >  1:500655/1‑93 (MQ=255)
                         gTTTGCGTAATGTCATGCGTAAACCGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCCATGTCGTTTTTACCtacc                                                                           >  1:87903/1‑96 (MQ=255)
                          tttGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCgcgc                                                                         <  2:1670141/100‑1 (MQ=255)
                           ttGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGGTTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCa                                                                        <  1:2175823/100‑1 (MQ=255)
                           ttGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcgc                                                                           >  2:78118/1‑97 (MQ=255)
                            tGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCGGAATCCGGCAATGTCGTTTTTACCTCGCGCAg                                                                       >  1:2458529/1‑100 (MQ=255)
                            tGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAg                                                                       <  2:635408/100‑1 (MQ=255)
                             gCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGt                                                                      >  1:1330788/1‑100 (MQ=255)
                                tAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAAc                                                                   <  2:876181/100‑1 (MQ=255)
                                    gTCATGCGTAAACAATGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCGTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGAcc                                                               >  1:2284729/1‑100 (MQ=255)
                                     tCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCg                                                              >  2:2529312/1‑100 (MQ=255)
                                      cATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGt                                                             >  2:1848307/1‑100 (MQ=255)
                                       cgacctAAACAGTGGTCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTg                                                            <  2:2552042/95‑1 (MQ=255)
                                        tGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGa                                                           >  1:2508150/1‑100 (MQ=255)
                                        tGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGa                                                           <  1:513498/100‑1 (MQ=255)
                                        tGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGa                                                           >  2:627523/1‑100 (MQ=255)
                                          cGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgg                                                         >  1:1872959/1‑100 (MQ=255)
                                           gcaaaCAGTGGGCGAAAGCGCCACGCACCCGGGCAAACTGCCTTACACTGACCGAATCCGGCAACGACGTTTTTACCTCGCCCAGTAACGAc                                                                >  1:1787073/3‑92 (MQ=255)
                                           gTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGc                                                        >  2:2256504/1‑100 (MQ=255)
                                            tAAACAGTGGGCGAAAGCGATTTGCAGCCGGGTAAACTCCCTTACACTTACAGAATCCGGCAATGTCGTTTGTAAATCGAGCAGTAAGGACCGTGAGGCa                                                       >  1:2260142/1‑100 (MQ=255)
                                             aaaCAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgatcg                                                      >  1:2552023/1‑96 (MQ=255)
                                               aCAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatc                                                    <  1:2357743/100‑1 (MQ=255)
                                               aCAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATCTCGTTTTTACCTCGCGCAGTAACGCCCGTGAGGCAatc                                                    <  1:1601693/100‑1 (MQ=255)
                                                   tGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAAtcatca                                                <  1:1963567/100‑1 (MQ=255)
                                                   tGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAAtcatca                                                <  2:1296773/100‑1 (MQ=255)
                                                    gggCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAAtcatca                                                >  1:223352/1‑99 (MQ=255)
                                                     ggCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAAc                                              <  1:823340/100‑1 (MQ=255)
                                                       cGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTg                                            >  2:2489255/1‑100 (MQ=255)
                                                         aaaGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGcc                                          <  2:200736/100‑1 (MQ=255)
                                                          aaGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGccc                                         >  2:882402/1‑100 (MQ=255)
                                                           aGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGcccc                                        >  2:195915/1‑100 (MQ=255)
                                                                 ttGCAGCGGGGGACACTCTCTTGCACTGGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAAGGACCGTGAGGCAATCATCAACTGCCCCCAAAGt                                  <  1:1780912/100‑1 (MQ=255)
                                                                  tggaGCCGGGCAAACTCCCTTACACTTGCGGAATCCTGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAg                                   <  1:915091/95‑1 (MQ=255)
                                                                  tGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtt                                 >  2:456747/1‑100 (MQ=255)
                                                                    cAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtttt                               >  1:1482888/1‑100 (MQ=255)
                                                                     aGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTc                              >  2:455499/1‑100 (MQ=255)
                                                                      gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                             >  1:324563/1‑100 (MQ=255)
                                                                      gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                             >  2:956515/1‑100 (MQ=255)
                                                                      gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                             <  1:1281878/100‑1 (MQ=255)
                                                                        cGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcgcg                           >  2:279264/1‑100 (MQ=255)
                                                                          ggCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCt                         <  2:1383746/100‑1 (MQ=255)
                                                                           gCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAAGCAATAATCAACTGCCCCCAAAGTTTTCGCGCt                         >  1:823116/1‑99 (MQ=255)
                                                                              aaCTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCaa                     >  2:2115163/1‑100 (MQ=255)
                                                                              aaCTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCaa                     <  1:627520/100‑1 (MQ=255)
                                                                                cTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTTCCCTCAAAGTTTTCGCGCTTCAACg                   <  1:144491/100‑1 (MQ=255)
                                                                                  cccTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTa                 >  2:645189/1‑100 (MQ=255)
                                                                                   aaacgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCGTCAACTGCCCCCAAAGTATTCGCGCTTCAACGTaa                >  1:1408120/9‑100 (MQ=255)
                                                                                    cTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGAAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAAc               >  1:1417938/1‑100 (MQ=255)
                                                                                    aacgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAAc               >  2:962486/8‑100 (MQ=255)
                                                                                    aacgacattGCCGGATTCGGAAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTaa                >  2:2439891/8‑99 (MQ=255)
                                                                                      cgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACgg             >  2:1414777/6‑100 (MQ=255)
                                                                                       gacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACggg            >  2:336850/5‑100 (MQ=255)
                                                                                       gacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACggg            >  2:1664494/5‑100 (MQ=255)
                                                                                       acacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACggg            <  2:1857045/100‑1 (MQ=255)
                                                                                        cacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACGGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGc           >  1:268113/1‑100 (MQ=255)
                                                                                        cacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGc           >  1:80019/1‑100 (MQ=255)
                                                                                        acattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGc           >  2:1853623/4‑100 (MQ=255)
                                                                                          cTTGCCGAATCCGGCAATGTCGTTTTGACCTAGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTCTTCGCGCTTCAACGTAACGGGCg          >  1:2472724/1‑99 (MQ=255)
                                                                                           ttGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCgag        >  1:865372/1‑98 (MQ=255)
                                                                                            tGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAc       <  2:2020604/100‑1 (MQ=255)
                                                                                             gacGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc      >  2:1747370/3‑100 (MQ=255)
                                                                                             gCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTCACGACCGTGAGGCAAATATCAACTGCCCCCAAAGTTATCGCGCTTCAACGTAACGGGCGTAcc      >  1:1101789/1‑100 (MQ=255)
                                                                                             gCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc      >  2:1992513/1‑100 (MQ=255)
                                                                                             gCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc      <  1:186748/100‑1 (MQ=255)
                                                                                             gCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc      >  2:1796797/1‑100 (MQ=255)
                                                                                              ccGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAccc     <  2:1369360/100‑1 (MQ=255)
                                                                                              ccGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAccc     <  1:1319659/100‑1 (MQ=255)
                                                                                              ccGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAccc     <  1:1319638/100‑1 (MQ=255)
                                                                                               cGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCg    >  2:1591698/1‑100 (MQ=255)
                                                                                                gAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCgg   >  1:1345836/1‑100 (MQ=255)
                                                                                                gAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGACAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCgg   >  2:753230/1‑100 (MQ=255)
                                                                                                 aaTCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGc  <  2:1968159/100‑1 (MQ=255)
                                                                                                 aaTCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGc  <  2:1659320/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
TGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGC  >  NC_000913/1608607‑1608803

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: