Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,178,810	T→C	6.6%	H15R (CAC→CGC) ‡	PP_4555 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,178,810	0	T	C	6.6%	213.7 / 4.6	75	H15R (CAC→CGC) ‡	PP_4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/40);  new base C (3/2);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

TGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGACCAGTGT  >  NC_002947/5178662‑5178958                                                                                                                                                      |                                                                                                                                                     tGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCtgt                                                                                                                                                      >  5:257014/1‑149 (MQ=255)  gTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCtgtg                                                                                                                                                     >  3:324032/1‑149 (MQ=255)   tCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCtgtg                                                                                                                                                     <  8:201670/148‑1 (MQ=255)   tCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCtgt                                                                                                                                                      <  1:325644/147‑1 (MQ=255)   tCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGc                                                                                                                                                    <  8:272040/149‑1 (MQ=255)       gCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACa                                                                                                                                                <  8:9539/149‑1 (MQ=255)       gCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCAGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACa                                                                                                                                                >  4:29273/1‑149 (MQ=255)           aTGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGc                                                                                                                                              <  5:293247/147‑1 (MQ=255)             gtgtCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCt                                                                                                                                          >  5:285513/1‑149 (MQ=255)              tgtCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGCTCTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCt                                                                                                                                          >  1:96912/1‑149 (MQ=255)              tgtCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTTGCACAGcctgcc                                                                                                                                         <  1:246819/149‑5 (MQ=255)              tgtCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTTGCACAGcctgcc                                                                                                                                         <  3:92346/149‑5 (MQ=255)              tgtCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTc                                                                                                                                         >  5:1618/1‑149 (MQ=255)                tCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGt                                                                                                                                       <  6:1618/149‑1 (MQ=255)                 ccGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGt                                                                                                                                       >  2:182867/1‑148 (MQ=255)                          aTGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGAc                                                                                                                             <  1:177231/149‑1 (MQ=255)                                aCCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGca                                                                                                                       >  1:359869/1‑149 (MQ=255)                                  cGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGcaca                                                                                                                     <  2:55334/149‑1 (MQ=255)                                   gtggtgGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGTACAg                                                                                                                    <  4:320706/149‑1 (MQ=255)                                   gtggtgGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAg                                                                                                                    <  2:5940/149‑1 (MQ=255)                                          ccGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTc                                                                                                             <  6:215788/149‑1 (MQ=255)                                           cGTCTGCACCGGCTTCAGCGA‑CTCGGCCGCCAGCATCCGCCCGGAACCGAGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCACGTCACAATAGACCAGGCACAGGTATTTCa                                                                                                            >  3:43592/1‑149 (MQ=255)                                             tCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATg                                                                                                          <  4:324032/149‑1 (MQ=255)                                             tCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATg                                                                                                          <  6:285513/149‑1 (MQ=255)                                                    cGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCt                                                                                                   >  3:300909/1‑149 (MQ=255)                                                       cTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTc                                                                                                  >  5:133393/1‑147 (MQ=255)                                                       cTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCa                                                                                                <  6:133393/149‑1 (MQ=255)                                                        ttCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAg                                                                                               >  8:150916/1‑149 (MQ=255)                                                          cAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAgg                                                                                              <  6:257014/148‑1 (MQ=255)                                                             cgc‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTAt                                                                                          <  2:31156/149‑1 (MQ=255)                                                              gc‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTAtt                                                                                         >  2:184162/1‑149 (MQ=255)                                                                 cTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTc                                                                                                             <  3:4439/127‑1 (MQ=255)                                                                 cTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTc                                                                                                             >  4:4439/1‑127 (MQ=255)                                                                     gccgccAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCaaa                                                                                   <  5:312055/149‑1 (MQ=255)                                                                      ccgccAGCATGCGCCCGGATCCGTGAATGGCATCGGCGTAGGGCATGCATTCGGCGTCTCCCGGGCAGTCGGGCAGGCTGTGCAACAGCCCCACGTCACAATAGACCAGGCACAGGTATTTCAAGGTCGCCTCCAGCGAATTGGCAAAg                                                                                  >  3:190664/1‑149 (MQ=255)                                                                      ccgccAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAg                                                                                  >  2:103891/1‑149 (MQ=255)                                                                            gCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGtt                                                                            >  4:324548/1‑149 (MQ=255)                                                                              aTGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTCTCGGGCAGGCTGTGCAACAGCCCATCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTAGTGGCAAAGTGCGTTTg                                                                          >  6:301586/1‑149 (MQ=255)                                                                               tGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTgg                                                                         <  5:23110/149‑1 (MQ=255)                                                                                  gcCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAg                                                                                                                    <  1:2336/103‑1 (MQ=255)                                                                                  gcCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAg                                                                                                                    >  2:2336/1‑103 (MQ=255)                                                                                  gcCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGa                                                                      >  3:77839/1‑149 (MQ=255)                                                                                    ccGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGAt                                                                     >  3:336766/1‑148 (MQ=255)                                                                                    ccGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTCGCGAt                                                                     >  5:371822/1‑148 (MQ=255)                                                                                        aaCCGTGAATGGATTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTTTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCTCCTCCAGGGTATTTGCAAAGTGCGTTTGGCGATAACgg                                                                <  1:350484/149‑1 (MQ=255)                                                                                            gTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCa                                                            <  7:150916/149‑1 (MQ=255)                                                                                                       cGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTgaga                                                  <  4:300909/148‑1 (MQ=255)                                                                                                          cGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCaa                                              <  2:96912/149‑1 (MQ=255)                                                                                                                ccATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAgcg                                        >  7:68845/1‑149 (MQ=255)                                                                                                                 cATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCtt                                                      >  3:99027/1‑134 (MQ=255)                                                                                                                 cATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCtt                                                      <  4:99027/134‑1 (MQ=255)                                                                                                                 cATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAgcgc                                       >  2:158397/1‑149 (MQ=255)                                                                                                                     cATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCttt                                    <  1:309306/148‑1 (MQ=255)                                                                                                                      aTTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCTATAACGGCTCAGGGCTTGAGATCAAACAGCGCttttt                                  <  1:182867/149‑1 (MQ=255)                                                                                                                      aTTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCttttt                                  >  4:199319/1‑149 (MQ=255)                                                                                                                      aTTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCttttt                                  <  2:359869/149‑1 (MQ=255)                                                                                                                       ttCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTg                                 <  4:77839/149‑1 (MQ=255)                                                                                                                         cGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGc                                <  3:29273/148‑1 (MQ=255)                                                                                                                         cGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGc                                <  3:324548/148‑1 (MQ=255)                                                                                                                          ggCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTAGTCTCAATAGACCAGGCAGTGGTATTTCATGGTCGCCTCCAGGGTACTGGCAAAGTGCGTTTGGCGATAACGGCTCTGGGCTTGAGATCAAACAGCGCTTTTTGCGt                              <  4:43592/149‑1 (MQ=255)                                                                                                                            cGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGttt                            >  7:108043/1‑149 (MQ=255)                                                                                                                              tCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTcc                          <  3:199319/149‑1 (MQ=255)                                                                                                                              tCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGCGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCATGGCTTGAGATCAAACAGCGCTTTTTGCGTTTcc                          <  6:371822/149‑1 (MQ=255)                                                                                                                                cgaggggCTGTTGCACAGCCTGCCCGACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCAt                        >  6:210511/5‑149 (MQ=255)                                                                                                                                 gaggggCTGTTGCACAGCCTGCCCGACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCAt                        <  3:198733/145‑1 (MQ=255)                                                                                                                                 gaggggCTGTTGCACAGCCTGCCCGACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCAt                        >  4:319604/4‑148 (MQ=255)                                                                                                                                 gaggggCTGTTGCACAGCCTGCCCGACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCAt                        >  4:198733/4‑148 (MQ=255)                                                                                                                                 gaggggCTGTTGCACAGCCTGCCCGACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCAt                        <  3:319604/145‑1 (MQ=255)                                                                                                                                 gCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTAtt                                                                                         <  4:42430/83‑1 (MQ=255)                                                                                                                                 gCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTAtt                                                                                         >  3:42430/1‑83 (MQ=255)                                                                                                                                  ccGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATg                       <  1:103891/148‑1 (MQ=255)                                                                                                                                   cGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTc                     >  8:346113/1‑149 (MQ=255)                                                                                                                                   cGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTc                     >  8:351142/1‑149 (MQ=255)                                                                                                                                       cTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGaaa                 <  4:336766/149‑1 (MQ=255)                                                                                                                                            gggCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGaaa                 >  6:96015/1‑144 (MQ=255)                                                                                                                                            gggCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGaaa                 <  5:96015/144‑1 (MQ=255)                                                                                                                                                aGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGAc        <  7:346113/149‑1 (MQ=255)                                                                                                                                                aGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGAc        <  7:351142/149‑1 (MQ=255)                                                                                                                                                  gCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGAc        >  2:137081/1‑147 (MQ=255)                                                                                                                                                    tgtgCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGACCAgt    >  5:130371/1‑149 (MQ=255)                                                                                                                                                     gtgCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGAc        <  7:30629/144‑1 (MQ=255)                                                                                                                                                     gtgCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGAc        >  8:30629/1‑144 (MQ=255)                                                                                                                                                      tgCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTTTTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTCCGTTTCCATGTCGAAAGGCGCCGACCAgggt  >  2:82694/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGACCAGTGT  >  NC_002947/5178662‑5178958

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