Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,380,722	G→T	7.7%	intergenic (‑110/+25)	PP_t67 ← / ← PP_t66	tRNA‑Met/tRNA‑Met

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,380,722	0	G	T	7.7%	194.4 / 3.0	63	intergenic (‑110/+25)	PP_t67/PP_t66	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base T (2/3);  total (29/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 = 1.00e+00

TGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATG  >  NC_002947/1380577‑1380859                                                                                                                                                    |                                                                                                                                          tGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGtt                                                                                                                                         <  7:22099/148‑1 (MQ=255) tGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGGTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGtt                                                                                                                                         <  6:20536/148‑1 (MQ=255)   aGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTaa                                                                                                                                       <  7:158952/148‑1 (MQ=255)    gCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACa                                                                                                                                     <  1:163276/149‑1 (MQ=255)        gacgaGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTgg                                                                                                                                 >  5:213857/1‑149 (MQ=255)          cgaGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGcc                                                                                                                               <  7:51459/149‑1 (MQ=255)            aGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGCCtt                                                                                                                             <  2:80256/149‑1 (MQ=255)              cTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCtttt                                                                                                                           <  5:251879/149‑1 (MQ=255)                 cctggCTGCTCCATCCCGCGTCAGTGGGGAGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTc                                                                                                                          <  1:285967/144‑1 (MQ=255)                 ccAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCtt                                                                                                                        >  5:148922/1‑149 (MQ=255)                  cAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCAACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCtt                                                                                                                        >  3:117318/1‑148 (MQ=255)                   aGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGCCtttttcttt                                                                                                                      <  7:185482/149‑6 (MQ=18)                       tgctCCCTCCCGCGTCAGTGGGGAGGATTCTACAGTGTGTGGTTTGTATGTCAAGCCGTTTGCGCTTTGATTATTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGCCtttttctttgc                                                                                                                    <  1:122357/149‑8 (MQ=25)                            cATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTg                                                                                                             >  8:109634/1‑149 (MQ=255)                                       aGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCgg                                                                                                          <  3:101004/141‑1 (MQ=255)                                       aGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCgg                                                                                                          >  4:101004/1‑141 (MQ=255)                                        gTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGAt                                                                                                 <  8:179178/149‑1 (MQ=255)                                        gTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGAt                                                                                                 >  7:179178/1‑149 (MQ=255)                                             gCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGc                                                                                                            >  3:356688/1‑133 (MQ=255)                                             gCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGc                                                                                                            <  4:356688/133‑1 (MQ=255)                                                  ttCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACcgac                                                                                       <  3:159929/149‑1 (MQ=255)                                                   tCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACcgac                                                                                       >  8:104138/1‑148 (MQ=255)                                                     tACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCt                                                                                                                              >  2:261902/1‑107 (MQ=255)                                                     tACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCt                                                                                                                              <  1:261902/107‑1 (MQ=255)                                                     tACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCgacgac                                                                                    >  6:118017/1‑149 (MQ=255)                                                       cAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCgacgac                                                                                    <  6:148922/147‑1 (MQ=255)                                                          tgtgtgGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCg                                                                               <  7:104138/149‑1 (MQ=255)                                                            tgtgGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCggg                                                                             >  4:150300/1‑149 (MQ=255)                                                            tgtgGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCgg                                                                              <  7:109634/148‑1 (MQ=255)                                                             gtgGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGt                                                                            >  4:174626/1‑149 (MQ=255)                                                                gTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATAAGCAGGCAAAGAAAAAGGCCACTGTTTTAACAGTGGCCTTTTCTTCAATCTGGTTGCGGGAGCCGGATTTGAAACGACGACATTCGGGTTAt                                                                         >  4:180678/1‑149 (MQ=17)                                                                     tATGTTTAG‑CGCTTG‑GGTTTGATTTTTCGACGTTTTTTCGGTTTTGGTTTTCAGCTTTGAAATTAAAATGGGAGTTTGTTAATTGTGTGATTTTCTTGAATATTTTTGCTGGAGCGTGATTTGAACCGACGACCTTCGGTTTTTGAGc                                                                     <  6:213857/148‑1 (MQ=255)                                                                      aTGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGccc                                                                   >  4:41088/1‑149 (MQ=255)                                                                           aaG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCgg                                                                                                          >  5:215594/1‑105 (MQ=255)                                                                           aaG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCgg                                                                                                          <  6:215594/105‑1 (MQ=255)                                                                           aaG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcgacga                                                              <  3:150300/149‑1 (MQ=255)                                                                                    g‑gCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCACGACGAGCTACCAg                                                      >  7:290504/1‑149 (MQ=255)                                                                                      gCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAgg                                                     >  2:371141/1‑149 (MQ=255)                                                                                        tttGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGgct                                                   >  5:330467/1‑149 (MQ=255)                                                                                        tttGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGgct                                                   >  6:202143/1‑149 (MQ=255)                                                                                            aTTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTTGCCTTTTCTTGAATCTGTTTGCGGGAGCCGGATTTGAACAGACGACCTTCGGGTTATGAGCCCGACGAGGTACCAGGCTGCTc                                               <  4:117318/149‑1 (MQ=255)                                                                                                  cGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATccc                                         >  5:145231/1‑149 (MQ=255)                                                                                                   gACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCcg                                        >  7:57407/1‑149 (MQ=255)                                                                                                   gACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCcg                                        >  6:210473/1‑149 (MQ=255)                                                                                                    aCGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTc                                                                                >  5:168558/1‑108 (MQ=255)                                                                                                    aCGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTc                                                                                <  6:168558/108‑1 (MQ=255)                                                                                                    aCGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCcgc                                       <  8:57407/149‑1 (MQ=255)                                                                                                       ttttttCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGt                                     >  2:350259/1‑148 (MQ=255)                                                                                                         ttttCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCa                                   <  1:371141/148‑1 (MQ=255)                                                                                                             cgcgAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTgg                               <  6:330467/148‑1 (MQ=255)                                                                                                              gcgAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTgggg                             >  3:19858/1‑149 (MQ=255)                                                                                                                     cTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTc                      <  5:118017/149‑1 (MQ=255)                                                                                                                           agcagGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAg                 <  8:290504/148‑1 (MQ=255)                                                                                                                           agcagGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGGTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGc                <  5:202143/149‑1 (MQ=255)                                                                                                                                 cAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAg          <  3:41088/149‑1 (MQ=255)                                                                                                                                    aGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATACCGCGTCAGTGGGGCGGATTCTACAGCGTCTAGtt        >  5:164342/1‑148 (MQ=255)                                                                                                                                     gAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTta      <  1:350259/149‑1 (MQ=255)                                                                                                                                      gaaaaGGCCACTGTTAACACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTta      >  2:88916/2‑148 (MQ=255)                                                                                                                                      aaacaGGCCACTGGCTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGACTCTACAGCGTGTAGTTtat     <  4:19858/149‑1 (MQ=255)                                                                                                                                      aaaaaGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTtat     <  3:174626/149‑1 (MQ=255)                                                                                                                                          aGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATg  >  7:66749/1‑148 (MQ=255)                                                                                                                                              cACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTcc                                              <  7:291256/100‑1 (MQ=37)                                                                                                                                              cACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTcc                                              >  8:291256/1‑100 (MQ=37)                                                                                                                                                    |                                                                                                                                          TGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAG‑CGCTAG‑GCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATG  >  NC_002947/1380577‑1380859

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑