Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,178,799	C→T	6.5%	D19N (GAC→AAC)	PP_4555 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,178,799	0	C	T	6.5%	227.4 / 4.0	77	D19N (GAC→AAC)	PP_4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/41);  new base T (3/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

AACGGGCCGTCTGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAA  >  NC_002947/5178651‑5178943                                                                                                                                                      |                                                                                                                                                 aaCGGGCCGTCTGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTc                                                                                                                                                  <  3:351394/149‑1 (MQ=255)    gggCCGTCTGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCggg                                                                                                                                               >  6:23110/1‑149 (MQ=255)    gggCCGTCTGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCggg                                                                                                                                               <  4:87598/149‑1 (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: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         <  3:319604/145‑1 (MQ=255)                                                                                                                                            gaggggCTGTTGCACAGCCTGCCCGACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCAt         >  4:198733/4‑148 (MQ=255)                                                                                                                                            gaggggCTGTTGCACAGCCTGCCCGACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCAt         >  4:319604/4‑148 (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)                                                                                                                                                      |                                                                                                                                                 AACGGGCCGTCTGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAA  >  NC_002947/5178651‑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: ‑