Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,550,578 C→G 12.6% intergenic (‑33/+63) yfeX ← / ← yfeY porphyrinogen peroxidase/DUF1131 domain‑containing lipoprotein YfeY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,550,5780CG12.6% 291.9 / 30.8 103intergenic (‑33/+63)yfeX/yfeYporphyrinogen peroxidase/DUF1131 domain‑containing lipoprotein YfeY
Reads supporting (aligned to +/- strand):  ref base C (52/37);  new base G (8/5);  total (60/43)
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.16e-01

Only 100 of 105 total aligned reads displayed.

GTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCC  >  NC_000913/2550429‑2550727
                                                                                                                                                     |                                                                                                                                                     
gTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATc                                                                                                                                                       <  2:1164086/150‑1 (MQ=255)
    gCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTg                                                                                                                                                   >  2:1266929/1‑150 (MQ=255)
     cAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTgg                                                                                                                                                  >  1:877583/1‑150 (MQ=255)
     cAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTgg                                                                                                                                                  <  2:21479/150‑1 (MQ=255)
     cAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTgg                                                                                                                                                  <  1:1373905/150‑1 (MQ=255)
      aaaTGTTTTTCTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAAGCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGc                                                                                                                                                 <  1:379439/150‑1 (MQ=255)
         tGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAtt                                                                                                                                              >  1:1483256/1‑150 (MQ=255)
         tGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATACCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGGATTGCGTGGCAtt                                                                                                                                              <  2:1531325/150‑1 (MQ=255)
          gTTTTACTGGCCGCACGCAGGGCGTCACCTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCGCCTGAAAATAATAATGCCACGCAATCCGTGGCAttg                                                                                                                                             >  1:762481/1‑150 (MQ=255)
           ttttACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttgt                                                                                                                                            <  1:417738/150‑1 (MQ=255)
           ttttACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAACAATGCCACGGATTGCGTGGCAttat                                                                                                                                            <  1:1069905/150‑3 (MQ=255)
            tttACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGGATTGCGTGGCattatt                                                                                                                                           <  2:1257540/150‑1 (MQ=255)
             ttACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttgttg                                                                                                                                          >  2:56109/1‑150 (MQ=255)
             ttACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttgttg                                                                                                                                          <  2:1521907/150‑1 (MQ=255)
               aCTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGCACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAg                                                                                                                                        >  2:1447911/1‑150 (MQ=255)
                cTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGGATTGCGTGGCattattattt                                                                                                                                       <  1:559226/150‑5 (MQ=255)
                cTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAgc                                                                                                                                       >  2:547156/1‑150 (MQ=255)
                         cgcaGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGCCATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTa                                                                                                                              >  2:1401955/1‑150 (MQ=255)
                           caGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACg                                                                                                                            <  1:11543/150‑1 (MQ=255)
                              ggCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACgtgt                                                                                                                         >  1:1126590/1‑150 (MQ=255)
                                cGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGttt                                                                                                                       >  2:449490/1‑150 (MQ=255)
                                 gTCAACTTCCCCTTTCAGGTTGGCTTCGATCCAAATCGCCGGGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGtttt                                                                                                                      <  2:1299344/150‑1 (MQ=255)
                                   cAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGtttttt                                                                                                                    >  1:1365357/1‑150 (MQ=255)
                                   cAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGtttttt                                                                                                                    <  1:359553/150‑1 (MQ=255)
                                    aaCTTCCGCTTTCAGGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGTCATTGTTTTTCATGAAAATAATAATGCCTCGAAATACGTGGAATTTTTGAGTGCTATTTTACGTGttttttt                                                                                                                   <  2:689396/150‑2 (MQ=255)
                                    aaCTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTc                                                                                                                   >  2:108548/1‑150 (MQ=255)
                                     aCTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAgcgc                                                                                                                                     >  2:1479169/1‑131 (MQ=255)
                                     aCTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAgcgc                                                                                                                                     <  1:1479169/131‑1 (MQ=255)
                                      cTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCCCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTcgc                                                                                                                 <  2:762481/150‑1 (MQ=255)
                                        tCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCgcgg                                                                                                               >  1:376557/1‑150 (MQ=255)
                                        tCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCgcgg                                                                                                               >  2:223932/1‑150 (MQ=255)
                                           ccTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCgcggcgg                                                                                                            >  1:666472/1‑150 (MQ=255)
                                           ccTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCgcggcgg                                                                                                            <  2:655620/150‑1 (MQ=255)
                                            cTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCgcggcggc                                                                                                           >  2:517875/1‑150 (MQ=255)
                                              ttCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGcggcggcgg                                                                                                         >  2:1492006/1‑150 (MQ=255)
                                                cACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGtt                                                                                                       <  2:622077/150‑1 (MQ=255)
                                                 aCGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGttt                                                                                                      >  1:1182085/1‑150 (MQ=255)
                                                 aCGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGttt                                                                                                      >  2:676562/1‑150 (MQ=255)
                                                    ttGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGtt                                                                                                   <  2:1094181/150‑1 (MQ=255)
                                                    ttGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGtt                                                                                                   >  1:1196426/1‑150 (MQ=255)
                                                     tGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTc                                                                                                  >  2:1402848/1‑150 (MQ=255)
                                                       gCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGa                                                                                               <  1:923596/150‑1 (MQ=255)
                                                               ccAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAAt                                                                                        >  2:510044/1‑150 (MQ=255)
                                                               ccAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAAt                                                                                        >  2:217754/1‑150 (MQ=255)
                                                               ccAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAAt                                                                                        >  1:94229/1‑150 (MQ=255)
                                                                   aTCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAAc                                                                                    >  1:70450/1‑150 (MQ=255)
                                                                         gcgcgGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACgccgcc                                                                              <  2:1365357/150‑1 (MQ=255)
                                                                         gcgcgGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACgccgcc                                                                              >  2:133967/1‑150 (MQ=255)
                                                                                  tGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAAttt                                                                     <  2:70450/150‑1 (MQ=255)
                                                                                     tCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTAc                                                                  <  2:666472/150‑1 (MQ=255)
                                                                                     tCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTAc                                                                  <  1:547156/150‑1 (MQ=255)
                                                                                         tcAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGAc                                                              <  2:1004446/149‑1 (MQ=255)
                                                                                         gCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGAc                                                              >  2:1490156/1‑150 (MQ=255)
                                                                                         gCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGAc                                                              >  1:244102/1‑150 (MQ=255)
                                                                                         gCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGAc                                                              >  2:1491088/1‑150 (MQ=255)
                                                                                             aaTGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTc                                                          >  2:1144634/1‑150 (MQ=255)
                                                                                                      ctGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTc                                                 >  1:383768/1‑150 (MQ=255)
                                                                                                      ctGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTc                                                 <  1:56109/150‑1 (MQ=255)
                                                                                                      ctGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCAAGttttta                                                 >  1:1048424/1‑149 (MQ=255)
                                                                                                       tGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCa                                                >  2:111242/1‑150 (MQ=255)
                                                                                                        gAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAg                                               <  1:547842/150‑1 (MQ=255)
                                                                                                          aCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAggg                                             >  2:71672/1‑150 (MQ=255)
                                                                                                          aCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAggg                                             <  1:979676/150‑1 (MQ=255)
                                                                                                          aCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAggg                                             <  1:200186/150‑1 (MQ=255)
                                                                                                           ccTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGt                                            >  1:1307067/1‑150 (MQ=255)
                                                                                                                gaCATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGt                                       >  2:937519/1‑150 (MQ=255)
                                                                                                                  cATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTcc                                     <  1:133967/150‑1 (MQ=255)
                                                                                                                  cATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTcc                                     >  1:704966/1‑150 (MQ=255)
                                                                                                                     tGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGaa                                  >  1:1248206/1‑150 (MQ=255)
                                                                                                                      gTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGGCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAg                                 >  2:1202971/1‑150 (MQ=255)
                                                                                                                      gTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAg                                 <  1:279588/150‑1 (MQ=255)
                                                                                                                       ttCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAgg                                <  1:499410/150‑1 (MQ=255)
                                                                                                                            cctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGa                                   >  2:702979/1‑142 (MQ=255)
                                                                                                                            cctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGa                                   <  1:702979/142‑1 (MQ=255)
                                                                                                                            cctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTa                           >  2:433652/1‑150 (MQ=255)
                                                                                                                               gAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAAcc                        <  1:223932/150‑1 (MQ=255)
                                                                                                                               gAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAAcc                        >  2:311961/1‑150 (MQ=255)
                                                                                                                                 gctcaacaatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCt                      >  1:882575/8‑150 (MQ=255)
                                                                                                                                 aaaTAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCt                      <  2:1182085/150‑1 (MQ=255)
                                                                                                                                    caacaatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCa                   >  1:1490340/5‑150 (MQ=255)
                                                                                                                                    caacaatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACACAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCa                   >  2:281579/5‑150 (MQ=255)
                                                                                                                                    caacaaagCCACGGATTGCGTTGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCAGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCa                   >  2:509198/8‑150 (MQ=255)
                                                                                                                                     aataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAg                  <  1:88715/150‑1 (MQ=255)
                                                                                                                                      ataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAgg                 >  1:1457433/1‑150 (MQ=255)
                                                                                                                                      acaatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAgg                 >  1:946574/3‑150 (MQ=255)
                                                                                                                                       caatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGa                >  1:803135/2‑150 (MQ=255)
                                                                                                                                       caatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGa                >  1:1440469/2‑150 (MQ=255)
                                                                                                                                       caatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGa                >  2:1613781/2‑150 (MQ=255)
                                                                                                                                        aatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGAc               >  1:1160172/1‑150 (MQ=255)
                                                                                                                                        aatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGAc               >  1:1420063/1‑150 (MQ=255)
                                                                                                                                          tGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCa             >  2:1266541/1‑150 (MQ=255)
                                                                                                                                           gCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCAc            <  2:1490340/150‑1 (MQ=255)
                                                                                                                                           gCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCa                   >  2:499987/1‑143 (MQ=255)
                                                                                                                                           gCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCa                   <  1:499987/143‑1 (MQ=255)
                                                                                                                                           gCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCAc            >  1:654090/1‑150 (MQ=255)
                                                                                                                                            ccACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACg           >  2:557645/1‑150 (MQ=255)
                                                                                                                                             cACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGc          <  1:1402848/150‑1 (MQ=255)
                                                                                                                                             cACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGc          >  2:1261757/1‑150 (MQ=255)
                                                                                                                                             cACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGc          >  2:899301/1‑150 (MQ=255)
                                                                                                                                                    tCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTccc   <  2:1296670/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCC  >  NC_000913/2550429‑2550727

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