Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 860,151 A→G 18.4% 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,1510AG18.4% 186.5 / 38.0 94intergenic (‑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 (32/44);  new base G (7/10);  total (39/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

AAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGG  >  NC_000913/860053‑860244
                                                                                                  |                                                                                             
aaaCCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGcgac                                                                                              >  2:1385676/1‑100 (MQ=255)
   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:75324/100‑2 (MQ=255)
              taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGct                                                                                <  2:5412/100‑2 (MQ=255)
              taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGct                                                                                <  2:1361499/100‑2 (MQ=255)
              taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGct                                                                                <  2:75257/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:60195/100‑1 (MQ=255)
                   aCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgagga                                                                           <  1:1443069/100‑1 (MQ=255)
                      aTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGgaga                                                                         <  2:2474462/99‑1 (MQ=255)
                       tgtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGAtt                                                                       >  2:2501438/1‑100 (MQ=255)
                       tgtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGAtt                                                                       <  1:2052394/100‑1 (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                                                                   <  1:2501419/100‑1 (MQ=255)
                           aCGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAg                                                                   <  2:1265238/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                                                <  2:1072833/100‑1 (MQ=255)
                                              tttCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCg                                                <  1:2337381/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)
                                                                                                  |                                                                                             
AAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGG  >  NC_000913/860053‑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: