Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,550,574 C→G 15.1% intergenic (‑29/+67) yfeX ← / ← yfeY porphyrinogen oxidase, cytoplasmic/RpoE‑regulated lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,550,5740CG15.1% 255.5 / 35.5 93intergenic (‑29/+67)yfeX/yfeYporphyrinogen oxidase, cytoplasmic/RpoE‑regulated lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (48/31);  new base G (8/6);  total (56/37)
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.83e-01

ACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCA  >  NC_000913/2550478‑2550673
                                                                                                |                                                                                                   
aCGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAAt                                                                                                  >  1:1971722/1‑100 (MQ=255)
 cGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATc                                                                                                 >  2:881227/1‑100 (MQ=255)
 cGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATc                                                                                                 <  2:1440524/100‑1 (MQ=255)
     ggCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTg                                                                                             >  2:1703983/1‑100 (MQ=255)
      gCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTgg                                                                                            <  2:364994/100‑1 (MQ=255)
      gCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTgg                                                                                            <  1:1835104/100‑1 (MQ=255)
       cTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGc                                                                                           <  1:783863/100‑1 (MQ=255)
       cTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGc                                                                                           >  1:1173654/1‑100 (MQ=255)
       cTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGc                                                                                           <  2:1177740/100‑1 (MQ=255)
       cTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGc                                                                                           <  1:2131402/100‑1 (MQ=255)
       cTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCAACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGc                                                                                           >  1:121408/1‑100 (MQ=255)
         tCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGGATTGCGTGGCAt                                                                                         <  1:1141890/100‑1 (MQ=255)
         tCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGGATTGCGTGGCAt                                                                                         <  1:197932/100‑1 (MQ=255)
          cGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGGATTGCGTGGCAtt                                                                                        <  1:90693/100‑1 (MQ=255)
            aTCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGGATTGCGTGGCAttat                                                                                      <  1:729263/100‑3 (MQ=255)
            aTCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttgt                                                                                      <  2:1232132/100‑1 (MQ=255)
             tCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttgtt                                                                                     >  2:2089635/1‑100 (MQ=255)
              ccAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttgttg                                                                                    <  2:1200609/100‑1 (MQ=255)
               cAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGa                                                                                   >  1:2332444/1‑100 (MQ=255)
                aaaTCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAg                                                                                  <  1:1379413/100‑1 (MQ=255)
                aaaTCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAg                                                                                  <  1:11820/100‑1 (MQ=255)
                aaaTCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAg                                                                                  >  2:2111938/1‑100 (MQ=255)
                 aaTCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAgc                                                                                 <  1:2540688/100‑1 (MQ=255)
                  aTCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAgcg                                                                                >  2:1390294/1‑100 (MQ=255)
                    cgccgcGCGGCCATGTTCTGGCAAAATGACACTCTGAACCTGGGACATTGTTCCTCCTGAAAATAATAATGCCACGGATTGCGTGGCAttattattttca                                                                              <  1:801083/100‑11 (MQ=25)
                    cgccgcGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGGATTGCGTGGCAttattattttca                                                                              <  2:68284/100‑11 (MQ=37)
                      ccgcGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTAt                                                                            >  1:2030773/1‑100 (MQ=255)
                       cgcgcgGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTAtt                                                                           >  1:2370445/1‑100 (MQ=255)
                          gcgGCAATGTTCTGGCAAAATGCCACTCTGAACCTGTGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTa                                                                        <  2:1822258/100‑1 (MQ=255)
                          gcgGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTa                                                                        <  1:1142368/100‑1 (MQ=255)
                           cgGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCTCTATTTTAc                                                                       <  1:20743/100‑1 (MQ=255)
                             gCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACgt                                                                     >  1:2522591/1‑100 (MQ=255)
                                 tGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGttt                                                                 >  1:2443740/1‑100 (MQ=255)
                                  gTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGtttt                                                                >  2:805069/1‑100 (MQ=255)
                                      tGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTcg                                                            >  2:2584504/1‑100 (MQ=255)
                                         cAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCgcgg                                                         >  2:274095/1‑100 (MQ=255)
                                         cAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCgcgg                                                         >  1:2497417/1‑100 (MQ=255)
                                         cAAAACCCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCgcgg                                                         >  2:751416/1‑100 (MQ=255)
                                             aTGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCgcggcggc                                                     <  2:312012/100‑1 (MQ=255)
                                                ccACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGcggcggcgg                                                   >  2:2504304/1‑99 (MQ=255)
                                                ccACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGt                                                  >  2:695257/1‑100 (MQ=255)
                                                ccACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGt                                                  >  1:1687283/1‑100 (MQ=255)
                                                ccACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGt                                                  >  2:601595/1‑100 (MQ=255)
                                                ccACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGt                                                  >  2:760844/1‑100 (MQ=255)
                                                ccACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGt                                                  <  1:1313782/100‑1 (MQ=255)
                                                ccACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGt                                                  >  2:2218614/1‑100 (MQ=255)
                                                 cACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGtt                                                 <  2:1023368/100‑1 (MQ=255)
                                                 cACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGtt                                                 <  2:2443757/100‑1 (MQ=255)
                                                 cACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGtt                                                 >  2:1271993/1‑100 (MQ=255)
                                                    tctGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGt                                              <  1:601592/100‑1 (MQ=255)
                                                    tctGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGt                                              <  1:2475042/100‑1 (MQ=255)
                                                      tGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTc                                            >  1:1249376/1‑100 (MQ=255)
                                                      tGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATAATACGTGTTTTTTCGCGGCGGCGGTTTGTTc                                            >  1:2581188/1‑100 (MQ=255)
                                                      tGAACCTGAGACATTGTTCCTACTGAAAATAATAATGCCACGCAATCCGTGGAATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGtt                                             >  1:1980/1‑99 (MQ=255)
                                                         aCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGa                                         >  2:1483792/1‑100 (MQ=255)
                                                          ccTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGan                                        <  2:283776/100‑2 (MQ=255)
                                                          ccTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGAc                                        <  1:1832302/100‑1 (MQ=255)
                                                           cTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACg                                       <  1:760841/100‑1 (MQ=255)
                                                            tGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGc                                      >  1:1895943/1‑100 (MQ=255)
                                                               gaCATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCaaa                                   <  2:1802726/100‑1 (MQ=255)
                                                               gaCATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTGACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCaaa                                   <  1:414154/100‑1 (MQ=255)
                                                               gaCATTGTTCCTCCTGAAAATAATAATGCCACGCAACCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCaaa                                   <  2:1818599/100‑1 (MQ=255)
                                                                 cATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAAtt                                 >  1:1425006/1‑100 (MQ=255)
                                                                 cATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAAtt                                 >  1:1013984/1‑100 (MQ=255)
                                                                  aTTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTa                                >  1:1124870/1‑100 (MQ=255)
                                                                       tcctccTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAAcgcc                           <  2:317255/100‑1 (MQ=255)
                                                                          tcctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAACTTAACgccgcc                        >  2:72854/1‑100 (MQ=255)
                                                                          tcctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACgccgcc                        >  1:1436830/1‑100 (MQ=255)
                                                                           cctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCa                       <  1:1268464/100‑1 (MQ=255)
                                                                            ctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAg                      >  1:1757224/1‑100 (MQ=255)
                                                                            ctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCACGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAg                      <  1:1592520/100‑1 (MQ=255)
                                                                             tGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGa                     >  1:705669/1‑100 (MQ=255)
                                                                              gcgctcaacaatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGGTGCGTCGGGTTGTTCAGACGCAAATTAACATCGCCAg                      >  1:251875/10‑98 (MQ=21)
                                                                              gAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGAt                    >  2:594674/1‑100 (MQ=255)
                                                                                aaaTAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATaa                  >  2:1345548/1‑100 (MQ=255)
                                                                                 aataataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAAt                 >  2:981063/1‑100 (MQ=255)
                                                                                  ataataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAAtt                <  1:2017223/100‑1 (MQ=255)
                                                                                   taataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGTAGCGTTTGTTTGTTCAGAGGCAAATTAAGTCCGCCAGATAAtt                >  1:1366744/1‑99 (MQ=255)
                                                                                   taataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAAttt               >  1:645275/1‑100 (MQ=255)
                                                                                   taataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAAttt               >  2:776227/1‑100 (MQ=255)
                                                                                   taataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAAttt               >  1:1631796/1‑100 (MQ=255)
                                                                                   caacaatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAAttt               >  2:317512/5‑100 (MQ=38)
                                                                                    aacaatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAAtttt              >  2:1010311/4‑100 (MQ=39)
                                                                                     ataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTa             >  2:732542/1‑100 (MQ=255)
                                                                                     acaatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTa             >  2:2339103/3‑100 (MQ=255)
                                                                                     acaatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTa             >  1:1914462/3‑100 (MQ=255)
                                                                                      taatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTAc            <  1:2215076/100‑1 (MQ=255)
                                                                                      caatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCTGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTa             >  1:2365854/2‑99 (MQ=255)
                                                                                        atGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTg          >  2:1019366/1‑100 (MQ=255)
                                                                                        atGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGGTTTTTCGCGGCGGCGGTTTGTTAAGACGAAAATTAACGCCGCCAGATACTTTTACTg          >  1:1420983/1‑100 (MQ=255)
                                                                                           ccACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGGGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACt       >  1:116799/1‑100 (MQ=255)
                                                                                              cGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTc    >  2:811498/1‑100 (MQ=255)
                                                                                               gCAATCCGTGGCATTGTTGAGCGCCATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTcc   >  1:1218668/1‑100 (MQ=255)
                                                                                                cAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATATTACTGACTTTCCa  >  1:1310640/1‑100 (MQ=255)
                                                                                                |                                                                                                   
ACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCA  >  NC_000913/2550478‑2550673

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