Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	870,231	A→C	7.1%	A64A (GCT→GCG)	PP_5456 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	870,231	0	A	C	7.1%	208.9 / 7.2	70	A64A (GCT→GCG)	PP_5456	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/30);  new base C (3/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

GGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGC  >  NC_002947/870094‑870377                                                                                                                                          |                                                                                                                                                   gggCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCg                                                                                                                                         >  6:151530/1‑149 (MQ=255)   gCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGAGCAAACCCggg                                                                                                                                       >  8:112467/1‑149 (MQ=255)      cGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCgg                                                                                                                                        <  6:32267/145‑1 (MQ=255)      cGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCgg                                                                                                                                        >  5:32267/1‑145 (MQ=255)       gATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTc                                                                                                                                   >  1:181422/1‑149 (MQ=255)       gATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTc                                                                                                                                   >  3:274522/1‑149 (MQ=255)              agaCGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGc                                                                                                                            <  3:349025/149‑1 (MQ=255)                   tCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTc                                                                                                                       >  3:265406/1‑149 (MQ=255)                     cTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTcgc                                                                                                                     >  5:127913/1‑149 (MQ=255)                      ttATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTcgcg                                                                                                                    >  1:64868/1‑149 (MQ=255)                      ttATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTcgcg                                                                                                                    >  3:236324/1‑149 (MQ=255)                              ggggggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccagcc                                                                                                            <  1:155352/149‑1 (MQ=255)                                       ttATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAAccc                                                                                                                                          <  8:29235/110‑1 (MQ=255)                                       ttATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAAccc                                                                                                                                          >  7:29235/1‑110 (MQ=255)                                                          aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCAc                                                                                >  2:258001/1‑149 (MQ=255)                                                          aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCAc                                                                                >  3:307274/1‑149 (MQ=255)                                                          aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCAc                                                                                >  6:266462/1‑149 (MQ=255)                                                          aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTATGTCAGGCAGCGAGGAGGCGCTCGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCAc                                                                                >  5:317623/1‑149 (MQ=255)                                                           gAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCAcc                                                                               >  5:177742/1‑149 (MQ=255)                                                           gAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCAcc                                                                               <  4:307274/149‑1 (MQ=255)                                                               tGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCTGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGt                                                                           <  5:269281/149‑1 (MQ=255)                                                                          gTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCCGCACCCAGTCGATGCCCCgt                                                                >  6:130536/1‑149 (MQ=255)                                                                           tGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtg                                                               <  5:304634/149‑1 (MQ=255)                                                                           tGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCCGCACCCAGTCGATGCCCCgtg                                                               <  5:130536/149‑1 (MQ=255)                                                                           tGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCCGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtg                                                               <  1:300815/149‑1 (MQ=255)                                                                           aggTATTTATTGCACCGTTTGTCAGGCAGCGAGGAGGCGCTGGTCCCAACGGCCGTTTGCGGGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTACGCTGATCTGCACCCACTCGATCCCCCgtg                                                               >  2:237551/2‑149 (MQ=255)                                                                               ttttATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCgg                                                                                                                                        >  2:318704/1‑72 (MQ=255)                                                                               ttttATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCgg                                                                                                                                        <  1:318704/72‑1 (MQ=255)                                                                               ttttATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGc                                                           >  3:54118/1‑149 (MQ=255)                                                                                 ttATTGCAGCGTTTGGCAGGCAGCGAGGCGGCGCGGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGcc                                                          <  4:220495/148‑1 (MQ=255)                                                                                            tttGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGc                                              <  4:274522/149‑1 (MQ=255)                                                                                               gTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCtt                                            >  5:238519/1‑148 (MQ=255)                                                                                               gTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCtt                                            >  7:30495/1‑148 (MQ=255)                                                                                                tCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTg                                           <  5:129102/148‑1 (MQ=255)                                                                                                  aGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGaa                                         <  4:265406/148‑1 (MQ=255)                                                                                                  aGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAAATTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAg                                        >  5:229804/1‑149 (MQ=255)                                                                                                    gCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAACCCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGc                                      >  6:1996/1‑149 (MQ=255)                                                                                                    gCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGc                                      >  4:166082/1‑149 (MQ=255)                                                                                                      agggaggagGCGCTGGTGCGAAGGGCGGTTTGCGCGAGCGCAAAGCGGGGCTTCGGCGGGCAGCTCGCGGCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCa                                    <  1:237551/146‑1 (MQ=255)                                                                                                       gCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGaa                                         >  4:281748/1‑143 (MQ=255)                                                                                                       gCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGaa                                         <  3:281748/143‑1 (MQ=255)                                                                                                       gCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAg                                   <  5:266462/149‑1 (MQ=255)                                                                                                          aggaggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCa                                >  1:11954/1‑149 (MQ=255)                                                                                                          aggaggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCa                                >  7:120991/1‑149 (MQ=255)                                                                                                            gaggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCAtt                              >  3:258158/1‑149 (MQ=255)                                                                                                             aggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCAttt                             <  2:181422/149‑1 (MQ=255)                                                                                                             aggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCAttt                             <  7:112467/149‑1 (MQ=255)                                                                                                             aggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCAttt                             <  4:258158/149‑1 (MQ=255)                                                                                                              ggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTg                            <  1:258001/149‑1 (MQ=255)                                                                                                               gcgcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTg                            <  5:151530/148‑1 (MQ=255)                                                                                                                 gcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCa                         >  1:331226/1‑149 (MQ=255)                                                                                                                     gTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATcc                      <  4:54118/148‑1 (MQ=255)                                                                                                                            cggcggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAg              >  8:154768/1‑149 (MQ=255)                                                                                                                             gggggTTTGCGGGAGCGCAAAGCGGGGGTTCGGCGGGCAGCTCGCGGCCCAGCCTGGCCTGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGGTGCATTGGTTTAAGGCCAGGCATTTGCCTTCCGCGTcttgt             <  2:308580/146‑5 (MQ=255)                                                                                                                             ccgggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGt                  <  3:65683/141‑1 (MQ=255)                                                                                                                             ccgggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGt                  >  4:65683/4‑144 (MQ=255)                                                                                                                             ccgggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAgg             >  4:30367/4‑149 (MQ=255)                                                                                                                              cgggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGt                  <  2:170984/141‑1 (MQ=255)                                                                                                                              cgggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGt                  >  1:170984/3‑143 (MQ=255)                                                                                                                                ggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg          <  6:127913/149‑1 (MQ=255)                                                                                                                                 gTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCGTTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCggg         >  5:317212/1‑149 (MQ=255)                                                                                                                                 gTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAg                                        <  6:219728/118‑1 (MQ=255)                                                                                                                                 gTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAg                                        >  5:219728/1‑118 (MQ=255)                                                                                                                                     gcgcGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg          >  6:124170/1‑144 (MQ=255)                                                                                                                                     gcgcGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg          <  5:124170/144‑1 (MQ=255)                                                                                                                                       gcgAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCa                                    >  6:296405/1‑116 (MQ=255)                                                                                                                                       gcgAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCa                                    <  5:296405/116‑1 (MQ=255)                                                                                                                                        cgAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCTCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCGAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGc  <  6:317623/149‑1 (MQ=255)                                                                                                                                         gAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg          >  1:335506/1‑140 (MQ=255)                                                                                                                                         gAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg          <  2:335506/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                   GGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGC  >  NC_002947/870094‑870377

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