Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,823,390	C→A	6.0%	L2615F (TTG→TTT)	pvdL ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,823,390	0	C	A	6.0%	220.4 / 3.0	67	L2615F (TTG→TTT)	pvdL	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base A (2/2);  total (36/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.70e-01

AGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTACCC  >  NC_002947/4823244‑4823534                                                                                                                                                   |                                                                                                                                                 agcagcACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAg                                                                                                                                               >  1:1285828/1‑150 (MQ=255)  gcagcaCCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGt                                                                                                                                              <  1:712617/150‑1 (MQ=255)  gcagcaCCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGt                                                                                                                                              <  2:1285828/150‑1 (MQ=255)     gcaCCAGATGCGCAGGCCCCATGTATTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCa                                                                                                                                           >  2:837659/1‑150 (MQ=255)       aCCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAgc                                                                                                                                         <  1:531007/150‑1 (MQ=255)       aCCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAgc                                                                                                                                         <  1:134184/150‑1 (MQ=255)         cAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAgcgc                                                                                                                                       >  2:113425/1‑150 (MQ=255)         cAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAgcgc                                                                                                                                       >  1:934400/1‑150 (MQ=255)                  aGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCacgac                                                                                                                              >  2:396164/1‑150 (MQ=255)                     cACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGc                                                                                                                           <  1:774869/150‑1 (MQ=255)                           gTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCtcgtc                                                                                                                     <  2:243832/150‑1 (MQ=255)                                    cAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTc                                                                                                            <  2:883995/150‑1 (MQ=255)                                       gTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTg                                                                                                         >  2:933721/1‑150 (MQ=255)                                               gggCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCa                                                                                                 <  1:113425/150‑1 (MQ=255)                                                  cATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAggc                                                                                              <  1:396164/150‑1 (MQ=255)                                                   aTCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCacgacg                                                                                                                             <  1:341819/118‑1 (MQ=255)                                                   aTCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCacgacg                                                                                                                             >  2:341819/1‑118 (MQ=255)                                                   aTCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAggcg                                                                                             >  2:1265347/1‑150 (MQ=255)                                                     cGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCTCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAggcggc                                                                                           >  1:487414/1‑150 (MQ=255)                                                     cGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAggcggc                                                                                           >  1:241658/1‑150 (MQ=255)                                                      ggCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGct                                                                                          >  2:843091/1‑150 (MQ=255)                                                       gCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGAAAGTCCAGCGCCAGCACGACGGcc                                                                                                                          >  2:1086902/1‑117 (MQ=255)                                                       gCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGAAAGTCCAGCGCCAGCACGACGGcc                                                                                                                          <  1:1086917/117‑1 (MQ=255)                                                       gCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGAAAGTCCAGCGCCAGCACGACGGcc                                                                                                                          <  1:1086902/117‑1 (MQ=255)                                                       gCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGAAAGTCCAGCGCCAGCACGACGGcc                                                                                                                          >  2:1086917/1‑117 (MQ=255)                                                         tttCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAgca                                                                                                                                             >  2:1000464/1‑94 (MQ=255)                                                         tttCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAgca                                                                                                                                             <  1:1000464/96‑3 (MQ=255)                                                         tttCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCtctc                                                                                       >  2:655281/1‑150 (MQ=255)                                                             aCCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGAtt                                                                                   >  1:1063270/1‑150 (MQ=255)                                                             aCCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGAtt                                                                                   >  1:577659/1‑150 (MQ=255)                                                                 cTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGc                                                                               <  1:655281/150‑1 (MQ=255)                                                                  ttCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGctc                                                                                         <  2:400127/139‑1 (MQ=255)                                                                  ttCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGctc                                                                                         >  1:400127/1‑139 (MQ=255)                                                                  ttCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCCGCACGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGcc                                                                              >  1:1095298/1‑150 (MQ=255)                                                                       gcagcaGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGtt                                                                         >  2:1297928/1‑150 (MQ=255)                                                                             gttgttGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGc                                                                   >  2:528230/1‑150 (MQ=255)                                                                                 ttgCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGaa                                                               <  1:1417313/150‑1 (MQ=255)                                                                                      cgtcgtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGctc                                                                                         <  2:195978/119‑1 (MQ=255)                                                                                      cgtcgtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGctc                                                                                         >  1:195978/1‑119 (MQ=255)                                                                                       gtcgtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAg                                                                                                                                          >  2:954999/1‑69 (MQ=255)                                                                                       gtcgtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAg                                                                                                                                          <  1:954999/69‑1 (MQ=255)                                                                                          gtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGAt                                                      >  2:477874/1‑150 (MQ=255)                                                                                                ggcggTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCAc                                                <  2:934400/150‑1 (MQ=255)                                                                                                    gTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACctgc                                            <  1:968052/150‑1 (MQ=255)                                                                                                    gTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACctgc                                            <  2:1291209/150‑1 (MQ=255)                                                                                                     tGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACctgct                                           >  2:362350/1‑150 (MQ=255)                                                                                                         ctgcctgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGc                                                                   <  2:1321816/122‑1 (MQ=255)                                                                                                         ctgcctgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGc                                                                   >  1:1321816/1‑122 (MQ=255)                                                                                                             ctgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCa                                   <  1:477874/150‑1 (MQ=255)                                                                                                              tgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCAc                                  >  2:483949/1‑150 (MQ=255)                                                                                                               gcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACg                                 <  1:110537/150‑1 (MQ=255)                                                                                                                   acCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCa                                                                            <  2:1243371/103‑1 (MQ=255)                                                                                                                   acCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCa                                                                            >  1:1243371/1‑103 (MQ=255)                                                                                                                           ggccaccaGCTGCTTGCCGCCCGGCAAGTCCAGCGGCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCGAGTTCGATGCGGAAGCCGCGGACCTTCACCTGCTGGTCGGCACGGCCGACATACTc                     <  2:1132706/148‑1 (MQ=255)                                                                                                                            cccACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCAc                                                >  1:625934/1‑122 (MQ=255)                                                                                                                            cccACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCAc                                                <  2:625934/122‑1 (MQ=255)                                                                                                                             ccaccaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGa                                                                                     <  2:199392/84‑1 (MQ=255)                                                                                                                             ccaccaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGa                                                                                     >  1:199392/1‑84 (MQ=255)                                                                                                                             ccaccaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTcca                   >  1:489188/1‑150 (MQ=255)                                                                                                                             ccaccaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTcca                   >  1:1099419/1‑150 (MQ=255)                                                                                                                              caccaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAATCGGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTccac                  >  1:661638/1‑150 (MQ=255)                                                                                                                              caccaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTccac                  <  1:528230/150‑1 (MQ=255)                                                                                                                                  ggctgctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAgg              >  2:355758/2‑150 (MQ=255)                                                                                                                                  ctggactTGCCGGGCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGACGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAgg              >  2:1225360/6‑150 (MQ=255)                                                                                                                                      gctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCAc                                                >  1:16174/1‑112 (MQ=255)                                                                                                                                      gctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCAc                                                <  2:16174/112‑1 (MQ=255)                                                                                                                                              cccGGCAAGTCCAGCGCCAGCAGGACGGCCTCGTCGACATCCGCGTGCGCTAGCAGGCGGCTCTCGACTGCGCCCATTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTAccc  <  1:601854/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTACCC  >  NC_002947/4823244‑4823534

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑