Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,178,813	A→G	6.7%	L14S (TTG→TCG)	PP_4555 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,178,813	0	A	G	6.7%	220.3 / 4.6	78	L14S (TTG→TCG)	PP_4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/38);  new base G (3/2);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

GCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGACCAGTGT  >  NC_002947/5178668‑5178958                                                                                                                                                   |                                                                                                                                                  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                                                                                                          <  6:285513/149‑1 (MQ=255)                                       tCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATg                                                                                                          <  4:324032/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                                                                                                             >  4:4439/1‑127 (MQ=255)                                                           cTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTc                                                                                                             <  3:4439/127‑1 (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                                  <  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                        <  3:319604/145‑1 (MQ=255)                                                                                                                           gaggggCTGTTGCACAGCCTGCCCGACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCAt                        >  4:198733/4‑148 (MQ=255)                                                                                                                           gaggggCTGTTGCACAGCCTGCCCGACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCAt                        <  3:198733/145‑1 (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)                                                                                                                                      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)                                                                                                                                                   |                                                                                                                                                  GCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGACCAGTGT  >  NC_002947/5178668‑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: ‑