Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 860,151 A→G 19.2% intergenic (‑123/+23) ybiV ← / ← ybiW sugar phosphatase; fructose‑1‑P/ribose‑5‑P/glucose‑6‑P phosphatase/putative pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913860,1510AG19.2% 349.9 / 74.7 171intergenic (‑123/+23)ybiV/ybiWsugar phosphatase; fructose‑1‑P/ribose‑5‑P/glucose‑6‑P phosphatase/putative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base A (64/74);  new base G (15/18);  total (79/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 182 total aligned reads displayed.

GAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAAT  >  NC_000913/860052‑860247
                                                                                                   |                                                                                                
  aaCCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGcgacg                                                                                                <  1:1876752/100‑1 (MQ=255)
  aaCCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGcgacg                                                                                                >  2:4001424/1‑100 (MQ=255)
      gTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGacgacg                                                                                             >  1:3997164/1‑99 (MQ=255)
          cAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACgccgcc                                                                                        >  1:4364499/1‑100 (MQ=255)
           aGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGccgccg                                                                                       <  2:3029888/100‑1 (MQ=255)
           aGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGccgccg                                                                                       <  1:4404099/100‑1 (MQ=255)
             tgtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCggg                                                                                     <  1:1335340/100‑1 (MQ=255)
             tgtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCggg                                                                                     <  2:356976/100‑1 (MQ=255)
             tgtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCggg                                                                                     <  2:261903/100‑1 (MQ=255)
              gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                    <  1:4210949/100‑1 (MQ=255)
              gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                    <  2:398060/100‑1 (MQ=255)
              gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                    <  2:1254110/100‑1 (MQ=255)
              gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                    <  1:313617/100‑1 (MQ=255)
              gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGc                                                                                    >  2:1822012/1‑100 (MQ=255)
              gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGCGAATGATTTAGCCCGGCGACGACGCCGCCGGGc                                                                                    <  2:2321909/100‑1 (MQ=255)
              gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGACGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                    <  1:2717754/100‑1 (MQ=255)
              gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATGGATGAGGTGAATGATTTAGCCCGGCGGCGTCGTCGCCggg                                                                                     <  1:3256072/99‑1 (MQ=255)
               taataataaTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                    <  1:567539/99‑1 (MQ=255)
               taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGct                                                                                   <  1:2090952/100‑2 (MQ=255)
               taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGcc                                                                                   <  2:4021223/100‑1 (MQ=255)
               taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGcc                                                                                   >  2:3699592/1‑100 (MQ=255)
                 ataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAAATGATGATGTGAATGATTTAGACGGGCGACGACGCCGCCTGGCCga                                                                                 >  1:1617608/1‑100 (MQ=255)
                  taaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAACTAAGGATGAAATTGATGATGTGAATTATTTAGCCCGGCGACGACGCCGCCGGGCCgag                                                                                <  2:502395/100‑1 (MQ=255)
                   aaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgagg                                                                               <  1:3316603/100‑1 (MQ=255)
                   aaCAATGTGACGCAGAGCACAAATTATATTGCGAATGAAAGTAAGGATGTAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGct                                                                                   <  1:3372232/96‑2 (MQ=255)
                    aCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgagga                                                                              >  2:2284501/1‑100 (MQ=255)
                      aaTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGgaga                                                                            >  2:2813631/1‑100 (MQ=255)
                       aTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAAGTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGAt                                                                           >  2:2624758/1‑100 (MQ=255)
                        tgtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGAtt                                                                          <  1:1428132/100‑1 (MQ=255)
                          tgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTAc                                                                        <  1:3242353/100‑1 (MQ=255)
                          tgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTAc                                                                        >  1:383028/1‑100 (MQ=255)
                           gACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACa                                                                       >  2:1038308/1‑100 (MQ=255)
                             cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGc                                                                     <  1:4393239/100‑1 (MQ=255)
                             cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGc                                                                     <  2:3422772/100‑1 (MQ=255)
                             cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGc                                                                     <  2:838520/100‑1 (MQ=255)
                             cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGc                                                                     >  1:1761023/1‑100 (MQ=255)
                                agagCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCata                                                                  <  2:2740267/100‑1 (MQ=255)
                                  agCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATg                                                                >  2:611600/1‑100 (MQ=255)
                                  agCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATg                                                                <  2:3321973/100‑1 (MQ=255)
                                  agCACAAATTATATTTCGAATGAAAGGAAGGATGAAATTGATGCTGTGAATGATTTAGCCCGGCGACGGCGGCGCCGGGCTGAGAAGATTACAGCATATg                                                                <  1:597709/100‑1 (MQ=255)
                                     acaAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTc                                                             <  1:1916649/100‑1 (MQ=255)
                                      caAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCa                                                            >  1:3614838/1‑100 (MQ=255)
                                       aaaTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAg                                                           >  2:1632811/1‑100 (MQ=255)
                                         aTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTa                                                         >  2:2492484/1‑100 (MQ=255)
                                             tatTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACggg                                                     <  2:447232/100‑1 (MQ=255)
                                               tttCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCg                                                   >  2:3948295/1‑100 (MQ=255)
                                               tttCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCg                                                   >  2:2245921/1‑100 (MQ=255)
                                                 tCGGATGAAAGTAAGGACGAAATTTATGATGTGAATGATTCAGCCGGGCCACGATGCTGCCGGGCCTAGGAGAGTACAGCATATGGTCAGTACGGGCgat                                                 <  1:155785/100‑1 (MQ=255)
                                                 tCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgat                                                 >  2:1453687/1‑100 (MQ=255)
                                                  cGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatg                                                <  2:3853325/100‑1 (MQ=255)
                                                  cGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatg                                                >  2:596383/1‑100 (MQ=255)
                                                  cGAATGAAAGTAAGGATGAAATTGATGAAGTGAATGATTTAGCCCGGTGACGACGCCGCAGGGCCGAGGAGATTACACCATATGTTCAGTACGGGCgat                                                 >  1:3929608/1‑99 (MQ=255)
                                                   gAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatga                                               <  2:3286347/100‑1 (MQ=255)
                                                      tGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGatat                                            <  1:4004230/100‑1 (MQ=255)
                                                      tGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGatat                                            >  2:1645689/1‑100 (MQ=255)
                                                         aaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGt                                         <  1:2196567/100‑1 (MQ=255)
                                                          aGTAATTATGAAATAGATGATGTTAATGATGTAGCATTGCGACGACGCCGCCTGGCCGAGGATATTATATCATATGTTCAGTACGGGCGATGATATCGTc                                        <  1:668638/100‑1 (MQ=255)
                                                          aGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGt                                         >  1:2520380/1‑99 (MQ=255)
                                                           gTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTcc                                       <  2:2348722/100‑1 (MQ=255)
                                                              aGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGa                                    <  2:587962/100‑1 (MQ=255)
                                                               ggATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAg                                   >  1:1987513/1‑100 (MQ=255)
                                                                 aTGAAATTGATGATGTGAATGATTTAGCCCGGCGACGGCGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCg                                 >  1:2352123/1‑100 (MQ=255)
                                                                 aTGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCg                                 <  2:1387618/100‑1 (MQ=255)
                                                                  tGAAATTGGTGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCAGATGTTCAGTACGCTCGATGATATCGTCCTg                                     >  1:260730/1‑95 (MQ=255)
                                                                   gAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTc                               <  2:4348613/100‑1 (MQ=255)
                                                                    aaaTGGATGATGTGAATGATGTCGCCCGGCTACGACGCCGCCGGGCCGAGTAGATTACAGCATATGTTCATTACGGGCGATTATATCGTCCTGAGCGTCt                              <  1:914820/100‑1 (MQ=255)
                                                                     aaTTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTc                               <  2:3048753/98‑1 (MQ=255)
                                                                        tgatgaTGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgag                          <  1:2196895/100‑1 (MQ=255)
                                                                        tgatgaTGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgag                          <  2:157512/100‑1 (MQ=255)
                                                                          atgatgTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgagag                        >  2:2926248/1‑100 (MQ=255)
                                                                             aggtgAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAgc                     <  1:2918626/98‑1 (MQ=255)
                                                                                tgAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGcgcg                  >  2:3148688/1‑100 (MQ=255)
                                                                                 gAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCgg                 <  1:2420099/100‑1 (MQ=255)
                                                                                  aaTGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGt                <  1:1011588/100‑1 (MQ=255)
                                                                                   aTGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTg               <  1:411827/100‑1 (MQ=255)
                                                                                    tGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTga              <  2:622342/100‑1 (MQ=255)
                                                                                    tGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTga              <  1:1412358/100‑1 (MQ=255)
                                                                                     gATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaa             <  2:1790466/100‑1 (MQ=255)
                                                                                      aTTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaag            <  1:1058713/100‑1 (MQ=255)
                                                                                      aTTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaag            <  1:1405662/100‑1 (MQ=255)
                                                                                      aTTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaag            <  1:1476873/100‑1 (MQ=255)
                                                                                       tttAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaaga           >  1:2317124/1‑100 (MQ=255)
                                                                                       tttAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaaga           >  1:1684756/1‑100 (MQ=255)
                                                                                       ntcggCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaaga           >  1:3802226/5‑100 (MQ=255)
                                                                                         cggCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAc         >  2:780328/3‑100 (MQ=255)
                                                                                         cggCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAc         >  2:278177/3‑100 (MQ=255)
                                                                                          ggCCCTGCGGCGTCGTCGCTGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGTGCGCGGTgaa             >  1:862306/2‑95 (MQ=255)
                                                                                          ggCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcg        >  2:3192717/2‑100 (MQ=255)
                                                                                          aGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcg        >  2:3266969/1‑100 (MQ=255)
                                                                                           gCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgc       >  2:175290/1‑100 (MQ=255)
                                                                                           gCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgc       >  2:3643303/1‑100 (MQ=255)
                                                                                           gCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCAGATGTTCAGGACGGGCGATGATACCGTGCTGGGCGTCTGGAGAGAGCGCGGCTAAGAAcgc       >  1:3357447/1‑100 (MQ=255)
                                                                                            cccGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGCGCGTCTGGAGAGAGCGCGGTGAAGACcgcg      >  1:3033406/1‑100 (MQ=255)
                                                                                            cccGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgcg      >  2:3594190/1‑100 (MQ=255)
                                                                                            cccGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgcg      >  2:2634459/1‑100 (MQ=255)
                                                                                            cccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgcg      >  2:2529441/1‑100 (MQ=255)
                                                                                             ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg     <  1:2133050/100‑1 (MQ=255)
                                                                                              cGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGa    <  1:3485573/100‑1 (MQ=255)
                                                                                              cGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGa    <  1:3518228/100‑1 (MQ=255)
                                                                                                gCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAAt  <  2:3933535/100‑1 (MQ=255)
                                                                                                   |                                                                                                
GAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAAT  >  NC_000913/860052‑860247

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: