Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,539,492	G→T	8.5%	A498A (GCC→GCA)	PP_4872 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,539,492	0	G	T	8.5%	195.1 / 5.8	66	A498A (GCC→GCA)	PP_4872	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/28);  new base T (3/3);  total (35/31)
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.27e-01

CACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCC  >  NC_002947/5539347‑5539635                                                                                                                                                    |                                                                                                                                                cACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGcc                                                                                                                                              >  2:143630/1‑149 (MQ=255)       cGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCgg                                                                                                                                        <  8:187964/149‑1 (MQ=255)        ggTGTGCGCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGCCGTACTCCGGGGATTCCACCGCGAGGTCGTAGCGCAGGCCGGAAGCACTCAGG‑ATCTT‑TTTCACACCCGGCAAGGCCCGTGCCTTGccgg                                                                                                                                       <  8:94280/149‑4 (MQ=255)               tCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGTCAAGGCACGGGCCTTGCGGTAAAGCTc                                                                                                                                >  5:75740/1‑149 (MQ=255)                 cGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGa                                                                                                                              >  1:270090/1‑149 (MQ=255)                      cAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAAc                                                                                                                         <  2:270090/149‑1 (MQ=255)                       aaTCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAAc                                                                                                                         >  6:186632/1‑148 (MQ=255)                           ttCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGTTGATCGAGCt                                                                                                                    >  1:281681/1‑149 (MQ=255)                                gTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGTTGATCGAGCtttacc                                                                                                               <  2:281681/149‑6 (MQ=255)                                gTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGt                                                                                                               >  5:242224/1‑149 (MQ=255)                                  aGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCg                                                                                                             >  4:317447/1‑149 (MQ=255)                                         aaCGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGa                                                                                                      <  5:362952/149‑1 (MQ=255)                                         aaCGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGa                                                                                                      <  2:303898/149‑1 (MQ=255)                                              ggtgggtgACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGttt                                                                                                 >  5:335868/1‑149 (MQ=255)                                               gtgggggACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGtttt                                                                                                <  4:60813/149‑1 (MQ=255)                                                tgggtgACCAGTTCCTTGACGTACTCGGGCGATTCCACAGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTc                                                                                               <  5:221332/149‑1 (MQ=255)                                                  ggtgACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGc                                                                                             <  8:353840/149‑1 (MQ=255)                                                   gtgACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCa                                                                                            <  6:335868/149‑1 (MQ=255)                                                     gACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCg                                                                                              >  7:302957/1‑145 (MQ=255)                                                          gTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTc                                                                                                              <  8:145456/124‑1 (MQ=255)                                                          gTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTCTGGTc                                                                                                              >  7:145456/1‑124 (MQ=255)                                                           ttCCTTGACGTACTCCGGCGATTCCTCGGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGT‑ATCTT‑TTTCTCACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCg                                                                                    <  5:313498/149‑1 (MQ=255)                                                             ccTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGGCATCTT‑TTTCACACCCGGCAAGGCCCGTGCCTTGcc                                                                                                                                         <  1:247202/95‑2 (MQ=255)                                                             ccTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCCCGTGCCTTGcc                                                                                                                                         >  2:247202/1‑93 (MQ=255)                                                             ccTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGCCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCgg                                                                                   <  4:111645/148‑1 (MQ=255)                                                             ccTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCCGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTCGTCGGTATTGAGGTTTTCGCAGATACCCgg                                                                                   >  5:51725/1‑148 (MQ=255)                                                              cTTGACGTACTTCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCggg                                                                                  >  8:27504/1‑148 (MQ=255)                                                                  aCGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATAc                                                                             <  3:317447/149‑1 (MQ=255)                                                                          cGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCacgacg                                                                     >  4:153577/1‑149 (MQ=255)                                                                          cGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCacgacg                                                                     <  6:242224/149‑1 (MQ=255)                                                                          cGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCacgac                                                                      <  1:252832/148‑1 (MQ=255)                                                                            gCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACACGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGCCGGc                                                                   >  6:352308/1‑149 (MQ=255)                                                                               aTTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTg                                                                >  2:56955/1‑149 (MQ=255)                                                                                   cACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCgg                                                             <  6:366427/148‑1 (MQ=255)                                                                                     ccGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAg                                                          >  6:116901/1‑149 (MQ=255)                                                                                     ccGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAg                                                          <  3:203047/149‑1 (MQ=255)                                                                                          aGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGGGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACtt                                                     >  2:332019/1‑149 (MQ=255)                                                                                           ggTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTTGTTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACtt                                                     >  1:242097/1‑149 (MQ=255)                                                                                           ggTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTtc                                                    >  6:164257/1‑149 (MQ=255)                                                                                                  gcgcAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGAttt                                             <  5:164257/149‑1 (MQ=255)                                                                                                        gCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGGAGTACTTCTCGATTTCGGGGc                                       >  8:104289/1‑149 (MQ=255)                                                                                                          cGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCggg                                         >  4:229471/1‑145 (MQ=255)                                                                                                          cGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCggg                                         <  3:229471/145‑1 (MQ=255)                                                                                                            gAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGTTACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCtttt                                   >  2:287752/1‑149 (MQ=255)                                                                                                               gCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCa                                <  7:27504/149‑1 (MQ=255)                                                                                                                  aTCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGc                             >  1:187490/1‑149 (MQ=255)                                                                                                                  aTCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGc                             >  4:281512/1‑149 (MQ=255)                                                                                                                       g‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCg                            <  8:199121/145‑1 (MQ=255)                                                                                                                       g‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGc                        >  1:258880/1‑149 (MQ=255)                                                                                                                       g‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACACGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCg                            >  7:199121/1‑145 (MQ=255)                                                                                                                        cgatct‑tttCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGc                        <  3:262338/144‑1 (MQ=255)                                                                                                                            tt‑tttCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTa                    >  1:90662/1‑149 (MQ=255)                                                                                                                                ttCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGc                        >  4:262338/1‑142 (MQ=255)                                                                                                                                   acacCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACTTGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGt               <  6:51725/148‑1 (MQ=255)                                                                                                                                   acacCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGa                                                <  8:76174/115‑1 (MQ=255)                                                                                                                                   acacCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGa                                                >  7:76174/1‑115 (MQ=255)                                                                                                                                    tttaccgCAAGGCCCGTGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTg             >  1:81462/7‑149 (MQ=255)                                                                                                                                    tttaccgCAAGGCCCGTGCCTTGCCGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGAt                          <  6:149239/130‑1 (MQ=255)                                                                                                                                    tttaccgCAAGGCCCGTGCCTTGCCGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGAt                          >  5:149239/7‑136 (MQ=255)                                                                                                                                    cacCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTg             <  1:143630/149‑1 (MQ=255)                                                                                                                                     acCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTgg            <  2:258880/149‑1 (MQ=255)                                                                                                                                       ccGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCa          <  1:56955/149‑1 (MQ=255)                                                                                                                                         ggCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGt        <  3:153577/149‑1 (MQ=255)                                                                                                                                            aaGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCgg     >  6:147411/1‑149 (MQ=255)                                                                                                                                               gCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGccc  >  4:323215/1‑149 (MQ=255)                                                                                                                                               gCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGccc  >  8:78355/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                CACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGCGATTCCACCGCCAGGTCGTAGCGCAGGCCGGAAGCAATCAGG‑ATCTT‑TTTCACACCCGGCAAGGCACGGGCCTTGCGGTAAAGCTCGATCAACGAGCTGTGGTCGGTATTGAGGTTTTCGCAGATACCCGGGAATACGCACGACGGCTTGCGGCAGTACTTCTCGATTTCGGGGCTTTTGCAGGCGATGCGGTACATGTTGGCAGTCGGCCC  >  NC_002947/5539347‑5539635

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