BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,380,731	A→C	5.5%	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,731	0	A	C	5.5%	413.4 / 11.2	127	intergenic (‑113/+22)	PP_t67/PP_t66	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand): ref base A (59/60); new base C (2/5); total (61/66)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 128 total aligned reads displayed.

TGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTCAAGCGCT  >  NC_002947/1380583‑1380874

tGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTaa                                                                                                                                                <  2:4439469/150‑1 (MQ=255)

 gAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAAc                                                                                                                                               >  2:2953051/1‑150 (MQ=255)

  aGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGTGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCTCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACa                                                                                                                                              <  1:2005352/150‑1 (MQ=255)

     ccGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTg                                                                                                                                           <  2:2717249/150‑1 (MQ=255)

      cgacgaGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTgg                                                                                                                                          <  2:5343553/150‑1 (MQ=255)

       gacgaGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAG‑AAAAGGCCACTGTTAACACAGTGGcc                                                                                                                                        <  1:4782056/150‑1 (MQ=18)

        acgaGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGcc                                                                                                                                        <  1:2341521/150‑1 (MQ=255)

           aGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGCCttt                                                                                                                                     <  1:3983409/150‑1 (MQ=255)

           aGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGCCttt                                                                                                                                     <  1:3983444/150‑1 (MQ=255)

             cTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGCATGTCAAGCGCTAGGCTTTGATTTTTCGAAGTTTTTTCGAGAATGCTTATAAGCTGGCAAAGAAATAGGACACGCTGGTAAGAGTGGCCttttg                                                                                                                                   >  2:3755129/1‑149 (MQ=17)

              tACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCt                                                                                                                                  <  2:2518453/150‑1 (MQ=255)

                  aGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGaa                                                                                                                              <  2:585362/150‑1 (MQ=255)

                  aGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGaa                                                                                                                              >  2:4278408/1‑150 (MQ=255)

                    gctgctCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATc                                                                                                                            >  1:4733928/1‑150 (MQ=255)

                    gctgctCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATc                                                                                                                            >  1:3346867/1‑150 (MQ=255)

                      tgctCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTg                                                                                                                          <  1:3953594/150‑1 (MQ=255)

                       gctCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTgg                                                                                                                         <  1:2542550/150‑1 (MQ=255)

                       gctCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTgg                                                                                                                         <  1:5107392/150‑1 (MQ=255)

                        ctCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGt                                                                                                                        >  2:4669068/1‑150 (MQ=255)

                         tCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGAGGTTTTTTAGCGAATGATTAACATAAGGCAAAGAAAAAGGACACTGTGTTAACAGTGGCCTTTACTTGAATCTGGtt                                                                                                                       >  1:2118486/1‑150 (MQ=18)

                         tCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGCCtttttcttt                                                                                                                               <  1:5217871/142‑6 (MQ=18)

                         tCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTTAACACAGTGGCCtttttcttt                                                                                                                               >  2:5217871/1‑137 (MQ=18)

                         tCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGtt                                                                                                                       >  2:3252418/1‑150 (MQ=255)

                           cATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGc                                                                                                                     <  1:4147822/150‑1 (MQ=255)

                           cATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGc                                                                                                                     >  1:3712962/1‑150 (MQ=255)

                            ctCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCg                                                                                                                    <  2:3080090/149‑1 (MQ=255)

                            aTCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAATGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCg                                                                                                                    >  1:815570/1‑150 (MQ=255)

                            aTCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCg                                                                                                                    <  1:2869466/150‑1 (MQ=255)

                            aTCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCg                                                                                                                    >  1:2469709/1‑150 (MQ=255)

                               ccGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGATTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGa                                                                                                                 <  2:3352120/150‑1 (MQ=255)

                               ccGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTGTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGATTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGa                                                                                                                 <  2:3352259/150‑1 (MQ=255)

                                  cgTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGcc                                                                                                              <  1:4658028/150‑1 (MQ=255)

                                    tCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCCGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCgg                                                                                                            <  2:3016829/150‑1 (MQ=255)

                                        tGGGGCTGATTCTACACTGTGTGTTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGATTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGAttt                                                                                                        <  2:3816442/150‑1 (MQ=255)

                                        tGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGAttt                                                                                                        <  1:4638796/150‑1 (MQ=255)

                                         ggggCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTg                                                                                                       <  1:67987/150‑1 (MQ=255)

                                         ggggCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTg                                                                                                       >  1:3318958/1‑150 (MQ=255)

                                          gggCGGATTCTACAGTGTGTGTTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGACACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGa                                                                                                      <  2:3972099/150‑1 (MQ=255)

                                          gggCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTgg                                                                                                      <  1:2953051/150‑2 (MQ=255)

                                          gggCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGa                                                                                                      <  1:4591235/150‑1 (MQ=255)

                                          gggCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGa                                                                                                      >  2:2734479/1‑150 (MQ=255)

                                           ggCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGaa                                                                                                     <  1:4504623/150‑1 (MQ=255)

                                           ggCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGaa                                                                                                     <  1:727378/150‑1 (MQ=255)

                                           ggCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGaa                                                                                                     <  1:5348606/150‑1 (MQ=255)

                                            gCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAAc                                                                                                    <  2:4322144/150‑1 (MQ=255)

                                            gCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATAAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAAc                                                                                                    >  2:5323442/1‑150 (MQ=255)

                                                aTTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACcgac                                                                                                <  1:4278408/150‑1 (MQ=255)

                                                  tCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACcgacga                                                                                              <  2:3681023/150‑1 (MQ=255)

                                                    tACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACa                                                                                                                                              >  1:3541000/1‑100 (MQ=255)

                                                    tACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACa                                                                                                                                              <  2:3541000/100‑1 (MQ=255)

                                                    tACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGAcc                                                                                            >  1:1359816/1‑150 (MQ=255)

                                                    tACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGAcc                                                                                            >  1:1895601/1‑150 (MQ=255)

                                                    tACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATAAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTACTTGAATCAGGTTGCGGTAGCAGGATTTGAACAGACGAcc                                                                                            >  2:2021033/1‑150 (MQ=18)

                                                     aCAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGTCCACTGTGTTAACAGTGGCCTTTTCTTGAATATGGTTGCGGGAGCCGGATTTGAACCGACGACCt                                                                                           >  2:1882465/1‑150 (MQ=255)

                                                      cAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCtt                                                                                          >  1:3419447/1‑150 (MQ=255)

                                                        gtgtgtGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCg                                                                                        <  1:511223/150‑1 (MQ=255)

                                                        gtgtgtGGTTTGTATGTCAAGCGCTAGGATTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGAATTTGAACCGACGACCTTCg                                                                                        <  1:2989129/150‑1 (MQ=255)

                                                           tgtgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGt                                                                                     <  2:3318958/150‑1 (MQ=255)

                                                           tgtgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGt                                                                                     <  2:3319007/150‑1 (MQ=255)

                                                            gtgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGATTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGtt                                                                                    <  2:1935089/150‑1 (MQ=255)

                                                             tgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCtt                                                                                                                                 >  2:5292283/1‑104 (MQ=255)

                                                             tgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCtt                                                                                                                                 <  1:5292283/104‑1 (MQ=255)

                                                              ggTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTAt                                                                                  >  1:3566926/1‑150 (MQ=255)

                                                                tttGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGa                                                                                <  1:4287771/150‑1 (MQ=255)

                                                                 ttGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAg                                                                               >  1:1841838/1‑150 (MQ=255)

                                                                    tATGTCAAGCGCTAGGCTTTGATTTTTCGCCGTTTTTTCGCGACTGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGccc                                                                            <  2:1553363/150‑1 (MQ=255)

                                                                     aTGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcg                                                                           <  1:2734479/150‑1 (MQ=255)

                                                                     aTGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCATTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcg                                                                           >  2:1494139/1‑150 (MQ=255)

                                                                      tGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTTTCTGTTGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcga                                                                          <  1:1724374/150‑1 (MQ=255)

                                                                      tGTCAAGCGCTAGGCTTTGATTGTTCGACGTTTTTTCGCTAATGGTTATCAGCAGGCAAAGTAATTGGCCACTGTGTTAACTGTTGCCTTTTCTCGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcga                                                                          <  2:2412268/150‑1 (MQ=17)

                                                                        tCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcgacg                                                                        <  2:1935057/150‑1 (MQ=255)

                                                                        tCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcgacg                                                                        <  2:4674753/150‑1 (MQ=255)

                                                                        tCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGAGGGAGCCGGATTTGAACCGACTACCTTCGGGTTATGAGCCcgacg                                                                        >  1:57615/1‑150 (MQ=255)

                                                                         cAAGCGCTAGTCTTTGATTTTTCGACGTTTTTTCGCGAATTCTTATCAGCAGGAAAAGAAAAAGGCCACTGTGTTAACAGTGTACTTTTCTAGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcgacga                                                                       <  2:4387485/150‑1 (MQ=17)

                                                                         cAAGCGCGAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcgacga                                                                       <  2:3683851/150‑1 (MQ=255)

                                                                                tAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCa                                                                >  1:1087125/1‑150 (MQ=255)

                                                                                        gATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTc                                                        >  2:5298467/1‑150 (MQ=255)

                                                                                             ttCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGTGTTAGGAGCCCGACGAGGTAAAAGGCTGCTccctac                                                   >  1:2333698/1‑146 (MQ=17)

                                                                                                gACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTAGAATCTGGTTGCGGGAGCCGGATATGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCCTCCcgc                                                >  1:5327909/1‑150 (MQ=255)

                                                                                                    ttttttCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCa                                            >  1:1200834/1‑150 (MQ=255)

                                                                                                      ttttCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGt                                          >  2:3728092/1‑150 (MQ=255)

                                                                                                        ttCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTgg                                        <  2:2082104/150‑1 (MQ=255)

                                                                                                                   ttATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCgg                                                                                       <  2:1845093/92‑1 (MQ=255)

                                                                                                                   ttATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCgg                                                                                       >  1:1845093/1‑92 (MQ=255)

                                                                                                                   ttATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTa                             >  2:2264027/1‑150 (MQ=255)

                                                                                                                   ttATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTa                             >  1:4491134/1‑150 (MQ=255)

                                                                                                                            ggCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACcgacg                                                                                               >  1:1774648/1‑75 (MQ=39)

                                                                                                                               aaaGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGAttt                                                                                                        >  2:712406/1‑63 (MQ=255)

                                                                                                                               aaaGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGAttt                                                                                                        >  2:712312/1‑63 (MQ=255)

                                                                                                                               aaaGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAgg                                                              <  2:1057007/105‑1 (MQ=35)

                                                                                                                               aaaGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAgg                                                              >  1:1057007/1‑105 (MQ=35)

                                                                                                                                  gAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcgac                                                                         >  1:2895910/1‑91 (MQ=37)

                                                                                                                                    aaaaGGCCACTGTTAACACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCcgcgc                                              <  1:287829/116‑2 (MQ=11)

                                                                                                                                    aaaaGGCCACTGTTAACACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCcgcgc                                              >  2:287829/1‑115 (MQ=11)

                                                                                                                                    aaaaGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGATGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTatat            <  1:3252418/150‑1 (MQ=255)

                                                                                                                                         gCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTAt                                                                                  >  1:2341840/1‑75 (MQ=39)

                                                                                                                                           cACTGTGTTAACAGTGGCCTTTTCTGGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGc                                                                     >  2:1882764/1‑86 (MQ=18)

                                                                                                                                                 gtTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTAt                                                                                  >  1:1618888/1‑67 (MQ=35)

                                                                                                                                                 gtTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTAt                                                                                  <  2:1618888/67‑1 (MQ=35)

                                                                                                                                                 gtTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCgt                       >  1:5237006/1‑126 (MQ=255)

TGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTAACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTCAAGCGCT  >  NC_002947/1380583‑1380874

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`