Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,384,784	A→C	7.6%	P211P (CCT→CCG)	PP_2986 ←	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,384,784	0	A	C	7.6%	197.5 / 7.3	66	P211P (CCT→CCG)	PP_2986	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base C (2/3);  total (31/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 = 8.12e-01

CGCCCGCCATACGTCCTAAAACAGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGGCATGCATCTGGGGGTGGCTCTCCAGGAACGACAATAC  >  NC_002947/3384636‑3384930                                                                                                                                                      |                                                                                                                                                   cgcccgccATACGTCCTAAAACAGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCTAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCa                                                                                                                                                    <  5:142868/149‑1 (MQ=255)  gcccgccATACGTCCTAAAACAGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAg                                                                                                                                                   >  7:178004/1‑149 (MQ=255)   cccgccATACGTCCTAAAACAGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAgg                                                                                                                                                  >  6:229939/1‑149 (MQ=255)    ccgccATACGTCCTAAAACAGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCG‑CCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGa                                                                                                                                                <  2:162651/149‑1 (MQ=255)                  aaaaCAGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCa                                                                                                                                    >  6:4360/1‑148 (MQ=255)                   aaaCAGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGTCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCCGGGTTGATCCCTGGCagt                                                                                                                                  <  5:5435/149‑3 (MQ=255)                   aaaCAGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCaa                                                                                                                                   <  5:229939/148‑1 (MQ=255)                   aaaCAGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAAc                                                                                                                                  <  2:279099/149‑1 (MQ=255)                      cAGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGCTGGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACtt                                                                                                                                <  1:262126/149‑1 (MQ=255)                      cAGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTg                                                                                                                               >  3:94626/1‑149 (MQ=255)                       aGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAg                                                                                                                                                   >  5:190547/1‑128 (MQ=255)                       aGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAg                                                                                                                                                   <  6:190547/128‑1 (MQ=255)                        ggCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGcc                                                                                                                             <  1:20254/149‑1 (MQ=255)                        ggCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCAAGGGATCAACCCGGGCAACTTGcc                                                                                                                             <  6:167080/149‑1 (MQ=255)                           aTCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAg                                                                                                                                                   <  4:350143/124‑1 (MQ=255)                           aTCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAg                                                                                                                                                   >  3:350143/1‑124 (MQ=255)                                       aTTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAg                                                                                                                                                   >  1:65294/1‑112 (MQ=255)                                       aTTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAg                                                                                                                                                   <  2:65294/112‑1 (MQ=255)                                       aTTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGtt                                                                                                              <  8:178004/149‑1 (MQ=255)                                       aTTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGtt                                                                                                              <  7:66522/149‑1 (MQ=255)                                          gTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGa                                                                                                            <  4:99580/148‑1 (MQ=255)                                           tCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCCgg                                                                                                                                                  >  1:319753/1‑109 (MQ=255)                                           tCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCCgg                                                                                                                                                  <  2:319753/109‑1 (MQ=255)                                           tCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGAc                                                                                                           >  1:113351/1‑148 (MQ=255)                                           tCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGAc                                                                                                           >  6:322403/1‑148 (MQ=255)                                             atatGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCaa                                                                                                        >  3:171784/1‑149 (MQ=255)                                             atatGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCaa                                                                                                        >  7:313467/1‑149 (MQ=255)                                                   cGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCCGGGTTGATCCCTGGCagttt                                                                                                                                <  5:305730/119‑5 (MQ=255)                                                   cGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCCGGGTTGATCCCTGGCagttt                                                                                                                                >  6:305730/1‑115 (MQ=255)                                                        cgccgcGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGc                                                                                                                                     >  7:145388/1‑109 (MQ=255)                                                        cgccgcGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGc                                                                                                                                     <  8:145388/109‑1 (MQ=255)                                                            gcgAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGt                                                                                         <  5:282771/149‑1 (MQ=255)                                                            gcgAGATAGAGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAgtt                                                                                         <  4:94626/149‑3 (MQ=255)                                                                        tCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACcac                                                                                                                    <  1:260019/110‑1 (MQ=255)                                                                        tCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACcac                                                                                                                    >  2:260019/1‑110 (MQ=255)                                                                        ccacagcGGCTCCCCTGTCTGGCTGCGGCAATGTGTGCCAAGGCCACATCTC‑GGGCATCTACAATGGAAAAACTGACATGGATCAACCCGGGCAACTTGCCCTGTACCACAGCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTc                                                                             <  8:248120/143‑1 (MQ=255)                                                                            cTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGcc                                                                          >  5:80999/1‑148 (MQ=255)                                                                            cTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGcc                                                                          <  6:80999/148‑1 (MQ=255)                                                                            cTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGc                                                                           >  1:287383/1‑147 (MQ=255)                                                                            cTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGc                                                                           >  1:143738/1‑147 (MQ=255)                                                                             tCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAg                                                                        >  2:350749/1‑149 (MQ=255)                                                                               ggCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGc                                                                      <  4:171784/149‑1 (MQ=255)                                                                                gcgcCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGcc                                                                     <  2:287383/149‑1 (MQ=255)                                                                                     cATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCca                                                                 >  2:52967/1‑148 (MQ=255)                                                                                         tCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATc                                                             >  3:148205/1‑148 (MQ=255)                                                                                         tCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATc                                                             >  3:142779/1‑148 (MQ=255)                                                                                                     aaTGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGGCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGa                                                <  5:322403/149‑1 (MQ=255)                                                                                                               aaTGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCATGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGgc                                      <  8:246396/149‑1 (MQ=255)                                                                                                               aaTGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGCCCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCCCATCCCGCCAGGCAGGACCATACAgg                                       >  7:160069/1‑148 (MQ=255)                                                                                                               aaTGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGgc                                      <  4:142779/149‑1 (MQ=255)                                                                                                               aaTGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGgc                                      <  4:148205/149‑1 (MQ=255)                                                                                                                 tGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGgcat                                    <  5:4360/149‑1 (MQ=255)                                                                                                                    cacaTCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCACCCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCCGGCAGGACCATACAGgcatgca                                 >  1:80656/1‑149 (MQ=255)                                                                                                                    cacaTCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTg                                                                                                     >  1:339649/1‑81 (MQ=255)                                                                                                                    cacaTCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTg                                                                                                     <  2:339649/81‑1 (MQ=255)                                                                                                                     acaTCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAAACCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGgcatgcat                                >  6:343413/1‑149 (MQ=255)                                                                                                                        tCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGGCATGCATCTg                             >  5:70102/1‑149 (MQ=255)                                                                                                                                     aCAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGGCATGCATCTGGGGGTGGCtctc                 <  6:70102/149‑1 (MQ=255)                                                                                                                                        aTGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGAc                                               >  1:142319/1‑116 (MQ=255)                                                                                                                                        aTGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGAc                                               <  2:142319/116‑1 (MQ=255)                                                                                                                                         tGGAAAAACTTCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGgcatg                                   <  6:89032/127‑1 (MQ=255)                                                                                                                                           taaaaaCTGCCTGGGATCAACCGGGCCAACTTGCCCTGTACCACATCGTTGACCAATTCCCCCGACGAGGTTGGGCCAAGGTCGCAAGGCCCCCACATCCAGCCAGGCTGGACAATACAGGCATGCATCTGGGGGTGGCTCTCCAGGaa           <  1:52967/148‑1 (MQ=255)                                                                                                                                             aaaaCTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGGCATGCATCTGGGGGTGGCTCTCCAGGAAc          >  3:19590/1‑148 (MQ=255)                                                                                                                                              aaaCTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGGCATGCATCTGGGGGTGGCTCTCCAGGAACGa        >  1:91816/1‑149 (MQ=255)                                                                                                                                                aCTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGGCATGCATCTGGGGGTGGCTCTCCAGGAACGACa      <  1:350749/149‑1 (MQ=255)                                                                                                                                                  tGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGGCATGCATCTGGGGGTGGCTCTCCAGGAACGACaa     >  5:164864/1‑148 (MQ=255)                                                                                                                                                    ccAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGGCATGCATCTGGGGGTGGCTCTCCAGGAACGACAATAc  <  4:19590/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CGCCCGCCATACGTCCTAAAACAGGCATCAGTTCGCGCATTGTCATATGTCGACCCGCCGCGAGATAGCGCTCACCTCGGCGCCCATGTCTGGCTGCGGCAATGTGTGCCAATGCCACATCGC‑GGGCATCGACAATGGAAAAACTGCCAGGGATCAACCCGGGCAACTTGCCCTGTACCACATCGTTGACCAATTGCCCCGACGAGGTTGGGCCAACGTCGCCAGGCCCCCACATCCAGCCAGGCAGGACCATACAGGCATGCATCTGGGGGTGGCTCTCCAGGAACGACAATAC  >  NC_002947/3384636‑3384930

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