Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,608,704 A→G 16.5% 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,7040AG16.5% 271.0 / 46.3 127S107P (TCG→CCG) yneEbestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (62/44);  new base G (10/11);  total (72/55)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-02

Only 100 of 127 total aligned reads displayed.

AGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCG  >  NC_000913/1608609‑1608801
                                                                                               |                                                                                                 
aGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATcc                                                                                               <  2:463973/100‑1 (MQ=255)
aGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGATAAACTCCCTTACACTTGCCGAATcc                                                                                               >  1:49884/1‑100 (MQ=255)
 gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCTTTACACTTGCCGAATc                                                                                                >  2:2538196/1‑98 (MQ=255)
  ccAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCgg                                                                                             <  1:2415896/100‑1 (MQ=255)
   cAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGc                                                                                            <  1:3711143/100‑1 (MQ=255)
    aGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCa                                                                                           <  1:406678/100‑1 (MQ=255)
    aGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCa                                                                                           <  1:2475551/100‑1 (MQ=255)
      cACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaag                                                                                         <  2:2060022/100‑2 (MQ=255)
      cACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaag                                                                                         <  2:1858134/100‑2 (MQ=255)
      cACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAAt                                                                                         <  2:2428626/100‑1 (MQ=255)
        cTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtn                                                                                       <  2:2341945/100‑4 (MQ=255)
        cTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGt                                                                                       >  1:1624725/1‑100 (MQ=255)
         ttCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTc                                                                                      >  2:3170169/1‑100 (MQ=255)
          tCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGCAAGCGATTTGCAGCCGGGCAAACTACCTTACACTTGCCGAATCCGGCAATg                                                                                        >  2:2285536/1‑97 (MQ=255)
          tCCGCCTGTGGCTGTTGGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtg                                                                                       <  2:2925795/98‑4 (MQ=255)
            cGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgtaag                                                                                   <  1:2711673/100‑8 (MQ=255)
            cGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgtaag                                                                                   <  2:1982742/100‑8 (MQ=255)
            cGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgtaag                                                                                   <  2:1310972/100‑8 (MQ=255)
              ccTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtttt                                                                                 <  2:1958379/100‑1 (MQ=255)
              ccTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGAGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTAACACTTGCCGAATCCGGCAATGTCGtgt                                                                                  >  2:111606/1‑97 (MQ=255)
               cTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGttttt                                                                                <  2:767943/100‑1 (MQ=255)
                tgtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTa                                                                               <  1:902307/100‑1 (MQ=255)
                tgtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTa                                                                               >  1:3314564/1‑100 (MQ=255)
                tgtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTa                                                                               >  1:3619763/1‑100 (MQ=255)
                tgtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTa                                                                               <  1:433492/100‑1 (MQ=255)
                 gtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAc                                                                              <  1:2894104/100‑1 (MQ=255)
                 gtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAc                                                                              <  1:1859329/100‑1 (MQ=255)
                 gtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAc                                                                              <  1:3498871/100‑1 (MQ=255)
                  tgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAcc                                                                             <  1:2376873/100‑1 (MQ=255)
                  tgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAcc                                                                             >  1:546131/1‑100 (MQ=255)
                   ggCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAcc                                                                             >  2:6816/1‑99 (MQ=255)
                   ggCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCt                                                                            >  1:2447652/1‑100 (MQ=255)
                   ggCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCt                                                                            >  1:473333/1‑100 (MQ=255)
                    gCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATCCGGCAATGTCGTTTTTACCTc                                                                           >  2:1870755/1‑100 (MQ=255)
                    gCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTc                                                                           <  1:3815943/100‑1 (MQ=255)
                    gCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTc                                                                           >  1:1410815/1‑100 (MQ=255)
                      ngTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcgc                                                                         >  1:348974/2‑100 (MQ=255)
                      gatctGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcgc                                                                         <  2:1923277/96‑1 (MQ=255)
                       gTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTAGACTTGCCGAATCCGGCAATGTCGTTTTTACCTcgcg                                                                        >  2:2247714/1‑100 (MQ=255)
                       gTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcgcg                                                                        >  2:2991166/1‑100 (MQ=255)
                        tttGCGTAATGTCATGCGTAAACAGTTGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCgcgc                                                                       <  2:2994433/100‑1 (MQ=255)
                        tttGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCgcgc                                                                       <  2:2538532/100‑1 (MQ=255)
                         ttGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCgcgc                                                                       <  1:2453382/99‑1 (MQ=255)
                          tGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAg                                                                     >  2:2479363/1‑100 (MQ=255)
                          tGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCAGAATCCGGCAATGTCGTTTTTACCTAGCGCa                                                                      >  2:803910/1‑99 (MQ=255)
                           gCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcgcagatc                                                                    >  1:1923288/1‑95 (MQ=255)
                             gTAATGTCATGCGTAAACAGTGGGCGACAGCGATTTGCAGCCGGGCAAACTCCCTTCCACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTa                                                                   <  2:1126410/99‑1 (MQ=255)
                              tAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTa                                                                   <  2:3615185/98‑1 (MQ=255)
                                     aTGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTg                                                          <  2:3090308/100‑1 (MQ=255)
                                      tGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGCAACGACCg                                                            >  2:530250/1‑97 (MQ=255)
                                       gCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAg                                                        <  2:803110/100‑1 (MQ=255)
                                        cGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgg                                                       <  2:899211/100‑1 (MQ=255)
                                          tAAACAGTGGGCGTAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCa                                                     >  1:2564593/1‑100 (MQ=255)
                                                  gggCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCaa                                             >  1:3408180/1‑100 (MQ=255)
                                                   ggCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTGTTTACCTCGCGCAGTAACGACCGTGAGGCAATCGTTAAc                                            >  2:2290615/1‑100 (MQ=255)
                                                    gCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACt                                           <  2:546125/100‑1 (MQ=255)
                                                     cGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTg                                          <  1:3578755/100‑1 (MQ=255)
                                                     cGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCCTGAGGCAATCATCAACTg                                          >  2:3261031/1‑100 (MQ=255)
                                                      gAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGc                                         >  1:3345102/1‑100 (MQ=255)
                                                        aaGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGcc                                        >  1:2260618/1‑99 (MQ=255)
                                                          gCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGccccc                                     >  1:3112429/1‑100 (MQ=255)
                                                              tttGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAg                                 <  1:720274/100‑1 (MQ=255)
                                                              tttGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAg                                 <  1:2652045/100‑1 (MQ=255)
                                                              tttGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAg                                 <  1:1756186/100‑1 (MQ=255)
                                                                tGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtt                               >  1:36682/1‑100 (MQ=255)
                                                                  cAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCGGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtttt                             >  1:3338844/1‑100 (MQ=255)
                                                                  cAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtttt                             >  1:2357436/1‑100 (MQ=255)
                                                                  cAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGttt                              >  2:1550366/1‑99 (MQ=255)
                                                                   aGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTc                            <  1:1453110/100‑1 (MQ=255)
                                                                    gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                           <  2:3433727/100‑1 (MQ=255)
                                                                    gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                           <  1:3154501/100‑1 (MQ=255)
                                                                    gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                           <  1:1506513/100‑1 (MQ=255)
                                                                    gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                           <  1:702286/100‑1 (MQ=255)
                                                                    gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTc                            <  2:3768736/99‑1 (MQ=255)
                                                                     ccGGGCAAACTCCCTTACACTTGCCTAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcgc                          <  2:1738591/100‑1 (MQ=255)
                                                                      cGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcgcg                         >  1:3541466/1‑100 (MQ=255)
                                                                        ggCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCt                       >  2:46543/1‑100 (MQ=255)
                                                                         gCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCtt                      <  2:278027/100‑1 (MQ=255)
                                                                           aaaCTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCa                    >  1:3046314/1‑100 (MQ=255)
                                                                              cTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACg                 <  1:2247728/100‑1 (MQ=255)
                                                                                 ccTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACAGTGAGGCAATCATCAACTGCCCACAAAGTTTTCGCGCTTCAACGTa               >  2:729438/1‑99 (MQ=255)
                                                                                 aaacgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTa               >  1:106682/9‑99 (MQ=255)
                                                                                  cTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCCGTCACGACCGTGAGGCCATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAAc             >  2:1525564/1‑100 (MQ=255)
                                                                                  cTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAGGTTTTCGCGCTTCAACGTAAc             >  1:3457621/1‑100 (MQ=255)
                                                                                   acgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACg            >  1:923471/7‑100 (MQ=255)
                                                                                   acgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACg            >  1:2469143/7‑100 (MQ=255)
                                                                                    cgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACgg           >  2:1924450/6‑100 (MQ=255)
                                                                                     gacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAAcag           >  2:338357/5‑97 (MQ=255)
                                                                                     acacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACggg          <  2:2188223/100‑1 (MQ=255)
                                                                                      cacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGc         >  1:3294145/1‑100 (MQ=255)
                                                                                      cacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGc         <  1:3058435/100‑1 (MQ=255)
                                                                                      cacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGc         >  1:6933/1‑100 (MQ=255)
                                                                                       acgtGCCGAATCCGGCAATGTCGTCTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGACCCCAAAGTTTTCGCGCTTCAACGTAACggg          >  2:2238204/4‑98 (MQ=255)
                                                                                         ttGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTa      >  1:1520616/1‑100 (MQ=255)
                                                                                          tGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGATTTCGCGCTTCAACGTAACGGGCGTAc     >  2:639805/1‑100 (MQ=255)
                                                                                           gCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc    >  2:1550337/1‑100 (MQ=255)
                                                                                           gCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc    >  1:2580429/1‑100 (MQ=255)
                                                                                             cGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCg  >  1:1918744/1‑100 (MQ=255)
                                                                                             cGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCg  <  1:1882704/100‑1 (MQ=255)
                                                                                             cGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCg  >  1:1419531/1‑100 (MQ=255)
                                                                                               |                                                                                                 
AGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCG  >  NC_000913/1608609‑1608801

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: