Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	938,369	C→G	5.2%	G2418A (GGG→GCG)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	938,369	0	C	G	5.2%	281.7 / 4.1	89	G2418A (GGG→GCG)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (47/36);  new base G (2/3);  total (49/40)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

GCATTGCCGCTGGCATCGGTCAGCACCACTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGTTGCCGGCGGCGTCGGTGGCGGT  >  NC_002947/938222‑938517                                                                                                                                                    |                                                                                                                                                     gCATTGCCGCTGGCATCGGTCAGCACCACTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGcccc                                                                                                                                                     <  7:182143/149‑1 (MQ=255)   aTTGCCGCTGGCATCGGTCAGCACCACTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCgg                                                                                                                                                   >  2:226878/1‑149 (MQ=255)        cGCTGGCATCGGTCAGCACCACTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCGCGGcgaag                                                                                                                                              >  3:192737/1‑145 (MQ=255)        cGCTGGCATCGGTCAGCACCACTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTc                                                                                                                                               <  7:188702/148‑1 (MQ=255)          cTGGCATCGGTCAGCACCACTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTcgcc                                                                                                                                            >  4:226355/1‑149 (MQ=255)                   gTCAGCACCACTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCAtt                                                                                                                                   >  5:122262/1‑149 (MQ=255)                    tCAGCACCACTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCGCGGCGAAGCCGGGGCCagcg                                                                                                                                  <  4:192737/149‑4 (MQ=255)                       gCACCACTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGc                                                                                                                               >  4:92220/1‑149 (MQ=255)                         accacTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTg                                                                                                                             >  1:107706/1‑149 (MQ=255)                          ccacTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTgg                                                                                                                            <  3:125429/149‑1 (MQ=255)                          ccacTTGCAGGTTGCTGCCATCTTTCACCGGGGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTgg                                                                                                                            <  1:167123/149‑1 (MQ=255)                           cacTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTggg                                                                                                                           >  7:138101/1‑149 (MQ=255)                           cacTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTAACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTggg                                                                                                                           >  7:176538/1‑149 (MQ=255)                                 cAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCg                                                                                                                      <  2:92569/148‑1 (MQ=255)                                  aGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCggc                                                                                                                    >  5:132654/1‑149 (MQ=255)                                     ttGCTGTCAGCTTTCCCCCCTCTGCCAACGGTGATCTCTAAGGTGCCGGCCTCCTGGACCGCGCCGTATGAGAGTACTTTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCggcc                                                                                                                   <  8:238750/147‑1 (MQ=255)                                     ttGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCggccag                                                                                                                 <  5:23836/149‑1 (MQ=255)                                         tgcCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCggccaggcc                                                                                                              >  4:292027/1‑148 (MQ=255)                                         tgcCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCggccaggcc                                                                                                              >  2:176189/1‑148 (MQ=255)                                         tgcCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCggccaggcc                                                                                                              >  7:129186/1‑148 (MQ=255)                                         tgcCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTAGGCGCCggccaggcc                                                                                                              >  1:144007/1‑148 (MQ=255)                                         tgcCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACATGCACGCTGTCCCCGGCTTCGCCGCGGCCATTGAGCTTGGCGCCggtcaggcc                                                                                                              >  5:211331/1‑148 (MQ=255)                                         tgcCATCTTTCACAGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCAAGTACATTTCCCTGTGCATCGCGCACCTTCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGTCGCCggccaggcc                                                                                                              >  1:209119/1‑148 (MQ=255)                                          gcCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCaa                                                                                                            >  4:283571/1‑149 (MQ=255)                                          gcCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCaa                                                                                                            >  1:281890/1‑149 (MQ=255)                                           ccATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAAc                                                                                                           >  6:353762/1‑149 (MQ=255)                                            cATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAAc                                                                                                           >  8:6444/1‑148 (MQ=255)                                              tCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCGCGGCGAAGCCGGGGCCagcgtgc                                                                                                                               <  6:153303/126‑7 (MQ=255)                                                tttCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTAcc                                                                                                       >  2:158329/1‑148 (MQ=255)                                                tttCACCGGTGGGTCAAGGGTGATCTGGAAGGCGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTAcc                                                                                                       >  5:382367/1‑148 (MQ=255)                                                 ttCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTAcc                                                                                                       >  5:300840/1‑147 (MQ=255)                                                  tCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCaa                                                                                                     <  6:122262/148‑1 (MQ=255)                                                   cACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGGCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAgt                                                                                                   >  1:311001/1‑149 (MQ=255)                                                   cACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAg                                                                                                    <  2:107706/148‑1 (MQ=255)                                                   cACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGGTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAg                                                                                                    <  7:306885/148‑1 (MQ=255)                                                    agcGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTAACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAgt                                                                                                   <  8:176538/146‑1 (MQ=255)                                                    aCCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAgtt                                                                                                  >  1:193688/1‑149 (MQ=255)                                                      cGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATTGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAgttgg                                                                                                >  6:56561/1‑149 (MQ=255)                                                      cGGTGGGTCAAGGGGGTTCTGGATGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGGTGGCCCCGGGTTCGCCGCGGCGATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAgtt                                                                                                  <  3:125627/147‑1 (MQ=255)                                                       ggTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTAGATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCTGCTTCGCCGCGGCCATTGTGCTGGGCGCCGGCGAGGCCAACTACCAAgttggt                                                                                               <  2:251509/149‑1 (MQ=255)                                                         tGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTTGGCGCCGGCCAGGCCAACTACCAAgttggttg                                                                                             <  2:144007/149‑1 (MQ=255)                                                           ggTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCg                                                                                           >  2:144684/1‑149 (MQ=255)                                                               aaGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCAt                                                                                       >  5:181661/1‑149 (MQ=255)                                                                 gggTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATcc                                                                                     >  8:311315/1‑149 (MQ=255)                                                                   gTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCgg                                                                                   >  7:200462/1‑149 (MQ=255)                                                                      aTCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGACAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCTCCTCTGCCATTGAGCTGGGCGCAGGACAGGCCAACTCCCAAGTTGGTTGAGGCATCCggtgg                                                                                >  3:333805/1‑149 (MQ=255)                                                                        cTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGgtggtg                                                                              >  8:6052/1‑149 (MQ=255)                                                                          ggAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGgtggtgtg                                                                            >  1:364686/1‑149 (MQ=255)                                                                           gAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACGAAGTTGGTTGCGGCATCCGgtggtgtg                                                                            <  1:176189/148‑1 (MQ=255)                                                                           gAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCggt                                                                                  <  6:161166/142‑1 (MQ=255)                                                                           gAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCggt                                                                                  >  5:161166/1‑142 (MQ=255)                                                                           gAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGgtggtgtg                                                                            <  3:292027/148‑1 (MQ=255)                                                                            aaGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGgtggtgtggt                                                                          <  5:353762/149‑1 (MQ=255)                                                                               gTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTc                                                                       <  3:182293/149‑1 (MQ=255)                                                                               gTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGACAGGCCAACTACCAAGTTGGTGTCTGCATCCGGTGGTGTGGTGTc                                                                       >  8:260546/1‑149 (MQ=255)                                                                                  ccGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGcc                                                                    >  3:326957/1‑149 (MQ=255)                                                                                    gTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCAc                                                                  >  5:99380/1‑149 (MQ=255)                                                                                    gTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCAc                                                                  <  3:107740/149‑1 (MQ=255)                                                                                           ttGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGAt                                                           <  2:193688/149‑1 (MQ=255)                                                                                             gACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCAcc                                                                 <  7:195482/141‑1 (MQ=255)                                                                                             gACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCAcc                                                                 >  8:195482/1‑141 (MQ=255)                                                                                             gACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGt                                                         >  1:163864/1‑149 (MQ=255)                                                                                               ccGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCACGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTcc                                                       <  2:209119/149‑1 (MQ=255)                                                                                                   gCCGTTGGCCTGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCTGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGgcg                                                   <  3:252835/149‑1 (MQ=255)                                                                                                     cGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGcgcg                                                 <  7:6052/149‑1 (MQ=255)                                                                                                      gTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCgg                                                <  3:92220/149‑1 (MQ=255)                                                                                                       ttGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGt                                               <  1:226878/149‑1 (MQ=255)                                                                                                            cAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGt                                          >  1:174588/1‑149 (MQ=255)                                                                                                            cAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGt                                          >  4:55039/1‑149 (MQ=255)                                                                                                             aGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTc                                         <  1:144684/149‑1 (MQ=255)                                                                                                             aGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTc                                         >  1:328484/1‑149 (MQ=255)                                                                                                                     gCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCgccggccg                                 >  8:3965/1‑149 (MQ=255)                                                                                                                      cccTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCgccggccg                                 <  6:300840/148‑1 (MQ=255)                                                                                                                       ccTGGGCATCGCGCACCTGCACGCTGTCCGCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGGCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTg                               <  4:333805/149‑1 (MQ=255)                                                                                                                       ccTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTg                               <  2:163864/149‑1 (MQ=255)                                                                                                                        cTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTgg                              <  2:311001/149‑1 (MQ=255)                                                                                                                        cTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTgg                              <  3:55039/149‑1 (MQ=255)                                                                                                                        cTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTgg                              <  4:326957/149‑1 (MQ=255)                                                                                                                         tGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTg                               >  6:114113/1‑147 (MQ=255)                                                                                                                               tCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTTCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGtt                       >  7:242573/1‑149 (MQ=255)                                                                                                                               tCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGtt                       <  8:138101/149‑1 (MQ=255)                                                                                                                               tCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGtt                       <  3:226355/149‑1 (MQ=255)                                                                                                                                     aCCTGCACGCTGGTCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGTTGCcggc                 >  2:86242/1‑149 (MQ=255)                                                                                                                                     aCCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGTTGCcggc                 >  6:186477/1‑149 (MQ=255)                                                                                                                                       cTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTg                               >  7:44096/1‑133 (MQ=255)                                                                                                                                       cTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTg                               <  8:44096/133‑1 (MQ=255)                                                                                                                                           aCGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGTTGCCGGCGGCGTc           <  8:129186/149‑1 (MQ=255)                                                                                                                                            cGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGTTGCCGGCGGCGTCg          <  6:99380/149‑1 (MQ=255)                                                                                                                                             gCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGTTGCCGGCGGCGTCgg         >  8:320737/1‑149 (MQ=255)                                                                                                                                              atGGCCGCGGCGAAGCCGGGGCCATTGAGCAGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGTTGCCGGCGGCGTCGGt        >  6:106190/2‑149 (MQ=255)                                                                                                                                                    ccGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGTTGCCGGCGGCGTCGGTGGCGGt  <  1:158329/149‑1 (MQ=255)                                                                                                                                                    ccGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGTTGCCGGCGGCGTCGGTGGCGGt  <  6:181661/149‑1 (MQ=255)                                                                                                                                                    ccGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGTTGCCGGCGGCGTCGGTGGCGGt  <  2:328484/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GCATTGCCGCTGGCATCGGTCAGCACCACTTGCAGGTTGCTGCCATCTTTCACCGGTGGGTCAAGGGTGATCTGGAAGGTGCCGTCCGGGTTGACCGAGCCGTTGGCCAGTACATTGCCCTGGGCATCGCGCACCTGCACGCTGGCCCCGGCTTCGCCGCGGCCATTGAGCTGGGCGCCGGCCAGGCCAACTACCAAGTTGGTTGCGGCATCCGGTGGTGTGGTGTCGCCACCGTCGATGTCCGGCGCGGTGATGTCGCCGGCCGTGGAGGTGTTGCCGGCGGCGTCGGTGGCGGT  >  NC_002947/938222‑938517

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