Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	133,970	G→C	5.3%	E250Q (GAG→CAG)	PP_0128 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	133,970	0	G	C	5.3%	243.8 / 2.3	75	E250Q (GAG→CAG)	PP_0128	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/34);  new base C (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

ATGAACACCCATGCCACTGACCTGCTGGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGTTTGCGGACCCGGGCATGACTGAAGAC  >  NC_002947/133826‑134116                                                                                                                                                 |                                                                                                                                                   aTGAACACCCATGCCACTGACCTGCTGGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCg                                                                                                                                                >  2:311022/1‑149 (MQ=255) aTGAACACCCATGCCACTGACCTGCTGGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTCCAGCg                                                                                                                                                >  3:42189/1‑149 (MQ=255)     tcacCCATGCCACTGACCTGCTGGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCgtgt                                                                                                                                             <  2:266832/147‑1 (MQ=255)        cccATGCCACTGACCTGCTGGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTCCAGCgtgaggct                                                                                                                                         <  4:42189/149‑6 (MQ=255)          cATGCCACTGACCTGCTGGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGt                                                                                                                                       >  3:169859/1‑149 (MQ=255)           aTGCCACTGACCTGCTGGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTg                                                                                                                                      >  3:260005/1‑149 (MQ=255)            tGCCACTGACCTGCTGGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGc                                                                                                                                     <  8:249378/149‑1 (MQ=255)                 cTGACCTGCTGGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGgcg                                                                                                                                <  5:32399/149‑1 (MQ=255)                    aCCTGCTGGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAg                                                                                                                             >  3:204832/1‑149 (MQ=255)                    aCCTGCTGGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAg                                                                                                                             >  5:224376/1‑149 (MQ=255)                       tgctgGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGgctgcgcca                                                                                                                          >  8:168770/1‑142 (MQ=255)                          tgGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAAtt                                                                                                                       >  4:224435/1‑149 (MQ=255)                            gAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGt                                                                                                                                       >  2:171711/1‑131 (MQ=255)                            gAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGt                                                                                                                                       <  1:171711/131‑1 (MQ=255)                                    gccgcTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGaccacc                                                                                                             >  6:343482/1‑149 (MQ=255)                                    gccgcTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGaccacc                                                                                                             >  5:179495/1‑149 (MQ=255)                                    gccgcTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGaccacc                                                                                                             <  5:343482/149‑1 (MQ=255)                                           cgcgATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCg                                                                                                       <  4:260005/148‑1 (MQ=255)                                           cgcgATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGt                                                                                                      <  6:179495/149‑1 (MQ=255)                                              gATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGc                                                                                                   >  2:145328/1‑149 (MQ=255)                                              gATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGc                                                                                                   >  6:75648/1‑149 (MQ=255)                                                  tGGTCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTACCTCGATCATATTCTGCTCTGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCGAATTTCGGACCACCTTCCCGTTGCTGt                                                                                                <  5:253901/148‑1 (MQ=255)                                                  tGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGt                                                                                                <  2:113567/148‑1 (MQ=255)                                                   gggCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTc                                                                                               <  3:188682/148‑1 (MQ=255)                                                       cTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGAt                                                                                          <  6:342054/149‑1 (MQ=255)                                                          gTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGAtt                                                                                         <  4:106987/147‑1 (MQ=255)                                                          gTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTc                                                                                        >  3:149187/1‑148 (MQ=255)                                                           tCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGAGTTGCCTCGATCATATTCAGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTACCCGATGCTGACGCGATTCGa                                                                                      >  7:315661/1‑149 (MQ=255)                                                            cGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGAt                                                                                     >  1:112618/1‑149 (MQ=255)                                                             gCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGAtt                                                                                    <  6:224376/149‑1 (MQ=255)                                                               cccACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGc                                                                                  <  3:224435/149‑1 (MQ=255)                                                                ccACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGAGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGcc                                                                                 <  4:171793/149‑1 (MQ=255)                                                                  acaAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCt                                                                                >  1:161734/1‑148 (MQ=255)                                                                    aaGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGAt                                                                             <  1:311022/149‑1 (MQ=255)                                                                           ggCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGa                                                                      <  3:26392/149‑1 (MQ=255)                                                                            gCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCATCGATCATATTATGATCAGTAAAAGCCTCACGATCGAGCGTGTCTAGGTGATGTCGCATCCAATATCGGACCACATTCCCTTTGCTGGAGAGATTCGATAGCCTGATGCaattaa                                                                     >  6:224509/1‑144 (MQ=255)                                                                               aCCTTCCCGAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTGGAGGTGCTGGCGCAGGCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCTATTGCCTGATGCATTGACtgt                                                                  <  8:315661/149‑1 (MQ=255)                                                                               aCCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACtgt                                                                  >  8:163737/1‑149 (MQ=255)                                                                               aCCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACtgt                                                                  <  7:163737/149‑1 (MQ=255)                                                                               aCCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACtgt                                                                  <  4:149187/149‑1 (MQ=255)                                                                               aCCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTCCAGCgtgagg                                                                                                                                           >  2:120903/1‑73 (MQ=255)                                                                               aCCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTCCAGCgtgagg                                                                                                                                           <  1:120903/76‑4 (MQ=255)                                                                                 cTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGac                                                                                                                 >  3:107725/1‑100 (MQ=255)                                                                                 cTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGac                                                                                                                 <  4:107725/100‑1 (MQ=255)                                                                                   tCCCCACCTCGCGAACGCAGCGTTCCCTAGATCATAATCTGCTCAGTTCAACCCTCACGCTGGAGCGTGACGAGGTGCTGGAGCAGCCAATTTCGGACCACCATCCCTTTGCTCTCGAGATTCGATAGCCTGATGCATTGACTGTGGa                                                               >  6:165055/1‑148 (MQ=255)                                                                                    ccccAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGa                                                               <  4:204832/147‑1 (MQ=255)                                                                                     cccAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATAc                                                            >  7:30140/1‑149 (MQ=255)                                                                                     cccAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGACTGATGCATTGACTGTGGATAc                                                            >  4:85664/1‑149 (MQ=255)                                                                                      ccAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGAt                                                                             >  5:271980/1‑131 (MQ=255)                                                                                      ccAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGAt                                                                             <  6:271980/131‑1 (MQ=255)                                                                                           tGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGcc                                                      <  4:263568/149‑1 (MQ=255)                                                                                             gCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGaccac                                                                                                              >  8:111275/1‑91 (MQ=255)                                                                                             gCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGaccac                                                                                                              <  7:111275/91‑1 (MQ=255)                                                                                                 ccGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGtttt                                                <  2:161734/149‑1 (MQ=255)                                                                                                   gCAGCGTTTCCTCGATCATATTCTGCTCTGTTAAAGCGTAACGCTGGAGCGTGTCGTGGTGCTGGCGCAGCGAATTTCTGACCACCTTCCCGTTGCTGTCGATATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGa                                              <  5:224509/149‑1 (MQ=255)                                                                                                    cAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAg                                             >  5:93877/1‑149 (MQ=255)                                                                                                      gCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCt                                           <  2:112618/149‑1 (MQ=255)                                                                                                             cTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGcc                                    <  5:75648/149‑1 (MQ=255)                                                                                                                 aTCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCa                                >  5:251008/1‑149 (MQ=255)                                                                                                                  tCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCAc                               >  1:119413/1‑149 (MQ=255)                                                                                                                     tatTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGt                            >  7:317483/1‑149 (MQ=255)                                                                                                                       ttCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGttt                          >  4:269659/1‑149 (MQ=255)                                                                                                                          tGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGTTTGCg                       >  4:194768/1‑149 (MQ=255)                                                                                                                            cTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGTTTGCGGa                     >  1:313427/1‑149 (MQ=255)                                                                                                                             tCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCAGTGACTGTGGATACGCTGCCGGGTTTGAGCTAATTGCCATCACCGTTTGCGGAc                    >  4:188188/1‑149 (MQ=255)                                                                                                                               atttCAAGCCTCACGATTGAGCTTGTGGAGGTGCTGGCGCAGCCTATTTCGGACGAGCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTTCCATCACCTTTTGCGGAc                    <  5:165055/145‑1 (MQ=255)                                                                                                                               aGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGTTTGCGGAccc                  >  8:151099/1‑149 (MQ=255)                                                                                                                                 ttCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGTTTGCGGACCCgg                <  3:269659/149‑1 (MQ=255)                                                                                                                                  tCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGTTTGCGGACCCgg                <  3:194768/148‑1 (MQ=255)                                                                                                                                   cAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGTTTGCGGACCCggg               >  4:200675/1‑148 (MQ=255)                                                                                                                                    aaGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGCCCACCATCCCGTTGCTGTCGAGAATCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGATTTGAGCTAATTGCCATCACCGTTAGCGGACCCGGGCa             >  7:129850/1‑149 (MQ=255)                                                                                                                                    aaGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGTTTGCGGACCCGGGCa             >  2:312876/1‑149 (MQ=255)                                                                                                                                       ccTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGTTTGCGGACCCGGGCATGa          <  6:93877/149‑1 (MQ=255)                                                                                                                                             gCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGTTTGCGGACCCGGGCATGACTg       >  1:224122/1‑146 (MQ=255)                                                                                                                                               tGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGTTTGCGGACCCGGGCATGACTGAAGAc  <  3:85664/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   ATGAACACCCATGCCACTGACCTGCTGGAGCATTCGCCGCTACGCGATTTGGGCCTGGTCGCCCCACAAGTCGAGGCCACCTTCCCCAGCTGGCGACCGCAGCGTTGCCTCGATCATATTCTGCTCAGTTCAAGCCTCACGCTGGAGCGTGTCGAGGTGCTGGCGCAGCCAATTTCGGACCACCTTCCCGTTGCTGTCGAGATTCGATTGCCTGATGCATTGACTGTGGATACGCTGCCGGTTTTGAGCTAATTGCCATCACCGTTTGCGGACCCGGGCATGACTGAAGAC  >  NC_002947/133826‑134116

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