Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,305,604	G→A	6.1%	R151C (CGC→TGC)	PP_3777 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,305,604	0	G	A	6.1%	247.3 / 3.9	82	R151C (CGC→TGC)	PP_3777	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/40);  new base A (2/3);  total (39/43)
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.26e-01

ACCGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGACC  >  NC_002947/4305458‑4305749                                                                                                                                                    |                                                                                                                                                  acCGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTCCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTCGGCAGCCTC‑GATGCGCt                                                                                                                                                 >  8:355809/1‑149 (MQ=255)           ccGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATc                                                                                                                                       <  2:266472/149‑1 (MQ=255)           ccGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATc                                                                                                                                       <  2:261457/149‑1 (MQ=255)             gACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATc                                                                                                                                       >  4:356347/1‑147 (MQ=255)              aCCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAg                                                                                                                                    <  4:156676/149‑1 (MQ=255)              aCCAACGAATTCAGCCCGATGAGCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAg                                                                                                                                    <  4:68862/149‑1 (MQ=255)              aCCAACGAATTCAGCCCGATGAACAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCACTACAGCGCATCGAg                                                                                                                                    <  3:37914/149‑1 (MQ=255)                cAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAgg                                                                                                                                   >  5:375336/1‑148 (MQ=255)                cAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAgg                                                                                                                                   >  8:82904/1‑148 (MQ=255)                cAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAgg                                                                                                                                   >  1:260872/1‑148 (MQ=255)                       ttCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATc                                                                                                                           >  4:336346/1‑149 (MQ=255)                            cccGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAgg                                                                                                                      <  2:260872/149‑1 (MQ=255)                            cccGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAgg                                                                                                                      >  5:339454/1‑149 (MQ=255)                              cGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGc                                                                                                                    <  2:100474/149‑1 (MQ=255)                                   atCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCt                                                                                                                <  2:197850/148‑1 (MQ=255)                                        cTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCAcc                                                                                                          >  7:1773/1‑149 (MQ=255)                                         ttGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCAccc                                                                                                         >  6:375468/1‑149 (MQ=255)                                            ccTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAg                                                                                                                       >  7:215961/1‑132 (MQ=255)                                            ccTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAg                                                                                                                       <  8:215961/132‑1 (MQ=255)                                             cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTCGc                                                                                                     >  2:20432/1‑149 (MQ=255)                                              tCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCTTCCTC‑GATGCACTACAGCGCATCGa                                                                                                                                     <  7:10710/116‑1 (MQ=38)                                              tCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCTTCCTC‑GATGCACTACAGCGCATCGa                                                                                                                                     >  8:10710/1‑116 (MQ=38)                                              tCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCg                                                                                                    >  1:73202/1‑149 (MQ=255)                                                cGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCg                                                                                                    <  8:194261/147‑1 (MQ=255)                                                cGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGa                                                                                                  >  6:50084/1‑149 (MQ=255)                                                 gggCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAt                                                                                                 <  6:375336/149‑1 (MQ=255)                                                  ggCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                                <  3:356347/149‑1 (MQ=255)                                                  ggCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                                <  3:336346/149‑1 (MQ=255)                                                    cTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGa                                                                                              <  5:125242/149‑1 (MQ=255)                                                     ttCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                             >  1:345996/1‑149 (MQ=255)                                                       cAGGTATCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                             <  5:50084/147‑1 (MQ=255)                                                        aGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATa                                                                                           <  6:289166/148‑1 (MQ=255)                                                        aGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAg                                                                                          <  4:346562/149‑1 (MQ=255)                                                          gTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGa                                                                                        <  2:73202/149‑1 (MQ=255)                                                          gTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCACTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGa                                                                                        >  7:4339/1‑149 (MQ=255)                                                               tCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATcac                                                                                    >  4:190039/1‑148 (MQ=255)                                                               tCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATc                                                                                      >  5:247307/1‑146 (MQ=255)                                                               tCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGAGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATc                                                                                      <  6:247307/146‑1 (MQ=255)                                                                 ggCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACAt                                                                                  <  5:5091/148‑1 (MQ=255)                                                                 ggCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATc                                                                                 <  5:375468/149‑1 (MQ=255)                                                                      ccGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATc                                                                            >  3:46561/1‑149 (MQ=255)                                                                       cGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACTTAGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGTAAGCCATCCGATGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATcc                                                                           <  1:20432/149‑1 (MQ=255)                                                                            gAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCAGATGCGGTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACAGCGTATCCGTcc                                                                       >  1:295636/1‑149 (MQ=255)                                                                              aaaaCCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTgg                                                                                                      >  1:81338/1‑115 (MQ=255)                                                                              aaaaCCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTgg                                                                                                      <  2:81338/115‑1 (MQ=255)                                                                                 aCCTGCGCAAGGTCCGCTCGCCATAGAACATCACTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGcc                                                                 <  4:243008/149‑1 (MQ=255)                                                                                   cTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCg                                                                >  8:86076/1‑148 (MQ=255)                                                                                     gcgcAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTc                                                                        >  5:223339/1‑138 (MQ=255)                                                                                     gcgcAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTc                                                                        <  6:223339/138‑1 (MQ=255)                                                                                     gcgcAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAg                                                             >  8:251345/1‑149 (MQ=255)                                                                                       gcAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCATGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAgg                                                            >  8:285155/1‑148 (MQ=255)                                                                                         aaGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAt                                                         <  6:339454/149‑1 (MQ=255)                                                                                          aGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                        >  1:6599/1‑149 (MQ=255)                                                                                                    ggcATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGGGGTGTAGTGCATCGAGGAAGCCATGCGAGGGCGTCTTCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACa                                              <  7:355809/147‑1 (MQ=255)                                                                                                     ccATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCACTACAGCGCATCGAgg                                                                                                                                   <  4:212903/63‑1 (MQ=25)                                                                                                     ccATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCACTACAGCGCATCGAgg                                                                                                                                   >  3:212903/1‑63 (MQ=25)                                                                                                      cATAGAACATCGCTTCCAGCGACAGGCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGa                                            <  8:4339/149‑1 (MQ=255)                                                                                                           aaCATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAg                                                                                          >  5:4784/1‑98 (MQ=255)                                                                                                           aaCATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAg                                                                                          <  6:4784/98‑1 (MQ=255)                                                                                                           aaCATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGc                                                      >  1:189713/1‑134 (MQ=255)                                                                                                           aaCATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGc                                                      <  2:189713/134‑1 (MQ=255)                                                                                                             cATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCt                                     <  4:46561/149‑1 (MQ=255)                                                                                                               tCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGt                                   <  5:340824/149‑1 (MQ=255)                                                                                                               tCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGt                                   <  8:1773/149‑1 (MQ=255)                                                                                                                 gCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTg                                                       >  7:197751/1‑127 (MQ=255)                                                                                                                 gCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTg                                                       <  8:197751/127‑1 (MQ=255)                                                                                                                             aGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACa                     <  3:190039/149‑1 (MQ=255)                                                                                                                             aGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACa                     <  7:82904/149‑1 (MQ=255)                                                                                                                               cccTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTATGCCACCATGGCGGATTGAATAGGATAACAACTTATCCTGCAAATGCCGCAGGCATTGCTTCTAACAGATAACCATGTACGAATTCGATAAAg                    >  2:102850/1‑148 (MQ=255)                                                                                                                                ccTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACATCTCTCCCTGTACGAATTCGAGACAGGa                  >  4:270608/1‑149 (MQ=255)                                                                                                                                 cTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCa                                                              <  1:18975/104‑1 (MQ=255)                                                                                                                                 cTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCa                                                              >  2:18975/1‑104 (MQ=255)                                                                                                                                      gCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGcc            >  4:109008/1‑149 (MQ=255)                                                                                                                                         gCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAg         >  6:379483/1‑149 (MQ=255)                                                                                                                                          ccTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAgg        >  5:175289/1‑149 (MQ=255)                                                                                                                                          ccTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAgg        <  3:109008/149‑1 (MQ=255)                                                                                                                                           cTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGt       >  5:193372/1‑149 (MQ=255)                                                                                                                                            tC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGtt      >  8:302019/1‑149 (MQ=255)                                                                                                                                               gATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGtt      <  7:285155/147‑1 (MQ=255)                                                                                                                                               gATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGa    <  6:175289/149‑1 (MQ=255)                                                                                                                                                 tGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGAcc  <  7:251345/149‑1 (MQ=255)                                                                                                                                                 tGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGAc   >  1:100029/1‑148 (MQ=255)                                                                                                                                                  gcgcTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGTGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGAcc  <  7:86076/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  ACCGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGACC  >  NC_002947/4305458‑4305749

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