Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,023,499:1	+C	5.3%	coding (531/897 nt)	rmd →	oxidoreductase Rmd

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,023,499	1	.	C	5.3%	306.1 / 10.6	75	coding (531/897 nt)	rmd	oxidoreductase Rmd
Reads supporting (aligned to +/- strand):  ref base . (35/36);  new base C (2/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

GACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGATCGAGGCACGCCCGCTGGGCCAGATA  >  NC_002947/2023358‑2023643                                                                                                                                               |                                                                                                                                                  gACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGccc                                                                                                                                             >  1:42623/1‑148 (MQ=255)  aCCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCA‑GCAGGAAATTCTCTGCCTGCCCaa                                                                                                                                           <  5:187428/149‑1 (MQ=255)      cGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGAtt                                                                                                                                       >  1:217339/1‑149 (MQ=255)            cgcCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGGCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTc                                                                                                                                  <  8:80640/148‑1 (MQ=255)               cGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCAt                                                                                                                              >  1:367513/1‑149 (MQ=255)                 gCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCAttt                                                                                                                            >  4:189993/1‑149 (MQ=255)                   cAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTc                                                                                                                          >  7:272432/1‑149 (MQ=255)                       gATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTg                                                                                                                                      <  1:17893/133‑1 (MQ=255)                       gATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTg                                                                                                                                      >  2:17893/1‑133 (MQ=255)                        aTTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTcgcc‑gc                                                                                                                    <  2:42623/149‑1 (MQ=255)                          tACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTcgcc‑gc                                                                                                                    >  3:159405/1‑147 (MQ=255)                           aCGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCaa                                                                                                                  <  2:367513/148‑1 (MQ=255)                             gcgGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGc                                                                                                               >  7:286890/1‑149 (MQ=255)                               ggTCAGCAATTTGTCCTTGTAATACCTGGCGCGGCTTTTGTGGTACCGTTTTCCTTTCGTCTTCACGCTTCCTTTCTATTACACCGGTGTTGGGCAGGCAGAGATTTTCCTG‑CGGCGCAAGATTGTTTCGCATTTTCTCC‑TAAAGGCCg                                                                                                             <  6:167145/149‑1 (MQ=255)                                 tCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGa                                                                                                            >  7:83423/1‑148 (MQ=255)                                    gCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACAcc                                                                                                        <  4:350571/149‑1 (MQ=255)                                     cAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGcc                                                                                                              <  4:160070/142‑1 (MQ=255)                                     cAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGcc                                                                                                              >  3:160070/1‑142 (MQ=255)                                     cAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGaca                                                                                                          <  2:251603/146‑1 (MQ=255)                                     cAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGaca                                                                                                          >  1:251603/1‑146 (MQ=255)                                        gTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCg                                                                                                    <  8:281017/149‑1 (MQ=255)                                         ttGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGa                                                                                                   <  4:159405/149‑1 (MQ=255)                                         cgcgCCATGTATTCCATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATc                                                                                                     >  3:349806/4‑147 (MQ=255)                                           ggCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAg                                                                                                  <  4:75132/148‑1 (MQ=255)                                             ccATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTc                                                                                               >  3:106821/1‑149 (MQ=255)                                              cATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATc                                                                                                     >  7:339850/1‑142 (MQ=255)                                                   attccATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCg                                                                                              >  3:245852/5‑144 (MQ=255)                                                   attccATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCg                                                                                              <  4:245852/140‑1 (MQ=255)                                                   aaTACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAAc                                                                                         <  8:50105/149‑1 (MQ=255)                                                       cATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTa                                                                                      <  1:190398/148‑1 (MQ=255)                                                        aTGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCCCGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GAAAGGCCGACACCATCGAGCTCGGCAACCTAGa                                                                                    >  6:287532/1‑149 (MQ=255)                                                        aTGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGa                                                                                    >  5:169480/1‑149 (MQ=255)                                                        aTGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑CCAAGGCCGACACCATCGAGCTCGGCAACCTAGa                                                                                    >  7:352459/1‑149 (MQ=255)                                                         tGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGAc                                                                                   <  8:272432/149‑1 (MQ=255)                                                         tGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGAc                                                                                   >  2:305931/1‑149 (MQ=255)                                                          ggCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGAc                                                                                   <  4:349806/148‑1 (MQ=255)                                                           gcgcGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACg                                                                                  <  6:169480/148‑1 (MQ=255)                                                            cgcgGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTa                                                                                <  2:217339/149‑1 (MQ=255)                                                            cgcgGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTATCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTa                                                                                <  5:121661/149‑1 (MQ=255)                                                             gcgGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTAt                                                                               <  5:241384/149‑1 (MQ=255)                                                                      gtgtGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATc                                                                                                     >  6:20776/1‑118 (MQ=255)                                                                      gtgtGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATc                                                                                                     <  5:20776/118‑1 (MQ=255)                                                                        gtgACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGAttt                                                                     >  5:363627/1‑148 (MQ=255)                                                                              gTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGAt                                                              <  7:188992/149‑1 (MQ=255)                                                                               tttGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGAAACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATg                                                             >  7:55860/1‑149 (MQ=255)                                                                                       tCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTg                                                       <  8:361226/147‑1 (MQ=255)                                                                                       tCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGc                                                      <  8:83423/148‑1 (MQ=255)                                                                                       tCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGc                                                      <  6:363627/148‑1 (MQ=255)                                                                                       tCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCg                                                     >  6:44416/1‑149 (MQ=255)                                                                                            atcaCGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCTACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGt                                                <  7:259844/149‑1 (MQ=255)                                                                                            atcaCGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCg                                                 <  1:305931/148‑1 (MQ=255)                                                                                                cgcgTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAg                                            >  4:76978/1‑149 (MQ=255)                                                                                                 gcgTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAgg                                           >  4:231662/1‑149 (MQ=255)                                                                                                   gTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCg                                         >  8:147535/1‑149 (MQ=255)                                                                                                   gTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGAGATTTCAGCGATGTACGTGCGGTCGTGCAGGCg                                         >  4:384613/1‑149 (MQ=255)                                                                                                            aTTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGc                                <  3:189993/149‑1 (MQ=255)                                                                                                                 acCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGAt                           >  3:116659/1‑149 (MQ=255)                                                                                                                   cGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTAc                                                          <  3:350816/116‑1 (MQ=255)                                                                                                                   cGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTAGGGCGCGATTTCAGCGATGTAc                                                          >  4:350816/1‑116 (MQ=255)                                                                                                                    ggTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGATc                          >  7:97375/1‑147 (MQ=255)                                                                                                                     gtgtTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGaacccg                       >  8:358692/1‑144 (MQ=255)                                                                                                                     gtgtTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGATCGAg                       <  3:231662/149‑1 (MQ=255)                                                                                                                        ttGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGGCGTGCAGGCGTACCGCGGCTTGATCGAGGCa                    >  3:212220/1‑149 (MQ=255)                                                                                                                        ttGGGCA‑GCAGGAAATTCTCTGCCTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACcg                                    >  3:128469/1‑133 (MQ=255)                                                                                                                        ttGGGCA‑GCAGGAAATTCTCTGCCTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACcg                                    <  4:128469/133‑1 (MQ=255)                                                                                                                        ttGGGCA‑GCAGGAAATTCTCTGCCTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGATCGAGGCa                    >  3:304659/1‑149 (MQ=255)                                                                                                                          gggCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGATCGAGGCACg                  >  8:205274/1‑149 (MQ=255)                                                                                                                               ggcagAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGATCGAGGCACGCCCGc             <  8:97375/149‑1 (MQ=255)                                                                                                                                 cagAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGc                                              >  1:176883/1‑114 (MQ=255)                                                                                                                                 cagAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGc                                              <  2:176883/114‑1 (MQ=255)                                                                                                                                 cagAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGATCGAGGCACGCCCGCTg           <  4:106821/149‑1 (MQ=255)                                                                                                                                  agagAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGATCGAGGCACGCCCGCTgg          <  3:76978/149‑1 (MQ=255)                                                                                                                                    agaATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGATCGAGGCACGCCCGCTGGGc        <  5:44416/149‑1 (MQ=255)                                                                                                                                      aaTTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCCGGCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGATCGAGGCACGCCCGCTGGGc        >  1:277162/1‑148 (MQ=255)                                                                                                                                           ccTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGATCGAGGCACGCCCGCTGGGCCAGATa  <  8:286890/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                  GACCCCGCCAGCGCCGGCCAATGATTACGCGGTCAGCAAGTTGGCCATGGAATACATGGCGCGGCTGTGGTGTGACCGTTTGCCTATCGTCATCACGCGTCCTTTCAATTACACCGGTGTTGGGCAGGCAGAGAATTTCCTG‑CTGCCCAAGATTGTTTCGCATTTTCGCC‑GCAAGGCCGACACCATCGAGCTCGGCAACCTAGACGTATGGCGCGATTTCAGCGATGTACGTGCGGTCGTGCAGGCGTACCGCGGCTTGATCGAGGCACGCCCGCTGGGCCAGATA  >  NC_002947/2023358‑2023643

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑