Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	225,369	A→G	7.7%	S98P (TCG→CCG)	PP_0173 ←	transcriptional factor‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	225,369	0	A	G	7.7%	179.2 / 4.9	65	S98P (TCG→CCG)	PP_0173	transcriptional factor‑like protein
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base G (3/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

TGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCGGCGCTGGTTATGCCT  >  NC_002947/225223‑225517                                                                                                                                                   |                                                                                                                                                     tGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGAcc                                                                                                                                                    <  1:40262/149‑1 (MQ=255) tGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGAc                                                                                                                                                     >  4:348678/1‑148 (MQ=255)   gcgcATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCa                                                                                                                                                   <  8:128876/148‑1 (MQ=255)    cgcATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGt                                                                                                                                                 <  7:214984/149‑1 (MQ=255)       aTCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGt                                                                                                                                                 <  8:183617/146‑1 (MQ=255)       aTCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGt                                                                                                                                                 >  7:183617/1‑146 (MQ=255)          ccGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGTCGGCCAGTGccgtg                                                                                                                                           <  3:301257/149‑4 (MQ=255)          ccGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAAc                                                                                                                                           >  7:3859/1‑149 (MQ=255)           cGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGTAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCaa                                                                                                                                            >  1:159774/1‑147 (MQ=255)           cGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGTAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCaa                                                                                                                                            <  2:159774/147‑1 (MQ=255)            gCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAAcc                                                                                                                                          >  6:28175/1‑148 (MQ=255)              aaCTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCgaa                                                                                                                                       >  1:248106/1‑149 (MQ=255)                      cgcCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGc                                                                                                                               >  3:56777/1‑149 (MQ=255)                        ccAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGctc                                                                                                                             <  6:60737/149‑1 (MQ=255)                                tCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCAtt                                                                                                                     <  5:355994/149‑1 (MQ=255)                                 cGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCACCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCCCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAAAGGAATGCTCTAGACATTg                                                                                                                    >  1:341879/1‑149 (MQ=255)                                   cTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCg                                                                                                                  <  6:265261/149‑1 (MQ=255)                                          gcgcCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGc                                                                                                           >  7:253373/1‑149 (MQ=255)                                           cgcCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCAtt                                                                                                                     >  8:260448/1‑138 (MQ=255)                                           cgcCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCAtt                                                                                                                     <  7:260448/138‑1 (MQ=255)                                                  aGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAatat                                                                                                   >  4:360879/1‑149 (MQ=255)                                                        aCGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCg                                                                                              >  3:306176/1‑148 (MQ=255)                                                        aCGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATACTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGa                                                                                             <  2:54750/149‑1 (MQ=255)                                                            ccTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGc                                                                                         >  3:220703/1‑149 (MQ=255)                                                             cTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGa                                                                                                       <  1:327430/134‑1 (MQ=255)                                                             cTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGa                                                                                                       >  2:327430/1‑134 (MQ=255)                                                             cTCACCCTGCGATGGCTGCAACGCATCCACTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAAACGAATGGAATGCTATTGTCATTGCGTTCCAGCCGGAATATAATCGaaacc                                                                                         >  1:301066/1‑146 (MQ=255)                                                                aCCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCg                                                                                      <  4:220703/148‑1 (MQ=255)                                                                 cccTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGtt                                                                                    <  2:83229/149‑1 (MQ=255)                                                                      cGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAg                                                                                 <  4:56777/147‑1 (MQ=255)                                                                      cGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAg                                                                                 <  5:28175/147‑1 (MQ=255)                                                                           gCTTCAACGCATCCCGTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACAAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGc                                                                          <  2:341879/149‑1 (MQ=255)                                                                             tGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAgg                                                                                <  6:367937/141‑1 (MQ=255)                                                                             tGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAgg                                                                                >  5:367937/1‑141 (MQ=255)                                                                             tGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTc                                                                        >  2:200987/1‑149 (MQ=255)                                                                               cAACGAATCCCCTCGCAGATGTGTCAGCAAGTAGAAATTACCGCGCAATTCGCCCAGGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGa                                                                      <  2:301066/149‑1 (MQ=255)                                                                                  cGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCgaat                                                                                                                                      <  6:35451/82‑1 (MQ=255)                                                                                  cGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCgaat                                                                                                                                      >  5:35451/1‑82 (MQ=255)                                                                                   gCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCtt                                                                                                                           <  1:146098/92‑1 (MQ=255)                                                                                   gCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCtt                                                                                                                           >  2:146098/1‑92 (MQ=255)                                                                                     aTCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGc                                                                                               >  2:224206/1‑118 (MQ=255)                                                                                     aTCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGc                                                                                               <  1:224206/118‑1 (MQ=255)                                                                                       ccccTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGt                                                              >  5:316163/1‑149 (MQ=255)                                                                                         ccTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGttt                                                            <  3:360879/149‑1 (MQ=255)                                                                                                gATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGAcc                                                     <  8:3859/149‑1 (MQ=255)                                                                                                  tgtgTCAGCAAGTCGAAATTACCGCGCACTTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGt                                                   <  4:94784/149‑1 (MQ=255)                                                                                                   gtgtCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGAc                                                      <  8:192263/145‑1 (MQ=255)                                                                                                   gtgtCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGAc                                                      >  7:192263/1‑145 (MQ=255)                                                                                                      tCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGaa                                                <  5:351816/148‑1 (MQ=255)                                                                                                      tCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAg                                               >  3:11150/1‑149 (MQ=255)                                                                                                          cAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCaa                                           >  3:198331/1‑149 (MQ=255)                                                                                                             gTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGt                                        <  3:348678/149‑1 (MQ=255)                                                                                                              tCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGaa                                                <  1:353283/140‑1 (MQ=255)                                                                                                              tCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGaa                                                >  2:353283/1‑140 (MQ=255)                                                                                                                       ccGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGACTGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGgc                               <  8:253373/148‑1 (MQ=255)                                                                                                                             aaTTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGAc                        <  4:11150/149‑1 (MQ=255)                                                                                                                              aTTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGAcc                       <  7:140582/149‑1 (MQ=255)                                                                                                                              aTTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGAc                        >  8:140582/1‑148 (MQ=255)                                                                                                                                     cACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCgg                >  1:170044/1‑149 (MQ=255)                                                                                                                                      ttggCACTGGTCGGCCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCGGGCGCCAg                          >  1:350965/3‑138 (MQ=255)                                                                                                                                      ttggCACTGGTCGGCCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAg                          <  2:350965/136‑1 (MQ=255)                                                                                                                                      ttggCACTGGTCGGCCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCGgc               >  8:252476/3‑149 (MQ=255)                                                                                                                                        ggCACTGGTCGGCCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCGgcgc             >  7:164573/1‑149 (MQ=255)                                                                                                                                                 cGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGa                                                                      >  5:350077/1‑83 (MQ=255)                                                                                                                                                 cGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGa                                                                      <  6:350077/83‑1 (MQ=255)                                                                                                                                                   aCCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCGGCGCTGGTTATGCCt  >  1:97537/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCGGCGCTGGTTATGCCT  >  NC_002947/225223‑225517

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