Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,859,102	G→A	5.6%	intergenic (+22/+79)	PP_1661 → / ← PP_1662	dehydrogenase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,859,102	0	G	A	5.6%	220.5 / 2.1	71	intergenic (+22/+79)	PP_1661/PP_1662	dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/32);  new base A (2/2);  total (37/34)
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.45e-01

TTGCCGCGTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTA  >  NC_002947/1858955‑1859250                                                                                                                                                    |                                                                                                                                                      ttGCCGCGTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCTGTGTAGGCCATGTCATGGCAACCCCACGa                                                                                                                                                      <  6:423526/149‑1 (MQ=255)    ccGCGTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGc                                                                                                                                                   >  6:171335/1‑149 (MQ=255)         tGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCgg                                                                                                                                              <  7:43798/149‑1 (MQ=255)              aGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATg                                                                                                                                         >  7:8777/1‑149 (MQ=255)                         ttGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCAttt                                                                                                                              >  2:36196/1‑149 (MQ=255)                             aGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTggg                                                                                                                          >  2:373955/1‑149 (MQ=255)                              gCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTgggg                                                                                                                         <  3:179229/149‑1 (MQ=255)                                   gCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGc                                                                                                                     >  3:211375/1‑148 (MQ=255)                                     aTCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCa                                                                                                                  >  8:298565/1‑149 (MQ=255)                                          tGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAg                                                                                                             >  5:383348/1‑149 (MQ=255)                                          tGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAg                                                                                                             <  1:373955/149‑1 (MQ=255)                                             cgacCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGcc                                                                                                           >  3:176361/1‑148 (MQ=255)                                              gacCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGcc                                                                                                           >  4:208094/1‑147 (MQ=255)                                                ccTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATTCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAgg                                                                                                       >  2:294336/1‑149 (MQ=255)                                                 cTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGc                                                                                                      >  1:38425/1‑149 (MQ=255)                                                 cTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGc                                                                                                      >  1:176027/1‑149 (MQ=255)                                                    gACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTgggg                                                                                                                         >  6:10378/1‑127 (MQ=255)                                                    gACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTgggg                                                                                                                         <  5:10378/127‑1 (MQ=255)                                                    gACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAg                                                                                                   >  1:221102/1‑149 (MQ=255)                                                       gggCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGcc                                                                                                                                       <  2:288995/110‑1 (MQ=255)                                                       gggCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGcc                                                                                                                                       >  1:288995/1‑110 (MQ=255)                                                         gCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAgg                                                                                              >  4:379203/1‑149 (MQ=255)                                                          cGTTACGTGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAggg                                                                                             <  5:405359/149‑1 (MQ=255)                                                          cGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAggg                                                                                             <  5:410986/149‑1 (MQ=255)                                                                  ggCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTcgccg                                                                                     <  2:221102/149‑1 (MQ=255)                                                                  ggCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTcgccg                                                                                     >  4:5045/1‑149 (MQ=255)                                                                       gCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCATTg                                                                                >  7:190402/1‑149 (MQ=255)                                                                         cTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCAAAATGCTGCGGCATTGCCGGAGCATCGTggg                                                                                                                          <  5:88107/105‑1 (MQ=25)                                                                         cTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCAAAATGCTGCGGCATTGCCGGAGCATCGTggg                                                                                                                          >  6:88107/1‑105 (MQ=25)                                                                               gATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGt                                                                       >  5:140559/1‑149 (MQ=255)                                                                                    ccTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAgggg                                                                                            >  3:303671/1‑124 (MQ=255)                                                                                    ccTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAgggg                                                                                            <  4:303671/124‑1 (MQ=255)                                                                                    ccTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCAtt                                                                   <  5:105715/148‑1 (MQ=255)                                                                                    ccTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCAtt                                                                   >  6:105715/1‑148 (MQ=255)                                                                                    ccTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATtg                                                                  >  2:268587/1‑149 (MQ=255)                                                                                     cTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCATTGG‑GCAGGGTCATtgt                                                                 <  8:190402/149‑1 (MQ=255)                                                                                              cTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGgcgcgc                                                        <  6:383348/149‑1 (MQ=255)                                                                                                gTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTg                                                      >  6:138395/1‑149 (MQ=255)                                                                                                 gcGGTAGAGGTGCGTTTCGGGGCGGGGTAGGCCCTGTCATGGCAACCCCACGCTGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTgg                                                     <  7:287735/148‑1 (MQ=255)                                                                                                       gagGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑CCAGGGTCATTGTGTGGCGCGCTGGTTTTcc                                               <  1:36196/149‑1 (MQ=255)                                                                                                          gTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAgg                                            >  4:133963/1‑149 (MQ=255)                                                                                                          gTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAgg                                            >  2:327607/1‑149 (MQ=255)                                                                                                              gTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGATTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCtt                                        <  5:171335/149‑1 (MQ=255)                                                                                                               tttCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTc                                       <  3:5045/149‑1 (MQ=255)                                                                                                                   ggggCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCA‑GGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCt                                  <  2:38425/149‑1 (MQ=255)                                                                                                                     ggCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCt                                  >  4:273368/1‑148 (MQ=255)                                                                                                                         gtgtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGtt                                                   <  7:372828/127‑1 (MQ=255)                                                                                                                         gtgtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAAcg                             >  3:398836/1‑149 (MQ=255)                                                                                                                            tAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAAcgcg                           <  3:379203/148‑1 (MQ=255)                                                                                                                            tAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCgg                          >  6:138109/1‑149 (MQ=255)                                                                                                                              ggCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTcc                                               <  7:49277/126‑1 (MQ=255)                                                                                                                              ggCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTcc                                               >  8:49277/1‑126 (MQ=255)                                                                                                                               gCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTAt                       >  3:245773/1‑149 (MQ=255)                                                                                                                                cccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTAtt                      >  8:114041/1‑149 (MQ=255)                                                                                                                                 ccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCTTTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                     <  4:211375/149‑1 (MQ=255)                                                                                                                                 ccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                     <  8:8777/149‑1 (MQ=255)                                                                                                                                 ccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                     <  4:398836/149‑1 (MQ=255)                                                                                                                                 ccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                     <  2:176027/149‑1 (MQ=255)                                                                                                                                    gTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGTGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgtg                  >  5:98626/1‑149 (MQ=255)                                                                                                                                     tCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgtgt                 <  7:298565/149‑1 (MQ=255)                                                                                                                                     tCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgt                   <  5:138395/147‑1 (MQ=255)                                                                                                                                       tgagCAACCCCAAAATGCTGCGGCATTGCCGGAGCATCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGc               >  3:28417/4‑149 (MQ=37)                                                                                                                                        gagCAACCCCAAAATGCTGCGGCATTGCCGGAGCATCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAg                                    <  4:356054/125‑1 (MQ=25)                                                                                                                                         ggCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTcc                                               <  2:147513/115‑1 (MQ=255)                                                                                                                                         ggCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTcc                                               >  1:147513/1‑115 (MQ=255)                                                                                                                                           cAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGTTCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTc            <  1:327607/149‑1 (MQ=255)                                                                                                                                                ccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc      <  4:245773/149‑1 (MQ=255)                                                                                                                                                 cACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAAt                                                                                 >  4:100934/1‑74 (MQ=255)                                                                                                                                                 cACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAAt                                                                                 <  3:100934/74‑1 (MQ=255)                                                                                                                                                 cACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAg                                             >  7:357800/1‑109 (MQ=255)                                                                                                                                                 cACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAg                                             <  8:357800/109‑1 (MQ=255)                                                                                                                                                    gATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTa  <  3:208094/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TTGCCGCGTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGG‑GCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTA  >  NC_002947/1858955‑1859250

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