Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,200,096	G→T	6.5%	S314S (TCC→TCA)	PP_3690 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,200,096	0	G	T	6.5%	188.6 / 2.8	61	S314S (TCC→TCA)	PP_3690	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/31);  new base T (2/2);  total (28/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.37e-01

CCTTGGCTGACTTTTCCTGTGTCTGAGAGAACACGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAAC  >  NC_002947/4199953‑4200240                                                                                                                                                |                                                                                                                                                   ccTTGGCTGACTTTTCCTGTGTCTGAGAGAACACGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAAt                                                                                                                                               >  6:257093/1‑149 (MQ=255)  cTTGGCTGACTTTTCCTGTGTCTGAGAGAACACGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGATTGGAAAtt                                                                                                                                              >  3:93444/1‑149 (MQ=255)  cTTGGCTGACTTTTCCTGTGTCTGAGAGAACACGAAGAGCTGATTTGGGCTTGAAGTTGCCGGTGCCACAGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTAGCTTGGAAAtt                                                                                                                                              >  6:275380/1‑149 (MQ=255)        tGACTTTTCCTGTGTCTGAGAGAACACGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGAACCGCTTGCTTGGAAATTGCGGc                                                                                                                                         >  7:345731/1‑148 (MQ=255)              ttCCTGTGTCTGAGAGAACACGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATc                                                                                                                                   <  5:257093/148‑1 (MQ=255)                ccTGTGTCTGAGAGAACACGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGtt                                                                                                                                >  6:76412/1‑149 (MQ=255)                 cTGTGTCTGAGAGAACACGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGtt                                                                                                                                <  1:94855/148‑1 (MQ=255)                 cTGTGTCTGAGAGAACACGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGtt                                                                                                                                <  2:318710/148‑1 (MQ=255)                         gagagaACACGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTcc                                                                                                                       >  7:18729/1‑149 (MQ=255)                          agagaACACGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCTAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTc                                                                                                                        <  2:375998/147‑1 (MQ=255)                            agaACACGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTcc                                                                                                                       >  1:97666/1‑146 (MQ=255)                               acacGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATc                                                                                                                 >  7:295425/1‑149 (MQ=255)                                 acGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCt                                                                                                                >  4:263450/1‑148 (MQ=255)                                  cGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTc                                                                                                              <  5:76412/149‑1 (MQ=255)                                     agagCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCaaa                                                                                                           <  4:13141/149‑1 (MQ=255)                                      gagTTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAg                                                                                                          <  2:211288/149‑1 (MQ=255)                                         cTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATc                                                                                                                                   <  4:217158/121‑1 (MQ=255)                                         cTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATc                                                                                                                                   >  3:217158/1‑121 (MQ=255)                                         cTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCAACAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTAc                                                                                                       <  6:153694/149‑1 (MQ=255)                                            ttttGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGcc                                                                                                    >  8:117479/1‑149 (MQ=255)                                                gggCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTTCGTTCCGAGATCTTCAAAGTACGCCTTc                                                                                                 <  1:209555/148‑1 (MQ=255)                                                gggCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTTCGTTCCGAGATCTTCAAAGTACGCCTTc                                                                                                 <  1:215357/148‑1 (MQ=255)                                                 ggCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCa                                                                                                <  6:229058/148‑1 (MQ=255)                                                        aGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGccc                                                                                        <  3:263450/149‑1 (MQ=255)                                                          ttGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTc                                                                                                                        <  6:222694/115‑1 (MQ=255)                                                             ccGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAATTACGCCTTCAGAACGCCCTCGat                                                                                   <  6:66639/149‑1 (MQ=255)                                                             ccGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGa                                                                                    <  4:46768/148‑1 (MQ=255)                                                             ccGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCATCAAAATACGCCTTCAGAACGCACTCGa                                                                                    >  8:230973/1‑148 (MQ=255)                                                                     aCCGAGGCAGGAGCGGATTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGc                                                                         <  1:205158/149‑1 (MQ=255)                                                                     aCCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGc                                                                         >  1:130169/1‑149 (MQ=255)                                                                     aCCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGc                                                                         >  8:289692/1‑149 (MQ=255)                                                                      ccGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGat                                                                                   >  8:253828/1‑140 (MQ=255)                                                                      ccGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGat                                                                                   <  7:253828/140‑1 (MQ=255)                                                                      ccGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGc                                                                         >  8:176025/1‑148 (MQ=255)                                                                      ccGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCt                                                                        <  2:130169/149‑1 (MQ=255)                                                                       cGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCt                                                                        >  3:227661/1‑148 (MQ=255)                                                                            cAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGtt                                                                  <  1:10917/149‑1 (MQ=255)                                                                                   ggAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGtgct                                                           <  1:216613/149‑1 (MQ=255)                                                                                    gAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑GTACGCTAACGTTCAGtgctg                                                          >  6:320069/1‑149 (MQ=255)                                                                                    gAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGtgctg                                                          >  3:19643/1‑149 (MQ=255)                                                                                      gTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGcc                                                        <  7:117479/149‑1 (MQ=255)                                                                                      gTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑AAACGCTAACGTTCAGTGCTGcc                                                        >  7:355404/1‑149 (MQ=255)                                                                                               ttACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTAt                                                <  6:127663/148‑1 (MQ=255)                                                                                                 aCGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGt                                                              <  8:32293/132‑1 (MQ=255)                                                                                                 aCGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGt                                                              >  7:32293/1‑132 (MQ=255)                                                                                                   gcgcTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTAtgcctg                                           >  4:203418/1‑149 (MQ=255)                                                                                                   gcgcTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCCTTGTCCGTTCCCAGATCTTCAAAGTACGCCTTCAGAACGCCCACGATAC‑‑ATACGCTACCGTTCAGAGCTGCCGAATCTAtgccag                                           >  7:64150/1‑147 (MQ=255)                                                                                                         atcaACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTTGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTAtgcctgcctgc                                      >  5:155727/1‑148 (MQ=255)                                                                                                         atcaACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTTGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTAtgcctgcctgc                                      >  5:149900/1‑148 (MQ=255)                                                                                                            aaCCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAGCGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAg                                  <  7:289692/149‑1 (MQ=255)                                                                                                             aCCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGc                                 >  3:51296/1‑149 (MQ=255)                                                                                                                   aGCTGCATCACGTCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTAtt                           <  7:230973/149‑1 (MQ=255)                                                                                                                          tCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTg                    >  1:31784/1‑149 (MQ=255)                                                                                                                            aCGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCAATAGAGCTTggc                  <  8:345731/149‑3 (MQ=255)                                                                                                                             cGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTc                 >  2:41906/1‑149 (MQ=255)                                                                                                                              ggCTGGCACCGCTTGCTTTGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCAATAGAGCTTggcag                <  6:155727/149‑5 (MQ=255)                                                                                                                              ggCTGGCACCGCTTGCTTTGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCAATAGAGCTTggcag                <  6:149900/149‑5 (MQ=255)                                                                                                                                cTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAg              <  8:18729/149‑1 (MQ=255)                                                                                                                                   gCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCAGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGc             <  7:73547/147‑1 (MQ=255)                                                                                                                                   gCACCGCTTGCTTGGAAATTGCCGCAAAATTGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATACTAATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGc             <  1:302047/149‑1 (MQ=255)                                                                                                                                    cACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGt          <  4:227661/149‑1 (MQ=255)                                                                                                                                             cttgGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAac  <  2:31784/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   CCTTGGCTGACTTTTCCTGTGTCTGAGAGAACACGAAGAGCTGTTTTGGGCTTGAAGTTGCCGGTGCCACCGAGGCAGGAGCGGAGTAAATTTTTTACGCGCTCATCAACCTCAAGCTGCATCACGGCTGGCACCGCTTGCTTGGAAATTGCGGCAAAATCGTTGTCCGTTCCGAGATCTTCAAAGTACGCCTTCAGAACGCCCTCGATAC‑‑ATACGCTAACGTTCAGTGCTGCCGAATCTATGCCTGCCTGCCAAGCTCTATTGAGCTTGTTCCAGCTGTCATAAAAC  >  NC_002947/4199953‑4200240

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