Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,608,704 A→G 15.2% 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,7040AG15.2% 222.7 / 32.4 99S107P (TCG→CCG) yneEbestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (48/36);  new base G (7/8);  total (55/44)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

TGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGG  >  NC_000913/1608607‑1608802
                                                                                                 |                                                                                                  
tGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGa                                                                                                    <  2:58509/98‑1 (MQ=255)
tGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAAt                                                                                                  <  2:2444904/100‑1 (MQ=255)
 gAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATc                                                                                                 >  1:1101846/1‑100 (MQ=255)
  aGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATcc                                                                                                >  2:2937157/1‑100 (MQ=255)
   gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCg                                                                                               <  2:1689713/100‑1 (MQ=255)
   gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCg                                                                                               <  2:4050179/100‑1 (MQ=255)
    ccAGCACTTCCGCCTGTGGCTGTTTGCGTACTGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCgg                                                                                              >  2:935817/1‑100 (MQ=255)
    ccAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCgg                                                                                              >  1:1499751/1‑100 (MQ=255)
     cAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGc                                                                                             <  1:3136516/100‑1 (MQ=255)
       gCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCCGCCGTGCAAACTCCCTTACACTTGCCGAATCCGGCa                                                                                            <  2:2100327/99‑1 (MQ=255)
        cACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaan                                                                                          <  1:1963528/100‑2 (MQ=255)
        cACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaag                                                                                          <  2:825163/100‑2 (MQ=255)
        cACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAAt                                                                                          <  1:2364161/100‑1 (MQ=255)
         aCTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagt                                                                                         <  2:937736/100‑3 (MQ=255)
         aCTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagt                                                                                         <  1:2322449/100‑3 (MQ=255)
          cTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtg                                                                                        <  1:3338177/100‑4 (MQ=255)
          cTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGt                                                                                        <  1:2762531/100‑1 (MQ=255)
           ttCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgt                                                                                       <  1:969074/100‑5 (MQ=255)
           ttCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtg                                                                                        <  2:3191770/99‑4 (MQ=255)
            tCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgta                                                                                      >  1:361509/1‑95 (MQ=255)
            tCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCg                                                                                      >  2:724178/1‑100 (MQ=255)
             ccGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCg                                                                                      >  2:2525724/1‑99 (MQ=255)
             ccGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGt                                                                                     <  2:300386/100‑1 (MQ=255)
              cGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagtgtaag                                                                                    <  2:361506/100‑8 (MQ=255)
              cGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtt                                                                                    >  1:3768291/1‑100 (MQ=255)
              cGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGt                                                                                     >  1:3182437/1‑99 (MQ=255)
              cGCCTGTGGCTGTTTGCGTAATGCCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtt                                                                                    <  1:1889211/100‑1 (MQ=255)
              cGCCTGGGGCTGTTAGCGTAATGCCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACACGCTTACACTTGCCGAATACGGCAATGTCg                                                                                      >  2:1023609/1‑98 (MQ=255)
                  tgtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTa                                                                                <  2:1408064/100‑1 (MQ=255)
                   gtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAc                                                                               <  2:1814157/100‑1 (MQ=255)
                    tgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAc                                                                               <  2:1028490/99‑1 (MQ=255)
                    tgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACACCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTAcc                                                                              >  2:295644/1‑100 (MQ=255)
                     ggCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCt                                                                             >  1:2985837/1‑100 (MQ=255)
                      gCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTc                                                                            >  2:3295418/1‑100 (MQ=255)
                       cTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcg                                                                           <  2:69245/100‑1 (MQ=255)
                       cTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcg                                                                           <  1:69040/100‑1 (MQ=255)
                        tGTTTGCGTAATGTCATGCGTAAACAGTGGGCGGAAGCGATTTGGAGCCGAGCAAACTCCATTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcgc                                                                          >  2:1995190/1‑100 (MQ=255)
                         gTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcgcg                                                                         >  1:2237886/1‑100 (MQ=255)
                            tGCGTAATGTCATGCGTGAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAg                                                                      >  2:2966545/1‑100 (MQ=255)
                            tGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAg                                                                      <  1:1247039/100‑1 (MQ=255)
                             gCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGt                                                                     >  1:274004/1‑100 (MQ=255)
                             gCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGGCGTTTTTACCTCGCGCAGt                                                                     >  2:4065061/1‑100 (MQ=255)
                             gCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACACCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGt                                                                     >  2:428781/1‑100 (MQ=255)
                             gCGTAATGTCATGCGTAAACAGTGGGCGAAACCGATTGGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGt                                                                     <  2:2795726/100‑1 (MQ=255)
                               gtagtggtaTGCGTGAACAGTGGGTGAAAGCGAGTTGCAGCCGTGCAGACTCCCTTACGCTGGGCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTa                                                                    <  2:1796543/91‑1 (MQ=255)
                                tAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAAc                                                                  >  1:288440/1‑100 (MQ=255)
                                 aaTGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACg                                                                 <  2:1303191/100‑1 (MQ=255)
                                          cGTAAACAGTGGGCGAAAGCGATTTGCAGCCGAGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgg                                                        <  2:1450175/100‑1 (MQ=255)
                                           gTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGc                                                       <  1:3295438/100‑1 (MQ=255)
                                             aaaCAGTGGGCGAAAGCGATTTGCAGCCGGGCAACCTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGc                                                       <  1:262366/98‑1 (MQ=255)
                                                cAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatca                                                  <  1:1364914/100‑1 (MQ=255)
                                                cAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatca                                                  <  2:2802592/100‑1 (MQ=255)
                                                cAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatca                                                  <  1:2937230/100‑1 (MQ=255)
                                                aaGTGGGCGGAAGAGATTTGCAGCCGGGCCAACTTACTTACACATGCACACTCCACCAAGGCCGTTTACACCTCGCGAAGTAAAGACCggcgggca                                                      >  2:2889932/2‑89 (MQ=255)
                                                 aGTGGGCGAAAGCTATTTGCAGCCGGGCAAACTCCCTTGCACTTGCAGAATCCGGCAATGTCGTTTTAACCTCGCGCAGTGACGACAGTGAGGCAatcat                                                 >  2:2546710/1‑100 (MQ=255)
                                                 aGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatcat                                                 <  2:1286724/100‑1 (MQ=255)
                                                  gTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatcatc                                                >  1:1576172/1‑100 (MQ=255)
                                                  gTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGCCAatcatc                                                >  1:4093319/1‑100 (MQ=255)
                                                   tGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAAtcatca                                               <  1:2534735/100‑1 (MQ=255)
                                                     ggCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAAc                                             <  1:1738157/100‑1 (MQ=255)
                                                      gCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACt                                            >  2:1522725/1‑100 (MQ=255)
                                                       cGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTg                                           >  1:3863989/1‑100 (MQ=255)
                                                          aaGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGccc                                        >  2:2172940/1‑100 (MQ=255)
                                                            gCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGccccc                                      <  2:1057602/100‑1 (MQ=255)
                                                              gATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCaa                                    >  2:2783672/1‑100 (MQ=255)
                                                               aattGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCaaa                                   >  1:3826007/3‑100 (MQ=255)
                                                                  tGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtt                                >  1:1366580/1‑100 (MQ=255)
                                                                   gCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGttt                               >  1:1629497/1‑100 (MQ=255)
                                                                   gCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGttt                               >  1:2906220/1‑100 (MQ=255)
                                                                    cAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtttt                              <  2:65242/100‑1 (MQ=255)
                                                                    cAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGttt                               >  2:3013730/1‑99 (MQ=255)
                                                                      gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                            <  1:1198661/100‑1 (MQ=255)
                                                                      gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                            >  2:3592071/1‑100 (MQ=255)
                                                                      gCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                            <  1:3522530/100‑1 (MQ=255)
                                                                         gggCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTcg                            >  2:3007429/1‑97 (MQ=255)
                                                                          ggCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCt                        <  1:3769627/100‑1 (MQ=255)
                                                                           gCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCtt                       >  1:2770901/1‑100 (MQ=255)
                                                                              aaCTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCaa                    >  1:943251/1‑100 (MQ=255)
                                                                              aaCTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCaa                    <  2:1868240/100‑1 (MQ=255)
                                                                              aaCTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCaa                    >  1:243267/1‑100 (MQ=255)
                                                                                cTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGTAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACg                  >  1:2922977/1‑100 (MQ=255)
                                                                                cTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACg                  >  2:3023128/1‑100 (MQ=255)
                                                                                  cccTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTa                >  1:2466674/1‑100 (MQ=255)
                                                                                     acgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACg             >  1:1654750/7‑100 (MQ=255)
                                                                                      tACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGACAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACgg            <  2:4093292/100‑1 (MQ=255)
                                                                                      cgacattGCCGGATTCGGCAATGTCGTTTTTTTTTTTTTTAGTAACGACCATAAGGCAATCATAAAATACCTCCAAAGTTTTAGCTCTTtaacgtaacg             >  2:1892039/6‑99 (MQ=255)
                                                                                      cgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACgg            >  2:2420512/6‑100 (MQ=255)
                                                                                      cgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACgg            >  1:9905/6‑100 (MQ=255)
                                                                                      cgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACgg            >  1:1232897/6‑100 (MQ=255)
                                                                                      cgacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCAAGAGGCAGTCATAAAATGCCCCCAAAGTATTCGCGCGTCAGAGTAAcag            >  2:1094580/6‑98 (MQ=255)
                                                                                       acacTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACggg           <  1:2615750/100‑1 (MQ=255)
                                                                                          cTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGt        <  1:2551308/100‑1 (MQ=255)
                                                                                          cTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGt        >  1:288341/1‑100 (MQ=255)
                                                                                           ttGCCGAATCCGGCACTGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATGAACTGCCCCCAAAGTTTTCGCGCTTCAACGTCACGGGCGt        >  2:926756/1‑99 (MQ=255)
                                                                                           ttGCCGAATCCGGCAATGTCGTTCTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTCACGGGc          >  2:3569205/1‑97 (MQ=255)
                                                                                            tGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAc      >  1:30697/1‑100 (MQ=255)
                                                                                            tGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAc      >  1:364207/1‑100 (MQ=255)
                                                                                            tGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATAAACTGCCCCCAAAGTTTACGCGCTTCAACGTAACGGGCGTAc      >  2:3515741/1‑100 (MQ=255)
                                                                                               cGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCg   >  2:1926363/1‑100 (MQ=255)
                                                                                                gAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCgg  >  2:3008801/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
TGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGG  >  NC_000913/1608607‑1608802

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: