Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,157,486	G→T	5.8%	G309V (GGA→GTA)	PP_5400 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,157,486	0	G	T	5.8%	254.6 / 3.7	79	G309V (GGA→GTA)	PP_5400	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base T (2/3);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

GCAAACGCAACCTACGGAAATGGCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTC  >  NC_002947/6157341‑6157634                                                                                                                                                     |                                                                                                                                                     gCAAACGCAACCTACGGAAATGGCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAAt                                                                                                                                                   >  3:62560/1‑149 (MQ=255) gCAAACGCAACCTACGGAAATGGCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAAt                                                                                                                                                   <  8:38788/149‑1 (MQ=255)     aCGCAACCTACGGAAATGGCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATggcg                                                                                                                                               >  7:18462/1‑149 (MQ=255)            cTACGGAAATGGCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCAGGTCATGTGGTACTGAATTTTGCCAAGAACAATGATGCCTGGGCGCTAGGAATGGCGGCGATc                                                                                                                                         >  1:267279/1‑148 (MQ=255)            cTAAGGAAATTGGACAGGCAACATTTCCGCATTTGGGCTTGCAGATGCGGCCATTCAATTTTGCAAGGCAAAACGTGCCCATGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGa                                                                                                                                           <  1:126324/149‑1 (MQ=255)             tACGGAAATGGCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCa                                                                                                                                        >  5:103767/1‑148 (MQ=255)               cGGAAATGGCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCaaa                                                                                                                                      <  3:206529/148‑1 (MQ=255)               cGGAAATGGCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAAc                                                                                                                                     <  8:18462/149‑1 (MQ=255)                ggAAATGGCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACt                                                                                                                                    <  5:227561/149‑1 (MQ=255)                ggAAATGGCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACt                                                                                                                                    >  6:227561/1‑149 (MQ=255)                 gAAATGGCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACtt                                                                                                                                   <  2:302437/149‑1 (MQ=255)                       gCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATggc                                                                                                                                                >  1:390565/1‑130 (MQ=255)                       gCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATggc                                                                                                                                                <  2:390565/130‑1 (MQ=255)                         acaGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCATGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGc                                                                                                                           >  7:62965/1‑149 (MQ=255)                          caGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCg                                                                                                                          >  3:202652/1‑149 (MQ=255)                           aGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCgg                                                                                                                         <  6:103767/149‑1 (MQ=255)                           aGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCgg                                                                                                                         <  1:16023/149‑1 (MQ=255)                                aCATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGa                                                                                                                     >  7:157575/1‑148 (MQ=255)                                     tCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATc                                                                                                                                         <  5:266410/123‑1 (MQ=255)                                     tCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATc                                                                                                                                         >  6:266410/1‑123 (MQ=255)                                     tCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAAc                                                                                                                >  5:173523/1‑148 (MQ=255)                                     tCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATAGCGGCGATCAAACTTCTACAGGCGGGTGAAAAAc                                                                                                                >  4:63767/1‑148 (MQ=255)                                       cGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAg                                                                                                             >  7:44954/1‑149 (MQ=255)                                       cGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAg                                                                                                             >  6:300512/1‑149 (MQ=255)                                          aTTGTGGCTTGCAGATGCGTCCATGCAATATTGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGTAATGGCGGCGATCATACTTCTACTGGCGGGTGAAAAACGAGCCg                                                                                                          <  6:127193/149‑1 (MQ=255)                                             ggggCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGc                                                                                                        <  4:62560/148‑1 (MQ=255)                                                    gCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGt                                                                                                >  2:102670/1‑149 (MQ=255)                                                    gCAGAAGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATCTGGTACTGAATATTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACATCTACCGGCGCGTCAAAAACCAGCCGGCATccaaga                                                                                                >  2:300365/1‑144 (MQ=255)                                                      aGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTc                                                                                               >  3:157923/1‑148 (MQ=255)                                                            ggCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCACCGATc                                                                                        >  6:95160/1‑149 (MQ=255)                                                            ggCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGAt                                                                                         >  2:280043/1‑148 (MQ=255)                                                            ggCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGAt                                                                                         <  1:280043/148‑1 (MQ=255)                                                               cATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGTAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGc                                                                                     >  3:333046/1‑149 (MQ=255)                                                               cATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGTAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGc                                                                                     >  3:338765/1‑149 (MQ=255)                                                                 tGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGctc                                                                                   <  2:271801/149‑1 (MQ=255)                                                                   cAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTg                                                                                 <  8:44954/149‑1 (MQ=255)                                                                    aaTATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGAACTGAATTTTGCAAAGCACAATCATGCCTGGCCGCAACGAATGGCGGCGATCAAACTTCTACAGGCGCCTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTgg                                                                                >  3:108767/1‑149 (MQ=255)                                                                        tGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGctct                                                                                  >  3:236000/1‑143 (MQ=255)                                                                        tGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGctct                                                                                  <  4:236000/143‑1 (MQ=255)                                                                        tGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGtt                                                                            >  1:54286/1‑149 (MQ=255)                                                                         ggcaaggcaaACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTg                                                                           >  1:157813/1‑149 (MQ=255)                                                                          gcaaggcaaACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTg                                                                           <  8:69776/148‑1 (MQ=255)                                                                          gcaaggcaaACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTAGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGAGAAAAACCAGCCGGCATCGTAGTCAACGATATGCTCTGGCGTTg                                                                           >  5:312075/1‑148 (MQ=255)                                                                           gaaggcaaACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGGACTATGATGCCTGGGCGCTAGGATTGTCGGGGATCATACTTCTACAGGCGGGTTAAAAACCAGCCGGCATCGTAGTCAACGATCTGGTCTGGCGTTGTa                                                                         <  6:312075/148‑1 (MQ=255)                                                                            aaggcaaACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAg                                                                                                             >  5:273465/1‑112 (MQ=255)                                                                            aaggcaaACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAg                                                                                                             <  6:273465/112‑1 (MQ=255)                                                                              ggcaaTCCGTG‑‑‑TGGATGCCGGTAATGTGGTACTGAATTTTGGCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGTGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTa                                                                      <  7:155107/149‑1 (MQ=255)                                                                               gcaaACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGt                                                                                                   <  1:96220/119‑1 (MQ=255)                                                                               gcaaACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGt                                                                                                   >  2:96220/1‑119 (MQ=255)                                                                                  aaCCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGcc                                                                  <  8:157575/149‑1 (MQ=255)                                                                                     cGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGa                                                               <  5:300512/149‑1 (MQ=255)                                                                                     cGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGa                                                               <  2:54286/149‑1 (MQ=255)                                                                                      gtg‑‑‑tGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGTAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGAt                                                              <  4:338765/149‑1 (MQ=255)                                                                                      gtg‑‑‑tGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGTAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGAt                                                              <  4:333046/149‑1 (MQ=255)                                                                                             ggATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGtt                                                           <  6:173523/148‑1 (MQ=255)                                                                                              gATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTtc                                                          <  8:62965/148‑1 (MQ=255)                                                                                                  ccGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGTCCGGATGTTCTCCCt                                                     >  3:20785/1‑149 (MQ=255)                                                                                                  ccGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCt                                                     <  4:202652/149‑1 (MQ=255)                                                                                                   cGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCt                                                     >  6:162983/1‑148 (MQ=255)                                                                                                              tACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCAt                                          >  7:53609/1‑148 (MQ=255)                                                                                                               aCTGAATTTTGCCAAGCACAATGATGCCTGGCCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGTATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATc                                         >  3:328913/1‑148 (MQ=255)                                                                                                                cTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGa                                       >  4:159045/1‑149 (MQ=255)                                                                                                                 tGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGAt                                      >  4:394489/1‑149 (MQ=255)                                                                                                                   aaTTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCTGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGAtc                                     <  7:67500/148‑1 (MQ=255)                                                                                                                   aaTTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGAtc                                     >  8:67500/1‑148 (MQ=255)                                                                                                                    aTTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGAtctc                                   <  3:394489/149‑1 (MQ=255)                                                                                                                    aTTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGAtctc                                   <  8:53609/149‑1 (MQ=255)                                                                                                                    aTTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGAtctc                                   <  5:162983/149‑1 (MQ=255)                                                                                                                    aTTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATAGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGAtctc                                   <  3:63767/149‑1 (MQ=255)                                                                                                                         gacAAGCACAATGATGCCTGGGCGCTAGGAATGGGGGCGATCAAACTTCTAGAGGGGGGGGAAATCCGAGCCGTCATCTTAGTCAACGATCTGCTGTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACTCCATCGATctcctgt                               <  4:108767/146‑3 (MQ=255)                                                                                                                                acaATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGtgtg                       <  5:84314/149‑1 (MQ=255)                                                                                                                                     gatgCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTc                   >  3:11017/1‑148 (MQ=255)                                                                                                                                        gCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGc               >  5:352008/1‑149 (MQ=255)                                                                                                                                           tGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGctct                                                                                  >  1:207744/1‑79 (MQ=255)                                                                                                                                           tGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGctct                                                                                  <  2:207744/79‑1 (MQ=255)                                                                                                                                             ggCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGaa           >  1:68423/1‑148 (MQ=255)                                                                                                                                              gcgcTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGAGGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATa         <  5:95160/149‑1 (MQ=255)                                                                                                                                               cgcTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAg        >  8:44712/1‑149 (MQ=255)                                                                                                                                                   aGGACTGGCGGCGATCAAACTTCTACAGCCGGGTGAAACCCCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGAtt    <  2:157813/149‑1 (MQ=255)                                                                                                                                                    ggAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGAttt   >  8:364004/1‑149 (MQ=255)                                                                                                                                                     gAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTc  <  7:44712/149‑1 (MQ=255)                                                                                                                                                     gAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAAGAGGAtt    >  8:265587/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCAAACGCAACCTACGGAAATGGCACAGGCAACATGTCCGCATTGGGGCTTGCAGATGCGGCCATGCAATATGGCAAGGCAAACCGTG‑‑‑TGGATCCCGGTCATGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTC  >  NC_002947/6157341‑6157634

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