Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,608,707 C→T 16.1% 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,7070CT16.1% 284.9 / 39.9 124D106N (GAT→AAT) yneEbestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (60/43);  new base T (9/11);  total (70/54)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

Only 100 of 129 total aligned reads displayed.

AGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATT  >  NC_000913/1608609‑1608806
                                                                                                  |                                                                                                   
aGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGATAAACTCCCTTACACTTGCCGAATcc                                                                                                    >  1:49884/1‑100 (MQ=255)
 gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCTTTACACTTGCCGAATc                                                                                                     >  2:2538196/1‑98 (MQ=255)
  ccAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCgg                                                                                                  <  1:2467567/100‑1 (MQ=255)
  ccAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCgg                                                                                                  <  1:2415896/100‑1 (MQ=255)
   cAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGc                                                                                                 <  1:3711143/100‑1 (MQ=255)
    aGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCa                                                                                                <  1:406678/100‑1 (MQ=255)
     gCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCaa                                                                                               >  1:2035449/1‑100 (MQ=255)
      cACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaag                                                                                              <  2:2060022/100‑2 (MQ=255)
      cACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaag                                                                                              <  2:1858134/100‑2 (MQ=255)
      cACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAAt                                                                                              <  2:2428626/100‑1 (MQ=255)
      cACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAAt                                                                                              >  1:3630569/1‑100 (MQ=255)
        cTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtn                                                                                            <  2:2341945/100‑4 (MQ=255)
         ttCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTc                                                                                           >  2:3170169/1‑100 (MQ=255)
          tCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGCAAGCGATTTGCAGCCGGGCAAACTACCTTACACTTGCCGAATCCGGCAATg                                                                                             >  2:2285536/1‑97 (MQ=255)
            cGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgtaag                                                                                        <  2:1310972/100‑8 (MQ=255)
              ccTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtttt                                                                                      <  2:1958379/100‑1 (MQ=255)
              ccTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGAGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTAACACTTGCCGAATCCGGCAATGTCGtgt                                                                                       >  2:111606/1‑97 (MQ=255)
              ccTGGGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgtaaggg                                                                                      <  2:2725912/100‑10 (MQ=255)
               cTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGttttt                                                                                     <  2:767943/100‑1 (MQ=255)
                tgtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTa                                                                                    <  1:902307/100‑1 (MQ=255)
                tgtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTa                                                                                    <  1:3168045/100‑1 (MQ=255)
                tgtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTa                                                                                    >  1:3619763/1‑100 (MQ=255)
                tgtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTa                                                                                    <  1:433492/100‑1 (MQ=255)
                tgtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACACCCTTACACTTGCCGAATCCGGCGATGTCGttttt                                                                                     >  2:453822/1‑99 (MQ=255)
                 gtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAc                                                                                   <  1:3498871/100‑1 (MQ=255)
                 gtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAc                                                                                   <  1:2894104/100‑1 (MQ=255)
                 gtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAc                                                                                   <  1:1859329/100‑1 (MQ=255)
                  tgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAcc                                                                                  >  1:546131/1‑100 (MQ=255)
                  tgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAcc                                                                                  <  1:2376873/100‑1 (MQ=255)
                   ggCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAcc                                                                                  >  2:6816/1‑99 (MQ=255)
                   ggCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCt                                                                                 >  1:473333/1‑100 (MQ=255)
                    gCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATCCGGCAATGTCGTTTTTACCTc                                                                                >  2:1870755/1‑100 (MQ=255)
                    gCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTc                                                                                >  1:1410815/1‑100 (MQ=255)
                    gCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTc                                                                                <  1:3815943/100‑1 (MQ=255)
                      tGTTTGCGTGGTGTAATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTGCCTTACAATTGCGGAATCCGGCAATGTCGTGTTTACCTc                                                                                >  2:2257234/1‑98 (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                                                                            >  1:2363923/1‑100 (MQ=255)
                        tttGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCgcgc                                                                            <  2:2538532/100‑1 (MQ=255)
                         ttGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCgcgc                                                                            <  1:2453382/99‑1 (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)
                                      tGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGCAACGACCg                                                                 >  2:530250/1‑97 (MQ=255)
                                       gCGTACACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAg                                                             >  2:20031/1‑100 (MQ=255)
                                       gCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAg                                                             <  2:803110/100‑1 (MQ=255)
                                        cGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgg                                                            <  2:899211/100‑1 (MQ=255)
                                        cGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgg                                                            <  1:1480173/100‑1 (MQ=255)
                                          tAAACAGTGGGCGTAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCa                                                          >  1:2564593/1‑100 (MQ=255)
                                               aGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatcat                                                     >  1:3357583/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)
                                                       aaaGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGcc                                             >  2:3086015/1‑100 (MQ=255)
                                                          gCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGccccc                                          >  1:3112429/1‑100 (MQ=255)
                                                              tttGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAg                                      <  1:2991188/100‑1 (MQ=255)
                                                              tttGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAg                                      <  1:1756186/100‑1 (MQ=255)
                                                              tttGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAg                                      <  1:2652045/100‑1 (MQ=255)
                                                              tttGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAg                                      <  1:720274/100‑1 (MQ=255)
                                                                  cAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCGGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtttt                                  >  1:3338844/1‑100 (MQ=255)
                                                                  cAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtttt                                  >  1:2357436/1‑100 (MQ=255)
                                                                  cAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtttt                                  >  1:1821940/1‑100 (MQ=255)
                                                                    gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                                <  1:702286/100‑1 (MQ=255)
                                                                    gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                                <  1:1506513/100‑1 (MQ=255)
                                                                    gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                                <  1:3154501/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)
                                                                                cccTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCACCAAAGATTTCGCGCTTCAACGTa                    >  2:564119/1‑100 (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)
                                                                                  aacggcattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACAGTGAGGCAATCATCAAATACGCACAAGGTTCTCGAGCGTCAACGTata                  >  2:3412068/8‑98 (MQ=255)
                                                                                   ttACACTTGCCGAATCCGGCCATGTCGTTTTTACTTCGCGCAGTACCGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTaa                   >  2:1613485/1‑98 (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)
                                                                                      cacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGc              >  1:3294145/1‑100 (MQ=255)
                                                                                      cacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGc              >  1:6933/1‑100 (MQ=255)
                                                                                       acgtGCCGAATCCGGCAATGTCGTCTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGACCCCAAAGTTTTCGCGCTTCAACGTAACggg               >  2:2238204/4‑98 (MQ=255)
                                                                                        cTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGt            >  1:2645583/1‑100 (MQ=255)
                                                                                        cTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTATTCGCGCTTCAACGTAACGGGCGt            >  2:2428999/1‑100 (MQ=255)
                                                                                         ttGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTa           >  1:1520616/1‑100 (MQ=255)
                                                                                           gCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc         >  1:2580429/1‑100 (MQ=255)
                                                                                           gCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc         >  2:1550337/1‑100 (MQ=255)
                                                                                             cGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCg       <  1:1882704/100‑1 (MQ=255)
                                                                                             cGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCg       >  1:1419531/1‑100 (MQ=255)
                                                                                                  ccGGCAATGTCGTTTTTTTTTTTCTTTGTAACGACCGTGAGGCAATAATCAACTGCCCCCAAAGTTTTCGCGCTTCATCGTAACGGGCGTACCCGGCat   >  2:2117125/1‑99 (MQ=255)
                                                                                                  ccGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCatt  >  1:1389235/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
AGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATT  >  NC_000913/1608609‑1608806

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