Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,178,801	G→C	6.7%	P18R (CCC→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,801	0	G	C	6.7%	216.4 / 6.9	75	P18R (CCC→CGC) ‡	PP_4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (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 = 3.04e-01

GGGCCGTCTGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAA  >  NC_002947/5178654‑5178943                                                                                                                                                     |                                                                                                                                               gggCCGTCTGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCggg                                                                                                                                               <  4:87598/149‑1 (MQ=255) gggCCGTCTGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCggg                                                                                                                                               >  6:23110/1‑149 (MQ=255)     cGTCTGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAgg                                                                                                                                           <  2:224135/149‑1 (MQ=255)         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                                       <  4:99027/134‑1 (MQ=255)                                                                                                                         cATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCtt                                       >  3:99027/1‑134 (MQ=255)                                                                                                                         cATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAgcgc                        >  2:158397/1‑149 (MQ=255)                                                                                                                             cATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCttt                     <  1:309306/148‑1 (MQ=255)                                                                                                                              aTTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCTATAACGGCTCAGGGCTTGAGATCAAACAGCGCttttt                   <  1:182867/149‑1 (MQ=255)                                                                                                                              aTTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCttttt                   <  2:359869/149‑1 (MQ=255)                                                                                                                              aTTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCttttt                   >  4:199319/1‑149 (MQ=255)                                                                                                                               ttCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTg                  <  4:77839/149‑1 (MQ=255)                                                                                                                                 cGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGc                 <  3:324548/148‑1 (MQ=255)                                                                                                                                 cGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGc                 <  3:29273/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         >  4:198733/4‑148 (MQ=255)                                                                                                                                         gaggggCTGTTGCACAGCCTGCCCGACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCAt         >  4:319604/4‑148 (MQ=255)                                                                                                                                         gaggggCTGTTGCACAGCCTGCCCGACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCAt         <  3:198733/145‑1 (MQ=255)                                                                                                                                         gaggggCTGTTGCACAGCCTGCCCGACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCAt         <  3:319604/145‑1 (MQ=255)                                                                                                                                         gCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTAtt                                                                          >  3:42430/1‑83 (MQ=255)                                                                                                                                         gCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTAtt                                                                          <  4:42430/83‑1 (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)                                                                                                                                                     |                                                                                                                                               GGGCCGTCTGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAA  >  NC_002947/5178654‑5178943

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