Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,605,501	C→G	9.1%	A306G (GCC→GGC)	PP_4077 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,605,501	0	C	G	9.1%	222.1 / 12.7	77	A306G (GCC→GGC)	PP_4077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (37/33);  new base G (2/5);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

CCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATC  >  NC_002947/4605353‑4605648                                                                                                                                                       |                                                                                                                                                     ccGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGc                                                                                                                                                      <  7:79634/149‑1 (MQ=255)    gCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAg                                                                                                                                                   >  7:18949/1‑149 (MQ=255)       caAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGa                                                                                                                                                <  6:370012/149‑1 (MQ=255)                ccACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGcc                                                                                                                                         >  7:117377/1‑147 (MQ=255)                ccACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCa                                                                                                                                       >  4:257173/1‑149 (MQ=255)                 cacaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAg                                                                                                                                      <  7:53676/149‑1 (MQ=255)                  acaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACGCTTCGCGTTGGCACAGTAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCTGCGATGCTGCGCAgg                                                                                                                                     <  4:36956/149‑1 (MQ=255)                   caGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAGGt                                                                                                                                    <  7:136649/149‑1 (MQ=255)                   caGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAGGt                                                                                                                                    <  8:18949/149‑1 (MQ=255)                       tGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAg                                                                                                                                <  3:195540/149‑1 (MQ=255)                         ggCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAGGTCGatg                                                                                                                               <  4:67136/148‑3 (MQ=255)                         ggCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCt                                                                                                                              >  3:368211/1‑149 (MQ=255)                         ggCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTGCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCt                                                                                                                              <  3:257173/149‑1 (MQ=255)                               gAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGcc                                                                                                                        >  6:224659/1‑149 (MQ=255)                               gAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGc                                                                                                                         <  1:242395/148‑1 (MQ=255)                                aGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGccc                                                                                                                       <  4:368211/149‑1 (MQ=255)                                  ttCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGcc                                                                                                                        >  3:368334/1‑146 (MQ=255)                                  ttCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATCTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAGGTCGat                                                                                                                                <  4:116580/138‑2 (MQ=255)                                  ttCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATATCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAGGTCGat                                                                                                                                >  3:116580/1‑137 (MQ=255)                                           aCGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGTTGATGACAACATGTTCAAGCACAGTCCCGAGCCCAACCACCTGGCCAGCGAGGCTACCAAGTTCGAGCTGAAGACCAGAACCGAg                                                                                                              >  1:299065/1‑149 (MQ=255)                                              cgTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCa                                                                                                         <  2:79504/149‑1 (MQ=255)                                               gTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAg                                                                                                        >  1:252130/1‑149 (MQ=255)                                                  aaCTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCAACCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCGCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTcc                                                                                                     >  7:337980/1‑149 (MQ=255)                                                       ccAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCt                                                                                                  <  6:44014/147‑1 (MQ=255)                                                          gACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCCGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACg                                                                                              >  8:115279/1‑148 (MQ=255)                                                           acacCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGa                                                                                             >  8:2776/1‑148 (MQ=255)                                                              cccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑c                                                                                        >  3:62568/1‑149 (MQ=255)                                                               ccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑Ct                                                                                       >  5:271327/1‑149 (MQ=255)                                                               ccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑Ct                                                                                       >  2:205097/1‑149 (MQ=255)                                                                cTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCACGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑Ctt                                                                                      >  7:100231/1‑149 (MQ=255)                                                                 ttCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑Ctt                                                                                      >  3:28000/1‑148 (MQ=255)                                                                    gcgTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAg                                                                                   <  8:66452/148‑1 (MQ=255)                                                                      gTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAgtgt                                                                                >  1:53494/1‑149 (MQ=255)                                                                       ttGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTc                                                                               <  2:53494/149‑1 (MQ=255)                                                                          gCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAAGCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGTAAT‑CTTCAGTGTCGa                                                                             <  2:299065/148‑1 (MQ=255)                                                                           cacaGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGTCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATg                                                                           <  8:141802/149‑1 (MQ=255)                                                                             caGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGTGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGaa                                                                         <  3:147453/149‑1 (MQ=255)                                                                                aagaTTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCAAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGt                                                                       >  5:331643/1‑148 (MQ=255)                                                                                 agaTTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGtc                                                                      <  4:62568/148‑1 (MQ=255)                                                                                    tttCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGtcatc                                                                   <  8:117377/148‑1 (MQ=255)                                                                                    cttCAGCCTGTTCTGCTCCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCATGTCGAGCTGAAGCCCAGAAGCCAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCTATGAATtcatc                                                                   <  8:337980/147‑1 (MQ=255)                                                                                          ccTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCAccc                                                            <  4:28000/149‑1 (MQ=255)                                                                                             gTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGaa                                                         >  8:54080/1‑149 (MQ=255)                                                                                                cTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGaaca                                                       <  4:368334/148‑1 (MQ=255)                                                                                                 tGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGaacaac                                                     >  4:310633/1‑149 (MQ=255)                                                                                                   cACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAAcc                                                    >  7:256245/1‑148 (MQ=255)                                                                                                      cccG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACt                                                 >  4:104517/1‑148 (MQ=255)                                                                                                            gCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACgg                                            <  3:104517/149‑1 (MQ=255)                                                                                                                tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCGCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAg                                         <  6:331643/148‑1 (MQ=255)                                                                                                                tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGtcat                                                                    <  8:162268/121‑1 (MQ=255)                                                                                                                tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGtcat                                                                    >  7:162268/1‑121 (MQ=255)                                                                                                                 cAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGATCTGAAGCCAAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CATCAGTGACGATGAAGTCATCAGAACCCCAACAACAACTGACCGAAGTa                                       >  4:33299/1‑149 (MQ=255)                                                                                                                  aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTAc                                      >  2:249110/1‑149 (MQ=255)                                                                                                                  aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTAc                                      <  7:2776/149‑1 (MQ=255)                                                                                                                  aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTACCCAGGTCGAGCTGAAGCCCAGAAGCGAAGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTaa                                      >  1:86488/1‑148 (MQ=255)                                                                                                                       gTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACg                                             <  4:147940/137‑1 (MQ=255)                                                                                                                       gTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACg                                             >  3:147940/1‑137 (MQ=255)                                                                                                                        ttCAATCACAGTGCCGAGCCGATGGACCTGGCCAGCGAGGCTGGCCAGGTCGAGCTGAATCCCATAAGCGAGGCGCAGTGCGGTTACGAAATCCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGt                                  <  1:102537/148‑1 (MQ=255)                                                                                                                        ttCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACACCCACTGACGGAAGTACCGGt                                  >  3:113674/1‑147 (MQ=255)                                                                                                                        ttCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGt                                  <  7:21671/147‑1 (MQ=255)                                                                                                                        ttCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGt                                  <  3:33299/147‑1 (MQ=255)                                                                                                                        ttCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGa                                <  6:271327/149‑1 (MQ=255)                                                                                                                            aGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGACGCCCAGAAGCGAGGCGCAGTCCGCTTACGAACT‑CGTCAGTGTCGATGAAGGCATCAGCACCCGAACAACCACGGACGGCAGTACGGGTGAGCAc                            >  8:130065/1‑149 (MQ=255)                                                                                                                               acaGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTg                         >  7:352400/1‑149 (MQ=255)                                                                                                                                   tGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGa                     >  2:47724/1‑149 (MQ=255)                                                                                                                                    gCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGAc                    >  5:57350/1‑149 (MQ=255)                                                                                                                                         gCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTc               <  2:252130/149‑1 (MQ=255)                                                                                                                                          cccATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCg              <  5:224659/149‑1 (MQ=255)                                                                                                                                           ccAGCGACCTGGCGAGCGAGGCGGACCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGAAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCTAACAACCACTGACGTAAGTACCGGTGAGCACCTGCGGACTTTTCGc             <  8:352400/149‑1 (MQ=38)                                                                                                                                                gACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCttt        <  7:115279/149‑1 (MQ=255)                                                                                                                                                 aCCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCtt         <  6:45245/147‑1 (MQ=255)                                                                                                                                                 aCCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCtt         >  5:45245/1‑147 (MQ=255)                                                                                                                                                 aCCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTc       >  3:122742/1‑149 (MQ=255)                                                                                                                                                 aCCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTc       >  7:310249/1‑149 (MQ=255)                                                                                                                                                 aCCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTc       >  2:257061/1‑149 (MQ=255)                                                                                                                                                   cTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGACCTTTCg      >  1:134672/1‑148 (MQ=255)                                                                                                                                                      gCCATCGAGGCTGACAAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCTCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGaa    <  8:256245/147‑1 (MQ=255)                                                                                                                                                      gCCAGCGAGGCTGCCCAGGACGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATc  >  8:7845/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                     CCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCG‑‑GCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAAT‑CTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATC  >  NC_002947/4605353‑4605648

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