Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,614,478 C→A 5.1% intergenic (+218/+202) yjjV → / ← yjjW putative DNase/putative pyruvate formate lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,614,4780CA 5.1% 251.3 / 2.1 79intergenic (+218/+202)yjjV/yjjWputative DNase/putative pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base C (44/31);  new base A (2/2);  total (46/33)
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.46e-01

ATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGAC  >  NC_000913/4614331‑4614626
                                                                                                                                                   |                                                                                                                                                    
aTCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCa                                                                                                                                                     >  1:24740/1‑149 (MQ=255)
aTCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCa                                                                                                                                                     >  6:80133/1‑149 (MQ=255)
   aTGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                                  >  7:329562/1‑149 (MQ=255)
   aTGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                                  >  7:158059/1‑149 (MQ=255)
   aTGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGAGCCGAACCGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                                  >  8:278098/1‑149 (MQ=255)
       cTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCAAAATGCCTGATGAGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                                >  6:22715/1‑147 (MQ=255)
          cAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                            <  3:127352/148‑1 (MQ=255)
          cAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                            <  8:196115/148‑1 (MQ=255)
            ggTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACACGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                         <  7:278098/149‑1 (MQ=255)
            ggTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                         <  4:77483/149‑1 (MQ=255)
            ggTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                         <  8:329562/149‑1 (MQ=255)
                  tgCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCg                                                                                                                                   >  1:220049/1‑149 (MQ=255)
                   gcagacccGTATGCCGGATTAGGCGTTCACGCCGCATCCGGCAGTTGGCGCTCAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCg                                                                                                                                   <  8:197209/142‑1 (MQ=14)
                   gCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCgg                                                                                                                                  >  4:16035/1‑149 (MQ=255)
                      gAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGgcg                                                                                                                                >  8:19042/1‑148 (MQ=255)
                       aaCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCa                                                                                                                              >  1:89439/1‑149 (MQ=255)
                          cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATaa                                                                                                                           <  5:80133/149‑1 (MQ=255)
                           gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAAt                                                                                                                          <  5:22715/149‑1 (MQ=255)
                                 cGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGAAGGCGCATAATGCCTGa                                                                                                                    >  4:118516/1‑149 (MQ=255)
                                        ggCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGGCGAACCGTAGGCCGGATAAGGCGGTCGCGCCGCATACGGGAGTCGGCGCATAATGACTGATGCGACg                                                                                                             >  6:248905/1‑149 (MQ=14)
                                        ggCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGAc                                                                                                              >  6:98147/1‑148 (MQ=255)
                                        ggCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGAc                                                                                                              <  5:98147/148‑1 (MQ=255)
                                              cACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCGTggc                                                                                                       >  5:87921/1‑147 (MQ=255)
                                                  ccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                       >  2:222115/1‑113 (MQ=255)
                                                  ccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                       <  1:222115/113‑1 (MQ=255)
                                                             aGTTGGCTCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGcc                                                                                        <  8:312456/149‑1 (MQ=255)
                                                             aGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGcc                                                                                        <  7:19042/149‑1 (MQ=255)
                                                                  gcgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAg                                                                                   >  6:116751/1‑149 (MQ=255)
                                                                   cgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                         <  6:334016/94‑1 (MQ=14)
                                                                   cgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                         <  6:328890/94‑1 (MQ=14)
                                                                   cgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                         >  5:328890/1‑94 (MQ=14)
                                                                   cgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                         >  5:334016/1‑94 (MQ=14)
                                                                    gcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGt                                                                                 >  7:19669/1‑149 (MQ=255)
                                                                      acaATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATAAGGCGTTCACGCCGAATCCGGCAGTCGGCGAATAAAGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGttt                                                                               >  8:142162/1‑149 (MQ=255)
                                                                      acaATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGttt                                                                               >  4:92884/1‑149 (MQ=255)
                                                                      acaATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGttt                                                                               <  3:16035/149‑1 (MQ=255)
                                                                        aaTGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTtgt                                                                             <  3:92884/149‑1 (MQ=255)
                                                                         aTGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCTCGTCTTATCATGCCTACAGGTTtgtg                                                                            <  8:158059/149‑1 (MQ=255)
                                                                         aTGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTtgtg                                                                            <  2:220049/149‑1 (MQ=255)
                                                                             cTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCTTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCg                                                                         >  6:40770/1‑148 (MQ=255)
                                                                                aTGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAAcc                                                                     >  4:85863/1‑149 (MQ=255)
                                                                                    gACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAg                                                                 >  1:66477/1‑149 (MQ=255)
                                                                                        cTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCTGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCg                                                             <  2:24740/149‑1 (MQ=255)
                                                                                        cTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCg                                                             <  2:89439/149‑1 (MQ=255)
                                                                                          tGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGa                                                           >  7:166859/1‑149 (MQ=255)
                                                                                             cgcgTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGAt                                                          >  7:175077/1‑147 (MQ=255)
                                                                                                  cTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCAGTCGGAGCATAAAGCCTGATGAGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGt                                                   >  4:118167/1‑149 (MQ=255)
                                                                                                   ttATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGtt                                                  >  1:273882/1‑149 (MQ=255)
                                                                                                     aTCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGGCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTc                                                 <  1:24384/148‑1 (MQ=255)
                                                                                                          gCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGc                                                                                                                        >  6:284330/1‑72 (MQ=255)
                                                                                                          gCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGc                                                                                                                        <  5:284330/72‑1 (MQ=255)
                                                                                                            cTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGt                                                                                 <  2:230464/109‑1 (MQ=255)
                                                                                                            cTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGt                                                                                 >  1:230464/1‑109 (MQ=255)
                                                                                                            cTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCa                                         >  5:18937/1‑149 (MQ=255)
                                                                                                            cTACAAGTCTGTGCCGAACAGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGACGCGTCTAATCATGCCAACAGCTATGTGACGAACCGAAGGCCCGATAAGGCCTTCGcgccgc                                          >  6:79099/1‑148 (MQ=14)
                                                                                                              aCAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCGTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATc                                       >  3:279557/1‑149 (MQ=255)
                                                                                                               cAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATcc                                      >  2:292156/1‑149 (MQ=255)
                                                                                                                  gTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAACGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCgg                                    >  8:322496/1‑148 (MQ=255)
                                                                                                                   tCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCa                                  <  1:292156/149‑1 (MQ=255)
                                                                                                                                gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAACGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAAt                     <  7:322496/149‑1 (MQ=255)
                                                                                                                                 tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCg                                     <  7:109266/132‑1 (MQ=255)
                                                                                                                                 tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATg                    >  3:211401/1‑149 (MQ=255)
                                                                                                                                 tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCg                                     >  8:109266/1‑132 (MQ=17)
                                                                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGcc                  >  6:215469/1‑149 (MQ=255)
                                                                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGc                   <  8:19669/148‑1 (MQ=255)
                                                                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACATGTTTGTTCCGAACCTTCCGCCTGATAAGGCGTTCTCGCCGCATCCGGCAGTTGGCGCACAATGCCt                 <  6:87921/149‑1 (MQ=14)
                                                                                                                                      cGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCCCAATGCCtta               >  7:281607/1‑147 (MQ=255)
                                                                                                                                       ggATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGAt              <  8:175077/149‑1 (MQ=255)
                                                                                                                                             ggCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGc                                                                                                            >  8:72809/1‑49 (MQ=255)
                                                                                                                                             ggCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGc                                                                                                            <  7:72809/49‑1 (MQ=255)
                                                                                                                                             ggCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATTCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGAGCCGCATCCGGCAGTTGGCGCACAATGCATGATGCg           >  5:279204/1‑146 (MQ=255)
                                                                                                                                             ggCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCCTTCGCGCCGCATCCGGCAGTTGGCGCACAAAGCCTGATGCGACg        >  6:17445/1‑149 (MQ=255)
                                                                                                                                              gCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGc       >  5:225879/1‑149 (MQ=255)
                                                                                                                                                gTTCACGCAGCATCCGGCAGTCGGCGCATAATTCCTGATGGGTCGCTTGTCGTGTATTATCATGCCTACAGGTTTTTGCCGATCCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCtt     <  8:323659/149‑1 (MQ=14)
                                                                                                                                                 ttCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACaa                      >  5:64064/1‑131 (MQ=255)
                                                                                                                                                 ttCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACaa                      <  6:64064/131‑1 (MQ=255)
                                                                                                                                                 ttCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGcc                  >  1:20367/1‑135 (MQ=255)
                                                                                                                                                 ttCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGcc                  <  2:20367/135‑1 (MQ=255)
                                                                                                                                                 ttCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCAACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTg    >  4:315844/1‑149 (MQ=255)
                                                                                                                                                   cACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGAc  >  3:85214/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGAC  >  NC_000913/4614331‑4614626

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: