Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	26,245	A→T	6.2%	Q122L (CAG→CTG)	PP_0021 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	26,245	0	A	T	6.2%	240.4 / 3.6	80	Q122L (CAG→CTG)	PP_0021	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (36/39);  new base T (3/2);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

CTAAGAGCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGATCTTTGCT  >  NC_002947/26101‑26393                                                                                                                                                 |                                                                                                                                                      ctaAGAGCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGtt                                                                                                                                                    <  2:312881/148‑1 (MQ=255) ctaAGAGCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTc                                                                                                                                                   >  3:300448/1‑149 (MQ=255)    agagCATTCGCGGGATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCCTACTTCCAGCCTCGCCAGCTCAAGAACACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCTAAGAGTTCGAt                                                                                                                                                >  8:353458/1‑149 (MQ=255)    agagCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGAt                                                                                                                                                >  1:141098/1‑149 (MQ=255)     gagCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATc                                                                                                                                               >  3:212738/1‑149 (MQ=255)      agCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATcg                                                                                                                                              >  8:197448/1‑149 (MQ=255)       gCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATcgc                                                                                                                                             <  3:207507/149‑1 (MQ=255)       gCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATcgc                                                                                                                                             >  6:94574/1‑149 (MQ=255)        cATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATcgcg                                                                                                                                            >  1:367514/1‑149 (MQ=255)         aTTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACTGTTCGATcaata                                                                                                                                           >  3:173985/1‑145 (MQ=255)         aTTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACTGTTCGATcaata                                                                                                                                           <  4:173985/149‑5 (MQ=255)          ttagcgGCATCGCTACAATTGATACTGAGGCCCGTAAAGATTTCCCGTACTTCCATCCTCGCCATCAAAATAAAACATAACGAAAAACCATGTTTGACCGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTg                                                                                                                                          <  3:142648/146‑1 (MQ=255)          ttCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACTGTTCGATcaatag                                                                                                                                          <  8:152990/149‑6 (MQ=255)            cgcgGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTgg                                                                                                                                         <  8:282266/148‑1 (MQ=255)             gcgGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGtt                                                                                                                                       >  2:57619/1‑149 (MQ=255)               ggCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑Ta                                                                                                                                    <  4:180364/149‑1 (MQ=255)                      cTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTg                                                                                                                             <  1:149379/149‑1 (MQ=255)                      cTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTg                                                                                                                             <  7:190079/149‑1 (MQ=255)                        ccATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTACCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTg                                                                                                                             >  5:121050/1‑147 (MQ=255)                         cATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTgg                                                                                                                            <  5:104099/147‑1 (MQ=255)                          aattGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGgaga                                                                                                                         <  2:158100/147‑1 (MQ=255)                               aTAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAaga                                                                                                                     <  1:240778/148‑1 (MQ=255)                                 aGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGAc                                                                                                                  >  1:150484/1‑149 (MQ=255)                                    gAGGCCCGTAACGATTGCCCGTACTTCCAGCCGCGGCAGCACAAGAAATCAGAACGAAAAACCAGGTTTGACGGCGTCCCACGTAAGATTCTCAAGCTATTGATCGAACAGTTCGTTCGCGTGGTT‑TACATCCTGGAGAAAGAGAccc                                                                                                                <  7:353458/148‑1 (MQ=255)                                     aGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCaa                                                                                                              >  2:279667/1‑149 (MQ=255)                                                gATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTc                                                                                                   >  1:227843/1‑149 (MQ=255)                                                gATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTc                                                                                                   <  1:57619/149‑1 (MQ=255)                                                 aTTGCCCGTACTTCCATCCTCGCCAGCACATGAAAACAGAACGAAAAACGCGGTTTGACGGCGTCCCACTGAAGATTCTCAAGCTATTGATCGACCAGCTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCACCTGGTCCTTTcc                                                                                                  <  2:322174/149‑1 (MQ=255)                                                  ttGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCg                                                                                                 <  2:136112/149‑1 (MQ=255)                                                    gCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCga                                                                                                <  4:127905/148‑1 (MQ=255)                                                      ccGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGaa                                                                                             <  5:147694/149‑1 (MQ=255)                                                            ttCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTAAATCGCGTGGTT‑TACATCCTGGAGAAAGCGACCCAACTGGTCCTTTCCGAGAACGCGta                                                                                       >  6:381445/1‑148 (MQ=255)                                                              ccAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGtt                                                                                       >  3:183956/1‑147 (MQ=255)                                                                aGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTga                                                                                     >  2:5818/1‑147 (MQ=255)                                                                 gCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTCTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTgaga                                                                                   <  1:279667/149‑1 (MQ=255)                                                                  ccTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACa                                                                                                                                  <  2:330557/100‑1 (MQ=255)                                                                  ccTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACa                                                                                                                                  >  1:330557/1‑100 (MQ=255)                                                                   cTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGc                                                                                <  2:367514/149‑1 (MQ=255)                                                                    tCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCa                                                                               >  3:215801/1‑149 (MQ=255)                                                                            acaAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAg                                                                       >  5:43892/1‑149 (MQ=255)                                                                             caAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGc                                                                      <  4:134105/149‑1 (MQ=255)                                                                                  aaaCAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTAc                                                                 <  4:212738/149‑1 (MQ=255)                                                                                     cAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGCCCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAg                                                              >  4:220751/1‑149 (MQ=255)                                                                                     cAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAg                                                              <  3:341512/149‑1 (MQ=255)                                                                                       gAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGTAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAggg                                                            <  6:121050/149‑1 (MQ=255)                                                                                       gAACGAAAAACCAGGTTTGACGGCGTCCCACGCAAGATACTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTTGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGACAGGAAAGCTGCAGCGTTACAAgg                                                             >  7:173312/1‑148 (MQ=255)                                                                                        aaCGTAAAACCATGTTTGACGGCTTACCTCTGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TTCATCCTGGAGAAAGAGACCCAACTGGTCGTTTCCGAGTACGCGTTGAGAGGCATGCTGCAGCGTTACTAGGga                                                           <  8:173312/149‑1 (MQ=255)                                                                                         aCGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGt                                                                    >  3:248286/1‑139 (MQ=255)                                                                                         aCGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGt                                                                    <  4:248286/139‑1 (MQ=255)                                                                                             aaaaCCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACa                                                                >  7:33995/1‑139 (MQ=255)                                                                                             aaaaCCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACa                                                                <  8:33995/139‑1 (MQ=255)                                                                                               acccAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTCCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCACCTGGTCCTTTCCGAGAACGCGGTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTg                                                    <  5:381445/147‑1 (MQ=255)                                                                                               aaCCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACTTCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTg                                                    <  4:300448/149‑1 (MQ=255)                                                                                                    ggTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTAc                                               <  4:73144/149‑1 (MQ=255)                                                                                                    ggTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTAc                                               <  4:183956/149‑1 (MQ=255)                                                                                                      tttGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACct                                             <  5:349662/149‑1 (MQ=255)                                                                                                      tttGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACct                                             <  5:94574/149‑1 (MQ=255)                                                                                                         gACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAGCTGGTCCTTTCCGAGAACGCGGTGAGAGGCATGCTGCCGCGTTAGACGGGAGAGCGTGGCGACctctg                                          >  8:249303/1‑148 (MQ=255)                                                                                                         gACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACctct                                           <  1:5818/148‑1 (MQ=255)                                                                                                            ggCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATAc                                       >  2:20482/1‑149 (MQ=255)                                                                                                              cGTCCCACGGAAGATTCTCAATCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCg                                     <  6:43892/149‑1 (MQ=255)                                                                                                                  ccACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGgcgc                                 >  6:298335/1‑149 (MQ=255)                                                                                                                  ccACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGgcgc                                 <  2:141098/149‑1 (MQ=255)                                                                                                                   cACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGgcgc                                 >  6:25329/1‑148 (MQ=255)                                                                                                                   cACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGgcgc                                 >  4:78507/1‑148 (MQ=255)                                                                                                                        aaGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCt                           >  5:360011/1‑149 (MQ=255)                                                                                                                        aaGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCt                           >  5:306272/1‑149 (MQ=255)                                                                                                                               tcAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAAt                    >  6:30046/1‑149 (MQ=255)                                                                                                                               tcAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAAt                    >  3:113322/1‑149 (MQ=255)                                                                                                                                cAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAAtg                   >  6:68500/1‑149 (MQ=255)                                                                                                                                 aaGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAAtgt                  <  5:68500/149‑1 (MQ=255)                                                                                                                                   gCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCCGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGc                >  3:351907/1‑149 (MQ=255)                                                                                                                                   ccacgcgATCGAACTGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGTCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGc                >  7:243818/7‑149 (MQ=255)                                                                                                                                      aTTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCAt             >  1:154531/1‑149 (MQ=255)                                                                                                                                       gcgATCGAACTGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATg            >  6:113327/3‑149 (MQ=255)                                                                                                                                            cGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGATCt       >  7:252436/1‑149 (MQ=255)                                                                                                                                             gAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGATCtt      >  1:115312/1‑149 (MQ=255)                                                                                                                                             gAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGATCtt      <  8:252436/149‑1 (MQ=255)                                                                                                                                              aaCAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCTCAATGTGCCATGGATCttt     <  6:306272/149‑1 (MQ=255)                                                                                                                                              aaCAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGATCttt     <  2:227843/149‑1 (MQ=255)                                                                                                                                                 aGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGATCTTTGCt  <  3:220751/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                      CTAAGAGCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTT‑TACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGATCTTTGCT  >  NC_002947/26101‑26393

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: ‑