BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,380,725	A→C	6.1%	intergenic (‑113/+22)	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,725	0	A	C	6.1%	309.7 / 9.9	99	intergenic (‑113/+22)	PP_t67/PP_t66	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand): ref base A (46/47); new base C (3/3); total (49/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

Only 100 of 102 total aligned reads displayed.

AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTCAAGCGCTAGG  >  NC_002947/1380579‑1380871

aGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGTCCACTGTGTTaa                                                                                                                                                     >  3:36528/1‑148 (MQ=255)

   ccGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAg                                                                                                                                                  <  5:357212/148‑1 (MQ=255)

    cgacgaGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTg                                                                                                                                                <  2:32263/149‑1 (MQ=255)

    cgacgaGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTg                                                                                                                                                <  4:266035/149‑1 (MQ=255)

    cgacgaGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTg                                                                                                                                                <  2:26191/149‑1 (MQ=255)

    cgacgaGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTg                                                                                                                                                <  1:306702/149‑1 (MQ=255)

       cgtgCTACCAGGCGTCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCATCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGcc                                                                                                                                             <  4:36528/146‑1 (MQ=255)

       cgaGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGcc                                                                                                                                             <  6:32104/149‑1 (MQ=255)

         aGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCtt                                                                                                                                           <  1:25847/149‑1 (MQ=255)

           cTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCtttt                                                                                                                                         <  4:349955/149‑1 (MQ=255)

             aCCTGGCTGGACCTTCCGGCTTCCGTGGGGGGGTTTCTACTGTGTGTGGTTTGTATGTCAAGCGCTTGGCTTTTATTTTTCTACGTTTTTTGGCTATTGCTTATCAGCAGGCAAATAAAACGCCCACTGTGTAAACAGGGGCCTTTTc                                                                                                                                        <  2:90564/148‑1 (MQ=255)

              ccAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTc                                                                                                                                        <  7:72303/147‑1 (MQ=255)

              ccAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCtt                                                                                                                                      <  5:203019/149‑1 (MQ=255)

                 ggCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑aa                                                                                                                                  <  8:46396/149‑1 (MQ=255)

                 ggCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑aa                                                                                                                                  >  6:82350/1‑149 (MQ=255)

                      ctCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCt                                                                                                                               <  5:126238/147‑1 (MQ=255)

                        ccATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCt                                                                                                                                       <  5:290088/138‑1 (MQ=255)

                        ccATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCt                                                                                                                                       >  6:290088/1‑138 (MQ=255)

                        ccATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGtt                                                                                                                           >  4:190939/1‑149 (MQ=255)

                         cATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTgg                                                                                                                             <  7:94278/146‑1 (MQ=255)

                         cATCCCGCGTCAGTGCGGCTGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGCCGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGtt                                                                                                                           <  3:48009/148‑1 (MQ=255)

                           tCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCg                                                                                                                        >  3:313565/1‑149 (MQ=255)

                           tCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCg                                                                                                                        >  3:2415/1‑149 (MQ=255)

                            cccGCGTCAGTCGGGCGGATTCTACAGTGTCTGGCCTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAACGAAAAAGGCCACTGTGTTAACTGTGGCCTTTTCTTG‑AATCTGGTTGCgg                                                                                                                       <  8:346659/149‑1 (MQ=18)

                             ccGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCgg                                                                                                                       >  6:187856/1‑148 (MQ=255)

                             ccGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTAATCAGCAGGCAAAGAAAAAGGCCACTCTGTTAACAGTGGCCTATTCTTG‑AATCTGGATGCggg                                                                                                                      >  2:199970/1‑149 (MQ=255)

                              cgcgTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGCCtttttcttt                                                                                                                                    <  2:262428/135‑6 (MQ=18)

                              cgcgTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGCCtttttcttt                                                                                                                                    >  1:262428/1‑130 (MQ=18)

                              cgcgTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCgg                                                                                                                       >  8:25823/1‑147 (MQ=255)

                                     gTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGAt                                                                                                              >  5:197515/1‑149 (MQ=255)

                                     gTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGAt                                                                                                              <  2:28194/149‑1 (MQ=255)

                                     gTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCACGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGAt                                                                                                              >  6:639/1‑149 (MQ=255)

                                        gggCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCACGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTg                                                                                                           <  5:639/149‑1 (MQ=255)

                                        gggCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTAATTATTAGACGTTTTTTAGAAAATGATTATAATCACTCAAAGAAAAAGTCAACTTTGTTAACAGTGGACTTTTCTTT‑AATCTGGTTGCGGGAGCCGGAttt                                                                                                            >  4:198218/1‑148 (MQ=17)

                                           cGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAAc                                                                                                        <  4:2415/149‑1 (MQ=255)

                                           cGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAAc                                                                                                        >  2:284454/1‑149 (MQ=255)

                                           cGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCAATGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAAc                                                                                                        >  1:238203/1‑149 (MQ=255)

                                                  tACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTTTCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACATTTGACTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCgccgac                                                                                                 >  3:319836/1‑149 (MQ=18)

                                                  tACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGc                                                                                                                                              <  8:361474/105‑1 (MQ=21)

                                                  tACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCgacgac                                                                                                 >  3:213455/1‑149 (MQ=255)

                                                    cAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCt                                                                                               >  5:241947/1‑149 (MQ=255)

                                                         tgtgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCgg                                                                                           <  6:241947/148‑1 (MQ=255)

                                                         tgtgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCgg                                                                                           <  4:313565/148‑1 (MQ=255)

                                                         tgtgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCgg                                                                                           <  3:190939/148‑1 (MQ=255)

                                                          gtgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGt                                                                                         >  3:338555/1‑149 (MQ=255)

                                                              tttGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGCAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTAt                                                                                      <  1:284454/149‑1 (MQ=255)

                                                                  tATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGt                                                                                         >  4:14666/1‑141 (MQ=255)

                                                                  tATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGt                                                                                         <  3:14666/141‑1 (MQ=255)

                                                                  tATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGc                                                                                 <  5:82350/148‑1 (MQ=255)

                                                                  tATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTAAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGcc                                                                                <  5:187856/149‑1 (MQ=255)

                                                                        aaGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCcgacg                                                                           <  6:197515/148‑1 (MQ=255)

                                                                        aaGCGCTAGGCTTTCATTATTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCcgacg                                                                           <  1:156276/148‑1 (MQ=255)

                                                                        aaGCGCTAGGCGTTGATTTTTCGACTTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGGCACTGTGTTAACAGTGGCCTTTTGTTG‑AATCTGGTTGCGGTAGCCTGATTTGAACCGACGACCTTCGGGTTATG‑AGCCcgacga                                                                          <  1:199970/149‑1 (MQ=17)

                                                                                   tttGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGTCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCACGACGAGCTACCAGgct                                                               >  5:55917/1‑149 (MQ=255)

                                                                                   tttGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCgg                                                                                           <  1:194955/122‑1 (MQ=255)

                                                                                   tttGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCgg                                                                                           >  2:194955/1‑122 (MQ=255)

                                                                                   tttGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATg                                                                                     <  1:326127/128‑1 (MQ=255)

                                                                                   tttGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATg                                                                                     >  2:326127/1‑128 (MQ=255)

                                                                                   tttGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGgct                                                               >  5:306993/1‑149 (MQ=255)

                                                                                     tGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGgctgc                                                             >  2:335821/1‑149 (MQ=255)

                                                                                      gATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAgg                                                                 >  3:337737/1‑144 (MQ=255)

                                                                                      gATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAgg                                                                 <  4:337737/144‑1 (MQ=255)

                                                                                        tttttCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑a                                                                                   <  4:369099/124‑1 (MQ=255)

                                                                                        tttttCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑a                                                                                   >  3:369099/1‑124 (MQ=255)

                                                                                             cGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATcc                                                      <  4:213455/148‑1 (MQ=255)

                                                                                              gACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCcg                                                    >  7:40847/1‑149 (MQ=255)

                                                                                              gACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGAGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCACGACGAGCTACCAGGCTGCTCCATCCcg                                                    >  3:237769/1‑149 (MQ=255)

                                                                                                 gTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGt                                                 <  7:25823/149‑1 (MQ=255)

                                                                                                        cgcgAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTGTCTTG‑AATCTGGTTCCGGGAGCCGGAGTTGAACCGAAGACCTTCGGGTTATG‑CGCCCGACGAGCTACCAGGCTGCTCCGTCCCGCGTCCGTggg                                          >  5:387182/1‑149 (MQ=17)

                                                                                                         gcgAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTgggg                                         >  2:379331/1‑149 (MQ=255)

                                                                                                             aTGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGa                                     <  3:224447/149‑1 (MQ=255)

                                                                                                             aTGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGa                                     <  6:306993/149‑1 (MQ=255)

                                                                                                               gCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCgacgac                                                                                                 >  1:166084/1‑88 (MQ=255)

                                                                                                               gCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCgacgac                                                                                                 <  2:166084/88‑1 (MQ=255)

                                                                                                               gCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTc                                                                                             <  5:24684/92‑1 (MQ=255)

                                                                                                               gCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTc                                                                                             >  6:24684/1‑92 (MQ=255)

                                                                                                                 ttATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTc                                  >  7:366916/1‑148 (MQ=255)

                                                                                                                      agcagGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAAGCCCGACGAGCTACCAGGCTGCTCCAt                                                        <  1:223660/122‑1 (MQ=255)

                                                                                                                      agcagGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCAt                                                        >  2:223660/1‑121 (MQ=39)

                                                                                                                       gcagGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCg                           >  8:39820/1‑149 (MQ=255)

                                                                                                                         agGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGTATTTGAACCGACGACCTTCGGGTTTTG‑AGCCCGACGATCTACCAGGCTGCTACATCCCGCGTCAGTGGGGCGGATTCTACAGCgtg                         <  2:238203/149‑1 (MQ=255)

                                                                                                                         agGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCgtg                         >  7:27595/1‑149 (MQ=255)

                                                                                                                           gCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCgtg                         >  3:93543/1‑147 (MQ=255)

                                                                                                                              aaGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGtt                    >  3:85171/1‑149 (MQ=255)

                                                                                                                                gAAAATGGCCACTGTGTTAACTGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTta                  <  4:237769/149‑1 (MQ=255)

                                                                                                                                 aaaaaGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTta                  >  3:323237/1‑148 (MQ=255)

                                                                                                                                   aaaGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCACGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTatat               >  7:108034/1‑149 (MQ=255)

                                                                                                                                    aaGGCCACTGTTAACACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATg              >  5:128352/1‑149 (MQ=255)

                                                                                                                                    aaGGCCACTGTGTTAACAGTGGCCTTTTGTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATg              <  4:338555/149‑1 (MQ=255)

                                                                                                                                    aaGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGc                            <  7:127160/135‑1 (MQ=255)

                                                                                                                                    aaGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGc                            >  8:127160/1‑135 (MQ=255)

                                                                                                                                      ggCCACTGTTAACACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGt             >  2:197149/1‑148 (MQ=255)

                                                                                                                                       gCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGCCGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTc            >  7:139409/1‑148 (MQ=255)

                                                                                                                                         cACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGc                                                                        >  4:85847/1‑86 (MQ=38)

                                                                                                                                         cACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGc                                                                        <  3:85847/86‑1 (MQ=38)

                                                                                                                                          aCTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTtat                 <  2:232789/140‑1 (MQ=255)

                                                                                                                                          aCTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTtat                 >  1:232789/1‑140 (MQ=255)

                                                                                                                                                ttAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTCAAGCGCTAgg  >  2:388954/1‑149 (MQ=255)

                                                                                                                                                  aaCAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTgggg                                         >  6:35691/1‑108 (MQ=255)

                                                                                                                                                  aaCAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTgggg                                         <  5:35691/108‑1 (MQ=255)

AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTG‑AATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATG‑AGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTCAAGCGCTAGG  >  NC_002947/1380579‑1380871

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