Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,877,871 A→T 6.7% intergenic (+253/+698) iap → / ← ygbF aminopeptidase in alkaline phosphatase isozyme conversion/CRISPR adaptation ssRNA endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,877,8710AT 6.7% 267.0 / 3.3 86intergenic (+253/+698)iap/ygbFaminopeptidase in alkaline phosphatase isozyme conversion/CRISPR adaptation ssRNA endonuclease
Reads supporting (aligned to +/- strand):  ref base A (40/40);  new base T (3/3);  total (43/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 = 1.00e+00

GGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCC  >  NC_000913/2877723‑2878015
                                                                                                                                                    |                                                                                                                                                
ggAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCACGGTTCAGGCGTTGCAAAACTGGCTACCGGGCGGTTTATACCCGCTAACGCGGGGAACTCGTATTCCATAATTCCACCTa                                                                                                                                                  >  8:151852/1‑149 (MQ=255)
 gAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTAt                                                                                                                                                 <  1:87132/149‑1 (MQ=255)
 gAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTAt                                                                                                                                                 <  5:273978/149‑1 (MQ=255)
 gAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCGGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTAt                                                                                                                                                 <  5:5138/149‑1 (MQ=255)
    cacCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTc                                                                                                                                              <  5:125221/149‑1 (MQ=255)
       cAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCAGCTGGCGCGGGGAACTCTCGGTTCAGTCGTTGCAAACCAGGCTACCGGGCGGATTATCCCCGCTAACGCGGGGAACTCGTATTCCATCATTCCACCTATGTCTg                                                                                                                                            >  6:143449/1‑148 (MQ=255)
              aGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCg                                                                                                                                    <  5:52111/149‑1 (MQ=255)
                gCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTTTGTCTGAACTCCCGGt                                                                                                                                  >  8:58462/1‑149 (MQ=255)
                gCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTTTGTCTGAACTCCCGGt                                                                                                                                  >  8:52027/1‑149 (MQ=255)
                gCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGt                                                                                                                                  >  4:339149/1‑149 (MQ=255)
                          tcACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCATGCGTTGCAAACCTGGCGACCGGGCGGTTTATCCCCGCTATCGCGGGGAACTCGTATTCCATCATTCCACCTATCTCTGAACTCCCGGTATATcccc                                                                                                                          <  2:313400/147‑1 (MQ=255)
                          tcACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTAt                                                                                                                              >  5:140346/1‑143 (MQ=255)
                          tcACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTAt                                                                                                                              <  6:140346/143‑1 (MQ=255)
                            aCCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTg                                                                                                                      >  7:311269/1‑149 (MQ=255)
                             ccGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCACCGCTAACGCGGGGAACTCGTAGTCCATCTTTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCggg                                                                                                                     >  5:61664/1‑149 (MQ=255)
                                 cgtTGCCGGCTTATCCCTGCTGGCGCGCGGAACTCTCGGTTCAGTCTTTGCAAACCTGGCGACCGGGCTGTTTATCCCCGCTAACGCGGGGAACTCTTAGTCCATCATTCCACATATTCCTTACCTGCCGGTTTATCCCCGCTGGcgcg                                                                                                                 <  6:61664/149‑1 (MQ=255)
                                   ttGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCggg                                                                                                               <  3:374185/149‑1 (MQ=255)
                                       cGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGATCATTCGTTTCAAGCCTGCCACCCGGGCGGTTTATCCCCGCTAAAGCGGGGAACTCGTAGTCGCTCATTACACCTAGGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAAc                                                                                                           >  4:54769/1‑149 (MQ=255)
                                       cGGTTTATCCCTGCTGGCGCGGGGAACTCGCGGTTCAGGCGTTTGAGAAAGGGGAACCCGCCGGTTTATCCCCGCTAGCGCGGGGAACTGGTAGACCAACATTCCAGCTATGGCTGTACTCCGGGGTTATCCCCGCTGGGGCGGGGacc                                                                                                           >  4:9103/1‑147 (MQ=255)
                                              tCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCggg                                                                                                    >  8:226400/1‑149 (MQ=255)
                                               cccTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCgggg                                                                                                              >  4:176786/1‑138 (MQ=255)
                                               cccTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCgggg                                                                                                              <  3:176786/138‑1 (MQ=255)
                                                ccTGCTGGCGCGGGGAACTCTCGTTTCAGGCGATGCAAAGATAGATACCGGCCGGTTTATCCCCGCGAACGCGGGGAACTCGTAGTGCATAATTACCCCTATGTCTGAACTCCCGGTTTATAACAGATGTCTAGGGGAACTCCCggggg                                                                                                  >  8:48478/1‑149 (MQ=255)
                                                ccTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGACCTCCCggggg                                                                                                  >  7:393884/1‑149 (MQ=255)
                                                      ggCGCGGGGAACTCTCGGTTCAGGCGTTGCAAAACTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATg                                                                                            >  5:390051/1‑149 (MQ=255)
                                                        cgctGGGAACTGTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACTCGGGGAACTCGTTGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGtt                                                                                          <  7:151852/149‑1 (MQ=255)
                                                               aaCTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTc                                                                                   <  1:233153/149‑1 (MQ=255)
                                                               aaCTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGATTATCCCCGATAACGCGGCGAACTCCTAGTCCATCATTCCACCTATGTCTGAACTCACGGTTTATCCCCTCTGGCGCGGGGAACTCCAGTGGGATAATGTTTACTGt                                                                                    >  8:105065/1‑148 (MQ=255)
                                                                aCTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCa                                                                                  <  5:57913/149‑1 (MQ=255)
                                                                 ctctCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCgg                                                                                                     >  7:186957/1‑129 (MQ=255)
                                                                 ctctCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCgg                                                                                                     <  8:186957/129‑1 (MQ=255)
                                                                  tctcGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTTTGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATg                                                                                >  5:181142/1‑149 (MQ=255)
                                                                    tcGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATgcg                                                                              >  6:4183/1‑149 (MQ=255)
                                                                           aGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCg                                                                       <  5:107469/149‑1 (MQ=255)
                                                                             gCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGt                                                                     <  3:339149/149‑1 (MQ=255)
                                                                              cGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGtt                                                                    >  8:110318/1‑149 (MQ=255)
                                                                                  gCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTAt                                                                 <  5:4183/148‑1 (MQ=255)
                                                                                  gCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATc                                                                >  5:106079/1‑149 (MQ=255)
                                                                                  gCAAACCTGGCTACCGGGCGGTTTATCCCCGCCAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTCTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATc                                                                >  7:226266/1‑149 (MQ=255)
                                                                                        cTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCt                                                          <  7:226400/149‑1 (MQ=255)
                                                                                        cTGGCGACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCt                                                          <  5:112950/149‑1 (MQ=255)
                                                                                           gCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGgc                                                       >  7:154186/1‑149 (MQ=255)
                                                                                             tACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCgg                                                                                                     >  8:178214/1‑101 (MQ=255)
                                                                                             tACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCgg                                                                                                     <  7:178214/101‑1 (MQ=255)
                                                                                             tACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGgcgc                                                     >  6:250768/1‑149 (MQ=255)
                                                                                               ccGGGCGTTTTATCCCCGCTAACGCGGTGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCATGTTTATCCCCGCTGTCGCGTGGAACTCCAGGGTGATAATGTTTACGGTGATGCGCCCCCCGTTTTATCCCCGCTGGcgcg                                                    <  7:105065/148‑1 (MQ=255)
                                                                                               ccGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGtt                                                                    >  4:125661/1‑132 (MQ=255)
                                                                                               ccGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGtt                                                                    <  3:125661/132‑1 (MQ=255)
                                                                                               ccGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGgcgc                                                     <  8:259135/147‑1 (MQ=255)
                                                                                               ccGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCgg                                                   >  3:259761/1‑149 (MQ=255)
                                                                                               ccGGGCGGGTAATCCCCGCTAACGCGGGGAACTCGGAGTCCCTCATTCCACCTATGTCTGAACTCCCGGGTTCTCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGCTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCgg                                                   <  3:9103/149‑1 (MQ=255)
                                                                                                cGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTTTGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCggg                                                  <  6:181142/149‑1 (MQ=255)
                                                                                                cGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCggg                                                  <  4:259761/149‑1 (MQ=255)
                                                                                                 gggCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCggg                                                  >  4:220795/1‑148 (MQ=255)
                                                                                                      gTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAActc                                            >  6:64219/1‑149 (MQ=255)
                                                                                                       tttATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAActct                                           >  3:11886/1‑149 (MQ=255)
                                                                                                        ttATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATTATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTg                                          >  3:270802/1‑149 (MQ=255)
                                                                                                          aTCCCCGCTAACGCGGGGAACTCGTAGGCCATCATTCCACCTATGTCTGCACTCCCGGTTTATCCCCGCTGCCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTgg                                         <  2:264714/148‑1 (MQ=255)
                                                                                                           tCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGggc                                       >  4:265251/1‑149 (MQ=255)
                                                                                                             cccGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGggcgg                                     >  1:92274/1‑149 (MQ=255)
                                                                                                               cGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGct                                   <  2:92274/149‑1 (MQ=255)
                                                                                                               cGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGct                                   >  2:244352/1‑149 (MQ=255)
                                                                                                                    aCGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGcttgcct                              >  8:291563/1‑149 (MQ=255)
                                                                                                                    aCGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGcttgc                                <  3:220795/147‑1 (MQ=255)
                                                                                                                      gcgGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGcttgccttg                            <  8:154186/149‑1 (MQ=255)
                                                                                                                      gcgGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGcttgccttg                            >  5:199127/1‑149 (MQ=255)
                                                                                                                        ggggAACTCGTAGTCCATCATTCCACCTTTGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGca                          <  7:52027/149‑1 (MQ=255)
                                                                                                                        ggggAACTCGTAGTCCATCATTCCACCTTTGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGca                          <  7:58462/149‑1 (MQ=255)
                                                                                                                        ggggAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTCTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGca                          <  8:226266/149‑1 (MQ=255)
                                                                                                                          ggAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagc                        <  6:199127/149‑1 (MQ=255)
                                                                                                                              cTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCAt                                                                                 <  8:16393/88‑1 (MQ=255)
                                                                                                                              cTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCAt                                                                                 >  7:16393/1‑88 (MQ=255)
                                                                                                                              cTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagccagc                    >  4:316555/1‑149 (MQ=255)
                                                                                                                                 gTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTcc                 <  1:244352/149‑1 (MQ=255)
                                                                                                                                  tAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTc                  >  1:208265/1‑147 (MQ=255)
                                                                                                                                  tAGGCCATCATTCCACCTATGTCTGAACTCCCGGTTTCTCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGGTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCca                <  3:54769/149‑1 (MQ=255)
                                                                                                                                     tCCATCATTCCAACTATGTCTGAACTCCCGGTGAAGCCCCGCTGGCGCGGTTAACTCCCGGGGGCTAATGTTTACGGTCATGCGCCCCCCGGGTTATCCCCGCTGGCGCGGGGAACTCTTggc                                       <  1:400159/123‑1 (MQ=255)
                                                                                                                                       catcatTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCagcagc           <  6:106079/149‑1 (MQ=255)
                                                                                                                                        atcatTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTggg                                        <  6:75981/119‑1 (MQ=255)
                                                                                                                                        atcatTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTggg                                        >  5:75981/1‑119 (MQ=255)
                                                                                                                                        atcatTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCg          <  2:208265/149‑1 (MQ=255)
                                                                                                                                        atcatTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCg          >  5:273254/1‑149 (MQ=255)
                                                                                                                                            ttCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGcttgcc                               <  6:55658/124‑1 (MQ=255)
                                                                                                                                            ttCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGcttgcc                               >  5:55658/1‑124 (MQ=255)
                                                                                                                                               cACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGTGAACTCCCGGGGGATAATGTATAAGGTCATGCGCCCCCCGGTTTATCCCCGCTGCCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGATTATc   >  3:150478/1‑149 (MQ=255)
                                                                                                                                               cACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATc   <  6:253144/149‑1 (MQ=255)
                                                                                                                                                aCCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATcc  >  8:106599/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCC  >  NC_000913/2877723‑2878015

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