Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 314,255 T→C 15.1% intergenic (‑450/‑102) ykgP ← / → eaeH pseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913314,2550TC15.1% 470.1 / 68.0 199intergenic (‑450/‑102)ykgP/eaeHpseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified
Reads supporting (aligned to +/- strand):  ref base T (87/82);  new base C (12/18);  total (99/100)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 202 total aligned reads displayed.

TAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCT  >  NC_000913/314156‑314354
                                                                                                   |                                                                                                   
taaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCt                                                                                                     >  1:1247866/1‑100 (MQ=255)
  aaaaTATAATCCCGCGCCCAACCACCTGATGATTGGCTATAGACACTGGATATACTAGGTGGCGGTGCACTTTCTTACAGAAAGGTATTTCCTTTTCTGc                                                                                                   >  1:2028378/1‑100 (MQ=255)
  aaaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGc                                                                                                   >  1:4335418/1‑100 (MQ=255)
   aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCg                                                                                                  >  1:2886620/1‑100 (MQ=255)
   aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCg                                                                                                  >  2:3218827/1‑100 (MQ=255)
   aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCg                                                                                                  <  2:313944/100‑1 (MQ=255)
     atatAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGa                                                                                                <  2:2385488/100‑1 (MQ=255)
       ataATCCCGCGCCCAGCCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAATGTATTTCCTTTTCCGCAGaaa                                                                                              <  1:475344/100‑1 (MQ=255)
        taATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                             <  1:87703/100‑1 (MQ=255)
        taATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                             >  2:3266341/1‑100 (MQ=255)
        taATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                             >  2:4325030/1‑100 (MQ=255)
        taATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                             >  2:2585269/1‑100 (MQ=255)
         aaTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAg                                                                                            <  1:3237862/100‑1 (MQ=255)
         aaTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAg                                                                                            <  2:301211/100‑1 (MQ=255)
          aTCCCGCGCCCAACCACCTGATGAGTTGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAgg                                                                                           <  1:1810257/100‑1 (MQ=255)
          aTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAgg                                                                                           <  2:1138548/100‑1 (MQ=255)
           tCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCCCTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGa                                                                                          >  1:3700332/1‑100 (MQ=255)
           tCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGa                                                                                          <  1:1470378/100‑1 (MQ=255)
           tCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGa                                                                                          <  2:2201260/100‑1 (MQ=255)
            cccGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGaa                                                                                         >  2:2957705/1‑100 (MQ=255)
            cccGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGaa                                                                                         >  1:2707749/1‑100 (MQ=255)
            cccGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGaa                                                                                         <  1:381853/100‑1 (MQ=255)
             ccGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGaaa                                                                                        >  1:1903945/1‑100 (MQ=255)
              cgcgCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGAAAt                                                                                       <  1:1919164/100‑1 (MQ=255)
               gcgcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGAAAta                                                                                      <  2:3660665/100‑2 (MQ=255)
               gcgcCCAACCACCTGATGAGTGGCAATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGAAAta                                                                                      <  1:3371559/100‑2 (MQ=39)
                cgcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGAAAtac                                                                                     <  2:15797/100‑3 (MQ=255)
                 gaccAACGACCTGATGAGTGGCTATATTCACTGGAGATATTAAGTCGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGAAAta                                                                                      <  1:1643093/96‑2 (MQ=37)
                  cccAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCggg                                                                                   >  2:1737465/1‑100 (MQ=255)
                  cccAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCggg                                                                                   >  2:3327418/1‑100 (MQ=255)
                  cccAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCggg                                                                                   >  2:361111/1‑100 (MQ=255)
                    cAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGTTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGaa                                                                                 <  2:1498952/100‑1 (MQ=255)
                      accaccTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAAt                                                                               >  2:920230/1‑100 (MQ=255)
                       ccaccTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATc                                                                              >  2:1179534/1‑100 (MQ=255)
                        caccTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATcc                                                                             >  1:2494045/1‑100 (MQ=255)
                         accTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATccc                                                                            >  2:1893158/1‑100 (MQ=255)
                         accTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATccc                                                                            >  1:4414326/1‑100 (MQ=255)
                          ccTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATcccc                                                                           >  2:3686714/1‑100 (MQ=255)
                           cTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCg                                                                          <  2:1363064/100‑1 (MQ=255)
                                gaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGATATCCCCGGtttt                                                                     <  2:2902941/100‑1 (MQ=255)
                                  gTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTc                                                                   <  1:4328746/100‑1 (MQ=255)
                                       tataGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACa                                                              <  2:1755736/100‑1 (MQ=255)
                                       tataGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACa                                                              >  2:3822513/1‑100 (MQ=255)
                                        ataGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACaa                                                             <  2:2427118/100‑1 (MQ=255)
                                         taGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAg                                                            <  2:261768/100‑1 (MQ=255)
                                          aGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGc                                                           >  2:32881/1‑100 (MQ=255)
                                              aCTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGAc                                                       <  2:3681525/100‑1 (MQ=255)
                                                 ggATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGcc                                                    <  2:1843388/100‑1 (MQ=255)
                                                 ggATATATTAGGTGGCGGTGCACTGTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGcc                                                    >  1:1769382/1‑100 (MQ=255)
                                                  gATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCa                                                   >  1:1313323/1‑100 (MQ=255)
                                                  gATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCa                                                   <  2:2687245/100‑1 (MQ=255)
                                                   ataaatTATGTGGCTTTGCACTTTATTACATAAAGGTAATTCCTTTTCTGCGGAAAAGGAAATCGGCAAATCCCCGGTTTTTCTACCAAGCAGACGCCat                                                  <  1:2052875/100‑1 (MQ=255)
                                                       atTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCattatt                                              >  2:3982983/1‑100 (MQ=255)
                                                         tAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtg                                            <  1:705900/100‑1 (MQ=255)
                                                         tAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtg                                            <  1:619642/100‑1 (MQ=255)
                                                            gTGGCGGTGCACTTTCTTACATAAAGGTATTTACTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTgtgt                                         >  2:4185356/1‑100 (MQ=255)
                                                             tGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTc                                        >  2:2511229/1‑100 (MQ=255)
                                                               gCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTg                                      <  1:2158051/100‑1 (MQ=255)
                                                                cGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGc                                     >  1:2297247/1‑100 (MQ=255)
                                                                 ggTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGcc                                    <  2:2729560/100‑1 (MQ=255)
                                                                    gCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTAt                                 >  1:4289586/1‑100 (MQ=255)
                                                                      aCTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGt                               >  2:1636889/1‑100 (MQ=255)
                                                                      aCTTTATTACATAAAGGTGTTTCCTTTTCTGCGGAAAAGGAAATCCGGAAATCCCCGGGTTATCTGATACAGAAACAACATTATTTACGTCTGcat                                   >  1:4059104/1‑94 (MQ=255)
                                                                        tttCTTACATAAAGGTATTTCCTTTTCTGCGGGAAATGAAATCGGGAAATCCCCTGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTc                             <  2:3094782/100‑1 (MQ=255)
                                                                         ttCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCg                            <  2:912578/100‑1 (MQ=255)
                                                                         ttCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCg                            <  1:4246993/100‑1 (MQ=255)
                                                                         ttCTTACATAAAGGTATTTCCTTTTATGAGGAAAAGAAAATATAAAAATAACAGATTTTTTTGACAAGCAAATACCATTATTTTTGTCTGCATATAtt                              >  1:358430/1‑98 (MQ=255)
                                                                          tCTTACATGAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGt                           <  1:1319587/100‑1 (MQ=255)
                                                                           cTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGtt                          <  2:4177131/100‑1 (MQ=255)
                                                                           cTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGtt                          <  1:2407482/100‑1 (MQ=255)
                                                                            ttACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTa                         <  2:4414358/100‑1 (MQ=255)
                                                                            ttACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTATGCCTATGTTCGTTa                         >  2:4339675/1‑100 (MQ=255)
                                                                             tACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATAAAGAAATCCCCGGTTTTTCTGACAAGCAGACGCCAATATTTGTGTCTGCCTactttccttc                         >  1:4382116/1‑90 (MQ=255)
                                                                              aCATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAAt                       >  2:331732/1‑100 (MQ=255)
                                                                               cATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAAtt                      >  2:4112666/1‑100 (MQ=255)
                                                                                aTAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCATACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTc                     <  1:1082766/100‑1 (MQ=255)
                                                                                aTAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTc                     <  2:142382/100‑1 (MQ=255)
                                                                                 tAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCg                    >  2:1149529/1‑100 (MQ=255)
                                                                                 tAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCg                    >  1:2755401/1‑100 (MQ=255)
                                                                                   ctccgaTTTCCTTTTCCGCAGAAAAGGAAATCTGTAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGtt                  <  1:4135558/95‑1 (MQ=38)
                                                                                    cccgaTTTCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtc                 >  1:3652270/5‑100 (MQ=38)
                                                                                     ggTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtca                <  2:2452130/100‑1 (MQ=255)
                                                                                     ccgaTTTCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTTTTTTTTTTTTACGCCAATATTTGTGTCTGCCTATTTTCGTTAATTCGTtca                >  1:1671956/4‑100 (MQ=16)
                                                                                     ccgaTTTCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtca                >  2:3390627/4‑100 (MQ=39)
                                                                                     ccgaTTTCCTTTTCCGCAGAAAAGGAAATCGGGAAATACCCGGTTTTTCTGACAAGCAGACGCCATGATTTGTGTCTGTCTATGTTCGTTAATTCggtag                >  2:4135589/4‑96 (MQ=37)
                                                                                      cgaTTTCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcat               >  2:1509007/3‑100 (MQ=255)
                                                                                       tATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatc              >  2:2930175/1‑100 (MQ=255)
                                                                                       tATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatc              >  1:2273734/1‑100 (MQ=255)
                                                                                       gaTTTCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatc              >  2:1035926/2‑100 (MQ=255)
                                                                                        aTTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatca             <  1:3930333/100‑1 (MQ=255)
                                                                                        aTTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatca             <  1:257505/100‑1 (MQ=255)
                                                                                        aTTTCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatca             >  2:2857208/1‑100 (MQ=255)
                                                                                         tttCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatca             >  1:3133095/1‑99 (MQ=255)
                                                                                            ccTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGaa         <  2:1550895/100‑1 (MQ=255)
                                                                                             cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGaaa        >  1:3455864/1‑100 (MQ=255)
                                                                                              ttttCTGCTGAAAAGTAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAg            <  1:950508/95‑1 (MQ=255)
                                                                                              ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAAt       <  2:3781580/100‑1 (MQ=255)
                                                                                               tttCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAAtt      >  2:619247/1‑100 (MQ=255)
                                                                                                  cTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtc   <  2:2273744/100‑1 (MQ=255)
                                                                                                   tGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTTTTTTTTTTTTTTGCAATTATTTGTGTATGCCTATTTTTGTTATTTCGTTTATCAGGAAAttgttt  >  1:1575951/1‑96 (MQ=255)
                                                                                                   |                                                                                                   
TAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCT  >  NC_000913/314156‑314354

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: