Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,305,611	T→C	5.8%	A148A (GCA→GCG)	PP_3777 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,305,611	0	T	C	5.8%	257.0 / 4.5	86	A148A (GCA→GCG)	PP_3777	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/44);  new base C (2/3);  total (39/47)
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.13e-01

CCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGACCAGCTCATT  >  NC_002947/4305468‑4305757                                                                                                                                                 |                                                                                                                                                   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                                                                                                                                           >  1:260872/1‑148 (MQ=255)      cAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAgg                                                                                                                                           >  5:375336/1‑148 (MQ=255)      cAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAgg                                                                                                                                           >  8:82904/1‑148 (MQ=255)             ttCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATc                                                                                                                                   >  4:336346/1‑149 (MQ=255)                  cccGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAgg                                                                                                                              >  5:339454/1‑149 (MQ=255)                  cccGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAgg                                                                                                                              <  2:260872/149‑1 (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                                                                                                                                             >  8:10710/1‑116 (MQ=38)                                    tCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCTTCCTC‑GATGCACTACAGCGCATCGa                                                                                                                                             <  7:10710/116‑1 (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:336346/149‑1 (MQ=255)                                        ggCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                                        <  3:356347/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                                                                                <  6:223339/138‑1 (MQ=255)                                                                           gcgcAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTc                                                                                >  5:223339/1‑138 (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                                                                                                                                           >  3:212903/1‑63 (MQ=25)                                                                                           ccATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCACTACAGCGCATCGAgg                                                                                                                                           <  4:212903/63‑1 (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                                           <  8:1773/149‑1 (MQ=255)                                                                                                     tCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGt                                           <  5:340824/149‑1 (MQ=255)                                                                                                       gCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTg                                                               <  8:197751/127‑1 (MQ=255)                                                                                                       gCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTg                                                               >  7:197751/1‑127 (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                                                                      >  2:18975/1‑104 (MQ=255)                                                                                                                       cTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCa                                                                      <  1:18975/104‑1 (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)                                                                                                                                            tGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                                >  5:255992/1‑90 (MQ=255)                                                                                                                                            tGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                                <  6:255992/90‑1 (MQ=255)                                                                                                                                              tAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCCTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGACCAGCTCa    <  3:270608/148‑1 (MQ=255)                                                                                                                                                gTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGACCAGCTCAtt  <  2:295636/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTC‑GATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGACCAGCTCATT  >  NC_002947/4305468‑4305757

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