Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	868,052	A→C	5.8%	intergenic (+150/+27)	PP_0750 → / ← mqo‑I	hypothetical protein/malate:quinone oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	868,052	0	A	C	5.8%	242.8 / 4.7	81	intergenic (+150/+27)	PP_0750/mqo‑I	hypothetical protein/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (43/33);  new base C (3/2);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

CAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGG  >  NC_002947/867911‑868196                                                                                                                                                 |                                                                                                                                                   cAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCGGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑Gt                                                                                                                                            >  5:191740/1‑149 (MQ=255)    tGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCCCAGGCGTGGCGCCGCCCCAGCGGGGAGTTCCGCCCCGGGGGGGGCGGGGCTGGCCGCGCAGAGGGGGGGGCAGGGG‑‑‑GGGGAAGCGCCCGGCACCAAAAACGCCCCGCGGGGG‑Gtttt                                                                                                                                         >  5:7957/1‑149 (MQ=255)    ggCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACCGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGTCACAAAAAACGCCCCTCGGGGC‑Gtttt                                                                                                                                         <  6:39763/148‑1 (MQ=255)      ccGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAGCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTg                                                                                                                                       <  6:311082/149‑1 (MQ=255)      ccGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCTCGGGGC‑Gtttttt                                                                                                                                       <  8:283113/149‑2 (MQ=255)           ctctTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGGGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGttt                                                                                                                                  >  1:305461/1‑149 (MQ=255)              ttCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCg                                                                                                                                     <  7:294688/143‑1 (MQ=255)              ttCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑c                                                                                                                                <  1:50873/147‑1 (MQ=255)              ttCGCGGGCGCGCCCCCTCCCACAGGAGGCGCGCCGCACCTGCGGGCATTACCGAACCGGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCG‑‑‑GGTCAAGCGCCCGGCACAAAAAAACCCCCGAGGGGCGGTTTTTGctagctc                                                                                                                                >  8:294688/1‑143 (MQ=255)                ggcgGGGGCGCCCGGCCCCCCCGGCGGCGCGCCGCACCTGGGGGCATTACCGAACCTGTGGGTGGGGGTATGCCCGCGAAGTGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGTCACAAAAATCGCCCCGAGGGGC‑GTTTTTGCGTTT‑Ca                                                                                                                               <  7:272505/145‑1 (MQ=255)                         gcCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGttt                                                                                                                                  >  2:392120/1‑135 (MQ=255)                         gcCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGttt                                                                                                                                  <  1:392120/135‑1 (MQ=255)                         gcCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCCGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAg                                                                                                                   >  2:90702/1‑149 (MQ=255)                             gCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTa                                                                                                               <  6:323620/149‑1 (MQ=255)                               tCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGt                                                                                                             >  1:70282/1‑149 (MQ=255)                                 ccACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGc                                                                                                                                      >  5:332421/1‑123 (MQ=255)                                 ccACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGc                                                                                                                                      <  6:332421/123‑1 (MQ=255)                                       gAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGgcc                                                                                                     >  2:207616/1‑149 (MQ=255)                                            gcgcCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGgccg                                                                                                    <  1:212559/145‑1 (MQ=255)                                            gcgcCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGgccg                                                                                                    >  2:212559/1‑145 (MQ=255)                                            gcgcCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCgg                                                                                                >  6:36524/1‑149 (MQ=255)                                              gccgcACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTa                                                                                                               <  1:251585/132‑1 (MQ=255)                                              gccgcACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTa                                                                                                               >  2:251585/1‑132 (MQ=255)                                               ccgcACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGc                                                                                                           <  2:321780/135‑1 (MQ=255)                                               ccgcACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGc                                                                                                           >  1:321780/1‑135 (MQ=255)                                               ccgcACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTg                                                                                                        >  7:182730/1‑138 (MQ=255)                                               ccgcACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTg                                                                                                        <  8:182730/138‑1 (MQ=255)                                                cgcACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGgctg                                                                                            <  1:83729/149‑1 (MQ=255)                                                       gCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTg                                                                                     <  1:207616/149‑1 (MQ=255)                                                       gCGGGAATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTg                                                                                     >  6:67665/1‑149 (MQ=255)                                                        cGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTTCCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGctctct                                                                                      <  1:83783/148‑1 (MQ=255)                                                        cGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGctctct                                                                                      <  6:170739/147‑1 (MQ=255)                                                          ggCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGc                                                                                    >  4:406674/1‑147 (MQ=255)                                                             aTTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGt                                                                               >  7:288913/1‑149 (MQ=255)                                                             aTTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGt                                                                               >  8:131844/1‑149 (MQ=255)                                                             aTTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGt                                                                               >  5:229086/1‑149 (MQ=255)                                                                 ccGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCCGGCACAAAAAACGCCCCGAGGGGG‑GTTTTTGCGTTA‑CAGACCGGCTTAAGCGTAGGGATTGGCCGCCGGGAGGCTATCTGCAAGGTCCAg                                                                           >  4:311615/1‑149 (MQ=38)                                                                       cTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAgg                                                                          <  4:57785/144‑1 (MQ=255)                                                                       cTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAgg                                                                          >  3:57785/1‑144 (MQ=255)                                                                        tgtgGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCcgcg                                                                    >  2:345931/1‑149 (MQ=255)                                                                           gggAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAg                                                                                 >  6:263173/1‑133 (MQ=255)                                                                           gggAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAg                                                                                 <  5:263173/133‑1 (MQ=255)                                                                             gAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGcgtcgt                                                                >  2:394390/1‑148 (MQ=255)                                                                              aGCGGGCGCGCCCGCGAAGAGGCCGGCACTGGCTACAGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGcgtcg                                                                 <  5:391454/149‑1 (MQ=37)                                                                                  ggCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGcc                                                                       >  3:72063/1‑136 (MQ=255)                                                                                  ggCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGcc                                                                       <  4:72063/136‑1 (MQ=255)                                                                                     aTGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCAc                                                       <  6:191740/149‑1 (MQ=255)                                                                                       acccGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCTCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACt                                                      >  4:234040/2‑148 (MQ=255)                                                                                          cgcgAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATc                                                  >  7:234535/1‑149 (MQ=255)                                                                                             gAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTtg                                               <  6:7957/149‑1 (MQ=255)                                                                                              aaGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTtgt                                              <  1:345931/149‑1 (MQ=255)                                                                                                  gggCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTa                                          <  5:36524/149‑1 (MQ=255)                                                                                                   ggCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATc                                                            <  7:294348/130‑1 (MQ=255)                                                                                                   ggCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATc                                                            >  8:294348/1‑130 (MQ=255)                                                                                                    gCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACa                                        <  2:305461/149‑1 (MQ=255)                                                                                                        tACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAgcg                                     >  2:226192/1‑148 (MQ=255)                                                                                                          cAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCgca                                   <  3:406674/148‑1 (MQ=255)                                                                                                             gCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGACGCGTCGTTATCTGCACTGATCTTGTGGTACAGCGCAGCAt                               >  6:417164/1‑149 (MQ=255)                                                                                                              cT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCATGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATc                              <  1:90702/149‑1 (MQ=255)                                                                                                               t‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGcc                                                                                   <  8:387946/95‑1 (MQ=255)                                                                                                               t‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGcc                                                                                   >  7:387946/1‑95 (MQ=255)                                                                                                               t‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCg                             >  7:108070/1‑149 (MQ=255)                                                                                                                   gATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCgg                            >  5:297754/1‑149 (MQ=255)                                                                                                                     tCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCGTCGGtt                          >  7:282206/1‑149 (MQ=255)                                                                                                                     tCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGtt                          <  6:162246/149‑1 (MQ=255)                                                                                                                       aaGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGc                        >  1:270190/1‑149 (MQ=255)                                                                                                                       aaGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGc                        >  2:269041/1‑149 (MQ=255)                                                                                                                           gcCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAgg                    >  2:339996/1‑149 (MQ=255)                                                                                                                            ccAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCgcagca                                <  2:151083/136‑1 (MQ=255)                                                                                                                                gCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGtctt               >  1:4526/1‑149 (MQ=255)                                                                                                                                  acaAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCGGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGtcttct             >  7:353588/1‑149 (MQ=255)                                                                                                                                   gcaaaaaCGCCCCTCGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGCCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTtctc            >  4:87380/3‑149 (MQ=255)                                                                                                                                    aaaaaaCGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCAt                               >  6:69463/1‑129 (MQ=255)                                                                                                                                    aaaaaaCGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCAt                               <  5:69463/129‑1 (MQ=255)                                                                                                                                     aaaaaCGCCCCTCGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTcgc          >  7:2445/1‑149 (MQ=255)                                                                                                                                      aaaaCGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCCCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTcgcg         <  3:234040/149‑1 (MQ=255)                                                                                                                                       aaaCGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGa        <  5:67665/149‑1 (MQ=255)                                                                                                                                            ccccGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCgg   >  8:385047/1‑149 (MQ=255)                                                                                                                                             cccGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTCCTTCGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCgg   >  4:300770/1‑148 (MQ=255)                                                                                                                                              ccGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCggg  >  4:291883/1‑148 (MQ=255)                                                                                                                                               cGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCggg  >  1:29374/1‑147 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGC‑GTTTTTGCGTTT‑CAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGG  >  NC_002947/867911‑868196

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