Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,178,806	G→C	6.4%	S16R (AGC→AGG)	PP_4555 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,178,806	0	G	C	6.4%	249.0 / 7.2	78	S16R (AGC→AGG)	PP_4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/42);  new base C (3/2);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

CGTCTGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGAC  >  NC_002947/5178658‑5178952                                                                                                                                                      |                                                                                                                                                   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                                                                                                                                   <  3:92346/149‑5 (MQ=255)                  tgtCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTTGCACAGcctgcc                                                                                                                                   <  1:246819/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                            <  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                  >  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:351142/149‑1 (MQ=255)                                                                                                                                                    aGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGAc  <  7:346113/149‑1 (MQ=255)                                                                                                                                                      gCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGAc  >  2:137081/1‑147 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CGTCTGTCAGGCTCATGTGTCCGTTGCGCATGCGTACCGTGGTGGCCGTCTGCACCGGCTTCAGCGC‑CTCGGCCGCCAGCATGCGCCCGGAACCGTGAATGGCTTCGGCGTAGGCCATGCATTCGGCGTCTGCCGGGCTGTCGGGCAGGCTGTGCAACAGCCCCTCGTCACAATAGACCAGGCACAGGTATTTCATGGTCGCCTCCAGGGTATTGGCAAAGTGCGTTTGGCGATAACGGCTCAGGGCTTGAGATCAAACAGCGCTTTTTGCGTTTCCATGTCGAAAGGCGCCGAC  >  NC_002947/5178658‑5178952

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