Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 860,154 A→T 5.7% intergenic (‑126/+20) 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,1540AT 5.7% 376.2 / 10.8 141intergenic (‑126/+20)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 (79/50);  new base T (4/4);  total (86/55)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-01

Only 100 of 151 total aligned reads displayed.

GTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCC  >  NC_000913/860006‑860303
                                                                                                                                                    |                                                                                                                                                     
 tgACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCgacgacc                                                                                                                                                     >  2:938657/1‑149 (MQ=255)
  gACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTTTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGAcgc                                                                                                                                                    >  2:1087817/1‑150 (MQ=255)
    cGATAACTTTTACGCTCATAGGTGTTCTCATGGCTCGAACATGAAACCGTAACAGAGTAATAACAATGTGACGCAGAGAACAAAGTATATGTCGAATCAAAGTAAGAATGAAATTGATGATGTGAATGATATAGCCCGGAGACGAcgccg                                                                                                                                                  >  1:119811/1‑150 (MQ=255)
         aCTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCggg                                                                                                                                             <  1:1349578/150‑1 (MQ=255)
         aCTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCggg                                                                                                                                             <  1:1620853/150‑1 (MQ=255)
          cTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                                                                            <  2:441649/150‑1 (MQ=255)
          cTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGc                                                                                                                                            >  1:1295051/1‑150 (MQ=255)
          cTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGc                                                                                                                                            >  2:1509018/1‑150 (MQ=255)
          cTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGAAGACGCCGCCGGGc                                                                                                                                            >  2:1157363/1‑150 (MQ=255)
           ttttACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGcc                                                                                                                                           >  1:1345121/1‑150 (MQ=255)
              tACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaa                                                                                                                                        <  2:6831/150‑5 (MQ=255)
              tACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgag                                                                                                                                        <  2:513481/150‑1 (MQ=255)
              tACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgag                                                                                                                                        <  2:780833/150‑1 (MQ=255)
                cGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCCCAAATTATATTTCGAATGAAAGTAAGGAGGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgagga                                                                                                                                      >  1:1296388/1‑150 (MQ=255)
                cGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaatc                                                                                                                                      <  1:534690/150‑7 (MQ=255)
                cGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAAAGATTTAGCCCGGCGACGACGCCGCCGGGCCgagga                                                                                                                                      >  2:1696995/1‑150 (MQ=255)
                  cTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGgaga                                                                                                                                    >  1:681825/1‑150 (MQ=255)
                    cATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGAtt                                                                                                                                  >  2:637272/1‑150 (MQ=255)
                      tAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTAc                                                                                                                                >  2:1636010/1‑150 (MQ=255)
                          tgtgCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCa                                                                                                                            >  1:272259/1‑150 (MQ=255)
                              cTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATg                                                                                                                        >  2:366159/1‑150 (MQ=255)
                                 cTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGtta                                                                                                                     >  1:1760401/1‑149 (MQ=255)
                                 cTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTc                                                                                                                     <  2:1481680/150‑1 (MQ=255)
                                 cTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTc                                                                                                                     >  1:1664579/1‑150 (MQ=255)
                                 cTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTc                                                                                                                     <  1:150125/150‑1 (MQ=255)
                                  tGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCa                                                                                                                    <  1:1745828/150‑1 (MQ=255)
                                    gCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGt                                                                                                                  >  2:232234/1‑150 (MQ=255)
                                     cTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTa                                                                                                                 >  2:1180418/1‑150 (MQ=255)
                                               aaaCCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatga                                                                                                       >  1:913205/1‑150 (MQ=255)
                                                aaCCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAATTTGGTGATGTGATTGAGTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGCGTACAGGATATGTTCAGGACGGGCgatgat                                                                                                      >  2:333778/1‑150 (MQ=255)
                                                       aCAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTcc                                                                                               <  2:1451733/150‑1 (MQ=255)
                                                        cAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCt                                                                                              <  2:720956/150‑1 (MQ=255)
                                                        cAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCt                                                                                              >  2:357975/1‑150 (MQ=255)
                                                          gtgtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGCCGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATAGTACTGa                                                                                            >  2:1400245/1‑150 (MQ=255)
                                                            gtAATAAAAATGTGACGCAGAGCACAAATTATATTTCCAATGAAAGTAAGGAGGAAATTCATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGAAAACAGCATATCTTCAGTACGGGCGATGATATCGTCCAGAGc                                                                                          >  2:188519/1‑150 (MQ=255)
                                                             taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCg                                                                                         >  1:360915/1‑150 (MQ=255)
                                                             taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACAGGCGATGATATCGTCCTGAGCg                                                                                         >  2:1304858/1‑150 (MQ=255)
                                                               ataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGATCAGTACGGGCGATGATATCGTCCTGAGCGta                                                                                       >  1:1798228/1‑149 (MQ=255)
                                                                taaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCt                                                                                      >  1:1124679/1‑150 (MQ=255)
                                                                 aaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTg                                                                                     <  2:1035510/150‑1 (MQ=255)
                                                                   cAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGga                                                                                   <  1:1986198/150‑1 (MQ=255)
                                                                    attgtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgag                                                                                  <  1:776689/148‑1 (MQ=255)
                                                                         gACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAgc                                                                             >  2:1679431/1‑150 (MQ=255)
                                                                           cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAgcgc                                                                           <  1:623971/150‑1 (MQ=255)
                                                                              agagCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGTGGt                                                                        <  2:1421717/150‑1 (MQ=255)
                                                                               gagCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTg                                                                       >  1:998786/1‑150 (MQ=255)
                                                                                agCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTga                                                                      >  1:1109003/1‑150 (MQ=255)
                                                                                 gCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGCTGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCTGTgaa                                                                     <  2:1346655/150‑1 (MQ=255)
                                                                                   acaAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACTACGACGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaaga                                                                   >  1:1790107/1‑150 (MQ=255)
                                                                                   acaAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTgaaga                                                                   <  2:525360/150‑1 (MQ=255)
                                                                                     aaaTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAc                                                                 >  1:1963139/1‑150 (MQ=255)
                                                                                         tataTTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg                                                             <  1:1338437/150‑1 (MQ=255)
                                                                                          atatTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGa                                                            >  2:2012387/1‑150 (MQ=255)
                                                                                            atTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAAt                                                          >  1:1680327/1‑150 (MQ=255)
                                                                                               tCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGc                                                       >  2:1421034/1‑150 (MQ=255)
                                                                                               tCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATAACAGCATATTTTAAGTACGGGCGATGATATCGTCCTGAGCGTCTGGACAGAGCGCGGTGAAGAACGCGGAATAGc                                                       >  1:764960/1‑150 (MQ=255)
                                                                                                 gAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGccc                                                     <  2:1760401/150‑1 (MQ=255)
                                                                                                 gAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGccc                                                     >  1:606214/1‑150 (MQ=255)
                                                                                                 gAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGccc                                                     >  2:642761/1‑150 (MQ=255)
                                                                                                   aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCcgc                                                   >  1:763707/1‑150 (MQ=255)
                                                                                                   aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCcgc                                                   >  1:251676/1‑150 (MQ=255)
                                                                                                    tGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCcgcg                                                  >  1:696675/1‑150 (MQ=255)
                                                                                                    tGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCcgcg                                                  >  1:445835/1‑150 (MQ=255)
                                                                                                    tGAAAGTAAGGAGGAAATTGATGATCTGAATGATTTAGCCCGGCGACGCCGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCcgcg                                                  <  2:386327/150‑1 (MQ=255)
                                                                                                     gAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGa                                                 >  2:665522/1‑150 (MQ=255)
                                                                                                      aaaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGac                                                >  2:1550602/1‑150 (MQ=255)
                                                                                                         gTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACg                                             <  2:1267235/150‑1 (MQ=255)
                                                                                                            aGGATGAAATTGATGATGGGAACGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCAc                                          <  1:1037290/150‑1 (MQ=255)
                                                                                                             ggATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACt                                         >  1:1252735/1‑150 (MQ=255)
                                                                                                               aTGAAATTGATGATGTGAATGATTTAGCCCGGCGACGCCGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGCCACGCCCTAc                                       >  1:277258/1‑150 (MQ=255)
                                                                                                                 gAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCa                                     <  1:539619/150‑1 (MQ=255)
                                                                                                                 gAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTAAGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCa                                     >  1:1393036/1‑150 (MQ=255)
                                                                                                                   aaTTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTGAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGa                                   >  1:1426747/1‑150 (MQ=255)
                                                                                                                   aaTTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGa                                   >  1:948794/1‑150 (MQ=255)
                                                                                                                       gatgatGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATcgc                               <  2:1688538/150‑1 (MQ=255)
                                                                                                                       gatgatGGGCATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATcgc                               <  2:310080/150‑1 (MQ=255)
                                                                                                                          gatgTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGAt                            <  1:1027740/150‑1 (MQ=255)
                                                                                                                           atgTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATa                           <  2:769748/150‑1 (MQ=255)
                                                                                                                           atgTGAATGATTTAGCCAGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATAGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCAGCGACACGCACTACCAGATCGCGATa                           <  2:552622/150‑1 (MQ=255)
                                                                                                                            tgtgAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATAc                          <  2:426051/150‑1 (MQ=255)
                                                                                                                            tgtgAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATAc                          >  2:624349/1‑150 (MQ=255)
                                                                                                                             gtgAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACt                         >  2:1877127/1‑150 (MQ=255)
                                                                                                                                  tGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCAGCGACACGCACTACCAGATCGCGATACTGATCg                    >  2:572217/1‑150 (MQ=255)
                                                                                                                                   gATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCgg                   <  1:1594698/150‑1 (MQ=255)
                                                                                                                                   gATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGAGCgg                   >  1:786126/1‑150 (MQ=255)
                                                                                                                                    aTTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCggg                  <  1:1873168/150‑1 (MQ=255)
                                                                                                                                      ttAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGAt                <  1:1823174/150‑1 (MQ=255)
                                                                                                                                       cggCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATg               >  2:1609398/3‑150 (MQ=255)
                                                                                                                                        aGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGt              >  1:58506/1‑150 (MQ=255)
                                                                                                                                         gCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGtt             >  1:570519/1‑150 (MQ=255)
                                                                                                                                           ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGtttt           >  1:1968897/1‑150 (MQ=255)
                                                                                                                                           ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGtttt           >  1:1545803/1‑150 (MQ=255)
                                                                                                                                           ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGtttt           >  1:1567267/1‑150 (MQ=255)
                                                                                                                                              gCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTc        <  1:357975/150‑1 (MQ=255)
                                                                                                                                               cgacgaCGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCg       <  1:143316/150‑1 (MQ=255)
                                                                                                                                               cgacgaCGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCg       <  1:445455/150‑1 (MQ=255)
                                                                                                                                               cgacgaCGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCg       <  1:1906316/150‑1 (MQ=255)
                                                                                                                                               cgacgaCGCCGCCGGGCCGAGGAGATTACAGCATATGTTAAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGTAATATAGAGCGACAGGCACTACCAGATCGAGATACTGATCTGGATGTTTTTTCg       <  2:272259/150‑1 (MQ=255)
                                                                                                                                                  cgacgCCGCCGGGACGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCATGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCAt    >  2:1398753/1‑150 (MQ=255)
                                                                                                                                                    acgCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATcc  <  2:1221291/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCC  >  NC_000913/860006‑860303

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: