Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 860,154 A→T 20.5% 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,1540AT20.5% 190.0 / 45.1 94intergenic (‑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 (32/42);  new base T (7/12);  total (39/55)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

CCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGG  >  NC_000913/860056‑860244
                                                                                                  |                                                                                          
ccGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGCCgacgac                                                                                           >  1:1391277/1‑100 (MQ=255)
 cGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGacgacg                                                                                          <  2:2284445/100‑1 (MQ=255)
     aCAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGAcgccgc                                                                                      >  1:2011088/1‑100 (MQ=255)
        gtgtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATGGATGATGTGAATGATTTAGCCCGGCGTCGTcgccgggct                                                                                   <  1:1988769/100‑5 (MQ=255)
         tgtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCggg                                                                                  >  2:78936/1‑100 (MQ=255)
          gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                 <  2:1457543/100‑1 (MQ=255)
          gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                 <  2:1896707/100‑1 (MQ=255)
          gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGc                                                                                 <  2:2238140/100‑1 (MQ=255)
          gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCggg                                                                                  <  2:1511367/99‑1 (MQ=255)
          gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGc                                                                                 >  2:977838/1‑100 (MQ=255)
          gtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGc                                                                                 <  1:921907/100‑1 (MQ=255)
           taataaCAATTTTACTCAGATCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGcc                                                                                <  2:1119985/100‑1 (MQ=255)
           taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGct                                                                                <  2:1361499/100‑2 (MQ=255)
           taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGct                                                                                <  2:5412/100‑2 (MQ=255)
           taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGct                                                                                <  2:75257/100‑2 (MQ=255)
           taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGct                                                                                <  2:75324/100‑2 (MQ=255)
            aataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGcta                                                                               <  2:1772104/100‑3 (MQ=255)
            aataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCg                                                                               <  1:449181/100‑1 (MQ=255)
             ataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCga                                                                              <  2:2542804/100‑1 (MQ=255)
              taaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaa                                                                             <  1:2383926/100‑5 (MQ=38)
              taaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgag                                                                             <  1:1447973/100‑1 (MQ=255)
               aaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgagg                                                                            <  1:1716362/100‑1 (MQ=255)
                aCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaatc                                                                           <  2:1262111/100‑7 (MQ=38)
                aCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgagga                                                                           <  1:1443069/100‑1 (MQ=255)
                aCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgagga                                                                           <  1:60195/100‑1 (MQ=255)
                   aTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGgaga                                                                         <  2:2474462/99‑1 (MQ=255)
                    tgtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGAtt                                                                       <  1:2052394/100‑1 (MQ=255)
                    tgtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGAtt                                                                       >  2:2501438/1‑100 (MQ=255)
                     gtgACGCAGAGCACAAATTATATTTCGAGTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTa                                                                      <  2:782846/100‑1 (MQ=255)
                     gtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTa                                                                      >  2:1304309/1‑100 (MQ=255)
                      tgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTAc                                                                     >  1:833735/1‑100 (MQ=255)
                       gACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACa                                                                    >  1:2575750/1‑100 (MQ=255)
                       gACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACa                                                                    <  2:1994242/100‑1 (MQ=255)
                        aCGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAg                                                                   <  2:1265238/100‑1 (MQ=255)
                        aCGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAg                                                                   <  1:2501419/100‑1 (MQ=255)
                         cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGc                                                                  <  2:2313735/100‑1 (MQ=255)
                         cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGc                                                                  >  1:1688333/1‑100 (MQ=255)
                           cAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGAAGCCGCCGGGCCGAGGAGATTACAGCat                                                                >  2:1148865/1‑100 (MQ=255)
                               gCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGt                                                            >  2:494157/1‑100 (MQ=255)
                               gCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGt                                                            >  1:1507921/1‑100 (MQ=255)
                                cacaAATTATATTTCGAATGAAAGTAAGGGTGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGtt                                                           >  2:848263/1‑100 (MQ=255)
                                 acaAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGtt                                                           >  1:1754770/1‑99 (MQ=255)
                                 acaAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTc                                                          <  1:1881160/100‑1 (MQ=255)
                                  caAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCa                                                         >  1:372234/1‑100 (MQ=255)
                                  caAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCa                                                         >  2:45200/1‑100 (MQ=255)
                                    aaTTATATTTCGAATGAAAGTAAGGATGAAATGGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGt                                                       <  1:1523169/100‑1 (MQ=255)
                                     aTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTa                                                      <  1:472825/100‑1 (MQ=255)
                                      ttATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTAc                                                     <  1:670403/100‑1 (MQ=255)
                                       tataTTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACg                                                    >  2:914391/1‑100 (MQ=255)
                                           tttCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCg                                                <  1:2337381/100‑1 (MQ=255)
                                           tttCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCg                                                <  2:1072833/100‑1 (MQ=255)
                                              cGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatg                                             >  1:1419260/1‑100 (MQ=255)
                                               gAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGTAGATTACAGCATATGTTCAGTACGGGCgatg                                             <  1:1228801/99‑1 (MQ=255)
                                               gAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatga                                            <  1:65722/100‑1 (MQ=255)
                                                aaTGAAAGTAAGGATAAAATTGATGATGTCAATGATTTAGCCCGCCGACGACTCCGCCGGGCCGAGCAGATTACAGCATATGTTCAGTACGGGCgatg                                             <  1:1214120/98‑1 (MQ=255)
                                                 atcaaaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGata                                          <  2:1906029/97‑1 (MQ=255)
                                                 aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGata                                          <  2:1904460/100‑1 (MQ=255)
                                                 aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGata                                          >  2:997074/1‑100 (MQ=255)
                                                 aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGata                                          <  1:1143756/100‑1 (MQ=255)
                                                   gAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATc                                        >  2:659553/1‑100 (MQ=255)
                                                     aaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGt                                      <  2:776923/100‑1 (MQ=255)
                                                      aGTAAGGATGAAATTGATGAGGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTc                                     <  2:948746/100‑1 (MQ=255)
                                                         aaGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTg                                  <  1:1009963/100‑1 (MQ=255)
                                                         aaGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTg                                  >  1:107195/1‑100 (MQ=255)
                                                           ggATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAg                                >  2:1721676/1‑100 (MQ=255)
                                                           ggATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAg                                >  1:1112888/1‑100 (MQ=255)
                                                            gATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGc                               >  1:408749/1‑100 (MQ=255)
                                                             aTGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCGGGGCCGAGGAGATTACAGGATATGTTCAGTACGGGCGATGATATCGTCCTGAGCg                              <  1:286570/100‑1 (MQ=255)
                                                             aTGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCg                              <  2:1507930/100‑1 (MQ=255)
                                                             aTGAAATGGAGGATGTGAATGATTTCGCCGGGCGACGTCGCTGCCGGGCCGAGGAGATTACAGCATATGTTTAGTACGGGCGATGATATCGTCCTGAGCg                              <  1:1714057/100‑1 (MQ=255)
                                                              tGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCg                              >  1:1262193/1‑99 (MQ=255)
                                                               gAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTc                            <  1:580062/100‑1 (MQ=255)
                                                                aaaTTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGt                             >  1:295563/1‑98 (MQ=255)
                                                                aaaTTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCt                           >  2:2441330/1‑100 (MQ=255)
                                                                    tgatgaTGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgag                       >  2:2410251/1‑100 (MQ=255)
                                                                     gatgatGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgaga                      >  2:232891/1‑100 (MQ=255)
                                                                      atgatgTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgagag                     >  2:2320324/1‑100 (MQ=255)
                                                                        gatgTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgagag                     <  1:1950358/98‑1 (MQ=255)
                                                                        gatgTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGAAGATTACAGCATATGTTCAGGACGGGCGATGATAGCGTCCTGAGCGTCTGgagaga                    <  1:1211508/99‑1 (MQ=255)
                                                                        gatgTGAATGATGTAGACCGGCGACGGCGCCGCCGGGCCGAGGCGAGTACAGCATATGTTCAGTACTGGCGATGATATCGTCCTGAGCGTCTGGagagag                   <  1:247047/100‑1 (MQ=255)
                                                                             gAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCgg              >  1:339814/1‑100 (MQ=255)
                                                                               aTGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTg            <  1:967009/100‑1 (MQ=255)
                                                                               aTGATTTAGCCCGGCGACGACGCCGCCGGACCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTg            <  1:2143883/100‑1 (MQ=255)
                                                                                       gCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGCTGATATCGTCCTGAGCGTCTGGAGCGCGCGCGGTGAAGAAcgc    >  1:1146050/1‑100 (MQ=255)
                                                                                       gCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgc    >  2:1691189/1‑100 (MQ=255)
                                                                                       gCCCGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgc    >  2:1709725/1‑100 (MQ=255)
                                                                                        cccGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgcg   >  2:538710/1‑100 (MQ=255)
                                                                                        cccGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgcg   >  2:47362/1‑100 (MQ=255)
                                                                                        cccGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgcg   >  2:1498753/1‑100 (MQ=255)
                                                                                        cccGGCGGCGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAc      >  1:2552852/1‑97 (MQ=255)
                                                                                        cccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgcg   >  2:2593915/1‑100 (MQ=255)
                                                                                        cccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgcg   >  2:2593175/1‑100 (MQ=255)
                                                                                        cccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgcg   <  2:125139/100‑1 (MQ=255)
                                                                                         ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAcgcg   <  2:125174/99‑1 (MQ=255)
                                                                                         ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCgg  <  2:833738/100‑1 (MQ=255)
                                                                                         ccGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGCGAATAACGCgg  <  1:1360531/100‑1 (MQ=255)
                                                                                                  |                                                                                          
CCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGG  >  NC_000913/860056‑860244

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: