BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,123,428	T→G	5.5%	noncoding (378/398 nt)	PP_mr13 ←	cspA regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,123,428	0	T	G	5.5%	426.7 / 10.3	127	noncoding (378/398 nt)	PP_mr13	cspA regulatory region
Reads supporting (aligned to +/- strand): ref base T (55/64); new base G (5/2); total (60/67)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 135 total aligned reads displayed.

AGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAAGTGAA  >  NC_002947/1123281‑1123577

    aaaaGTCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAgg                                                                                                                                                 <  1:1618824/150‑1 (MQ=255)

    aaaaGTCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAgg                                                                                                                                                 <  2:58592/150‑1 (MQ=255)

      aaGTCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAgggg                                                                                                                                               <  2:2428953/150‑1 (MQ=255)

        gTCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCt                                                                                                                                             <  2:1536040/150‑1 (MQ=255)

         tCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACGTGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCtt                                                                                                                                            >  2:1103748/1‑150 (MQ=255)

         tCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGGGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCtt                                                                                                                                            >  1:1779553/1‑150 (MQ=255)

         tCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCtt                                                                                                                                            <  2:3237410/150‑1 (MQ=255)

           aTAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCCAAGTGCAGGGGCtttt                                                                                                                                          <  1:487804/150‑1 (MQ=255)

           aTAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCtttt                                                                                                                                          >  1:4546382/1‑150 (MQ=255)

             aGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCa                                                                                                                                        >  1:4324812/1‑150 (MQ=255)

              gAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCaa                                                                                                                                       >  2:4687882/1‑150 (MQ=255)

              gAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCaa                                                                                                                                       >  2:1530211/1‑150 (MQ=255)

               aaCATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCaaa                                                                                                                                      <  1:4684740/150‑1 (MQ=255)

               aaCATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCaaa                                                                                                                                      <  1:4825998/150‑1 (MQ=255)

               aaCATCGCTGCCAGGCGCCGGC‑CCAGGCCCACCCGTGCAAAATCCCCCTGGGGACTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTCTTCTACGGGCAAAAATAAGCCCCTGCACTTGGCAGGGGCTTGTCaaa                                                                                                                                      >  2:2624883/1‑150 (MQ=255)

                   tCGCTGCCCGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCt                                                                                                                                  <  2:4066115/150‑1 (MQ=255)

                    cGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAATCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTg                                                                                                                                 <  1:240225/150‑1 (MQ=255)

                       tGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGGCGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGgcg                                                                                                                              >  1:2249506/1‑150 (MQ=255)

                       tGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGgcg                                                                                                                              >  1:2249331/1‑150 (MQ=255)

                        ggcctgggcCGGCGCCTGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGgcg                                                                                                                              >  1:1043454/8‑150 (MQ=255)

                        gCCAGGCGCCGCC‑CCAGCCCGACCAGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCACGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTAGGCAGGGGCTTTTCAAAGTCTGgcgc                                                                                                                             <  2:2707108/150‑1 (MQ=255)

                          cAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTg                                                                                                                           >  1:1287391/1‑150 (MQ=255)

                              cgcCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATc                                                                                                                       <  1:4723947/150‑1 (MQ=255)

                                   gC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGACTGGCGCTGGATCAGACa                                                                                                                  <  2:756133/150‑1 (MQ=255)

                                      ccAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAgg                                                                                                                >  1:1098300/1‑150 (MQ=255)

                                      ccAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAgg                                                                                                                <  1:3132744/150‑1 (MQ=255)

                                       cAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGt                                                                                                               >  1:2280431/1‑150 (MQ=255)

                                               ccGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCACAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCtt                                                                                                       <  2:2824189/150‑1 (MQ=255)

                                               ccGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCtt                                                                                                       <  1:771391/150‑1 (MQ=255)

                                               ccGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCtt                                                                                                       <  2:2146391/150‑1 (MQ=255)

                                                  gTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTc                                                                                                    >  1:2399034/1‑150 (MQ=255)

                                                   tGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCa                                                                                                   >  2:774414/1‑150 (MQ=255)

                                                     gAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGc                                                                                                 <  2:1858177/150‑1 (MQ=255)

                                                      aaacTCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGcc                                                                                                <  2:365834/150‑1 (MQ=255)

                                                       aaaTCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCt                                                                                               <  2:3810478/150‑1 (MQ=255)

                                                           cccccTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCAt                                                                                           <  2:236759/150‑1 (MQ=255)

                                                            ccccTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATg                                                                                          >  1:5082586/1‑150 (MQ=255)

                                                             cccTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGc                                                                                         <  1:248554/150‑1 (MQ=255)

                                                             cccTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGc                                                                                         >  1:4522682/1‑150 (MQ=255)

                                                             cccTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGc                                                                                         >  2:2412357/1‑150 (MQ=255)

                                                              ccTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGcc                                                                                        <  1:3599290/150‑1 (MQ=255)

                                                               cagcggAGAGACTACAGCCTTCCACATCCACGACCCACCAAAACACCACCCCACTATTCTACCGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTACAAAGACTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCt                                                                                       <  2:4863638/146‑1 (MQ=255)

                                                                  gggAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTc                                                                                    <  2:843763/150‑1 (MQ=255)

                                                                         agtAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCtt                                                                             >  2:2862876/1‑150 (MQ=255)

                                                                                gCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGagctg                                                                      >  2:2553667/1‑146 (MQ=255)

                                                                                gCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGc                                                                      >  2:1730761/1‑150 (MQ=255)

                                                                                  ttGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACTGGCAAAAAAAAGCCCCTGCACTTGTAAGGGGCTTTTCAAAGTCTGTCTCTGGATCAGACAGGTTGAACCTTTTCAGCCTGCATGACTTTATGGCCTTGGACTGCta                                                                    >  2:3837715/1‑149 (MQ=255)

                                                                                  ttGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCtt                                                                    >  2:377320/1‑150 (MQ=255)

                                                                                      aCATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGaa                                                                <  2:3157326/150‑1 (MQ=255)

                                                                                      aCATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGaa                                                                >  1:775029/1‑150 (MQ=255)

                                                                                      aCATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGaa                                                                <  2:1360328/150‑1 (MQ=255)

                                                                                         tGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGt                                                             <  2:4336273/150‑1 (MQ=255)

                                                                                                 gCTCCAAAAGACGAGCCCACTATTCTACGGGCACAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCt                                                     >  1:31961/1‑150 (MQ=255)

                                                                                                  cTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTg                                                    <  1:4695569/150‑1 (MQ=255)

                                                                                                     cAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCCGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGcc                                                 >  2:4330816/1‑150 (MQ=255)

                                                                                                     cAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGcc                                                 >  1:1497749/1‑150 (MQ=255)

                                                                                                     cAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACGTGGCAGGGGCTTTTCAAAGGCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGTCTGCTTCGAAGGTGACTTGCTGGcc                                                 >  2:1451678/1‑150 (MQ=255)

                                                                                                      aaaaGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCct                                                >  2:4900639/1‑150 (MQ=255)

                                                                                                           aCGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCttt                                           <  2:80581/150‑1 (MQ=255)

                                                                                                            cGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGAGCTGGATCAGACAGGTTGAACGTTGTCAGCCTGGATGGCTTTCCGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTc                                          <  2:889009/150‑1 (MQ=255)

                                                                                                             gAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCa                                         <  2:139062/150‑1 (MQ=255)

                                                                                                             gAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCa                                         <  2:4074749/150‑1 (MQ=255)

                                                                                                             gAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCa                                         <  2:3624930/150‑1 (MQ=255)

                                                                                                                cccACTATTCTACGGGCAAAACAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGc                                      >  1:1100093/1‑150 (MQ=255)

                                                                                                                 ccACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGct                                     >  2:2932805/1‑150 (MQ=255)

                                                                                                                  cACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCGTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTGCAGGctc                                    >  1:3485095/1‑150 (MQ=255)

                                                                                                                   aCTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGAc                                                          <  2:4882092/127‑1 (MQ=255)

                                                                                                                   aCTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGAc                                                          >  1:4882092/1‑127 (MQ=255)

                                                                                                                   aCTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGctct                                   <  1:2067610/150‑1 (MQ=255)

                                                                                                                       ttCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGaa                               <  1:4188146/150‑1 (MQ=255)

                                                                                                                        tCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAg                              >  2:4735855/1‑150 (MQ=255)

                                                                                                                        tCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAg                              >  2:4735799/1‑150 (MQ=255)

                                                                                                                             gggCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                         <  1:2175646/150‑1 (MQ=255)

                                                                                                                              gactttgaaaaGCCCCTGCCAAGTGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTg                        >  2:1090682/8‑150 (MQ=39)

                                                                                                                               gCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGc                       <  2:2368365/150‑1 (MQ=255)

                                                                                                                               gCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGc                       <  1:2685061/150‑1 (MQ=255)

                                                                                                                                cAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGcc                      >  2:3057011/1‑150 (MQ=255)

                                                                                                                                 tttgaaaaGCCCCTGCCAAGTGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTCGCTGGCCTTCTTTCa                                         >  2:612563/5‑130 (MQ=38)

                                                                                                                                 tttgaaaaGCCCCTGCCAAGTGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTCGCTGGCCTTCTTTCa                                         <  1:612563/126‑1 (MQ=38)

                                                                                                                                 aaaaaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCt                     >  1:2852824/1‑150 (MQ=255)

                                                                                                                                  aaaaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCtt                    <  2:3093818/150‑1 (MQ=255)

                                                                                                                                   aaaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTc                   <  1:4453563/150‑1 (MQ=255)

                                                                                                                                     aaaaGCCCCTGCCAAGTGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGa                 >  1:4703252/1‑150 (MQ=39)

                                                                                                                                     aaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGa                 <  2:3558746/150‑1 (MQ=255)

                                                                                                                                     aaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGa                 <  1:741176/150‑1 (MQ=255)

                                                                                                                                           cccTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAgcgc           >  1:340848/1‑150 (MQ=255)

                                                                                                                                           cccTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAgcgc           >  1:1377413/1‑150 (MQ=255)

                                                                                                                                            ccTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGTCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGcgcg          >  1:3211578/1‑150 (MQ=255)

                                                                                                                                            ccTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGcgcg          >  2:435069/1‑150 (MQ=255)

                                                                                                                                            ccTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGcgcg          >  2:203093/1‑150 (MQ=255)

                                                                                                                                            ccTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGcgcg          >  2:3015622/1‑150 (MQ=255)

                                                                                                                                             cTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCgg         >  1:3034365/1‑150 (MQ=255)

                                                                                                                                              tGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGa        >  2:623773/1‑150 (MQ=255)

                                                                                                                                                 aCTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAAGt     <  1:5089987/150‑1 (MQ=255)

                                                                                                                                                   ttGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAAGTGa   <  1:5086112/150‑1 (MQ=255)

                                                                                                                                                   ttGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAAGTGa   <  1:623773/150‑1 (MQ=255)

                                                                                                                                                   ttGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCATCGATAGCGCGGAAGTGa   >  1:2454732/1‑150 (MQ=255)

                                                                                                                                                    tGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAAGTGaa  <  2:2274905/150‑1 (MQ=255)

                                                                                                                                                    tGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCGGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAAGTGaa  <  2:2274781/150‑1 (MQ=255)

                                                                                                                                                    tGGCAGGGGCTTTACAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTACAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAAGTGaa  <  2:1098300/150‑1 (MQ=255)

AGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAAGTGAA  >  NC_002947/1123281‑1123577

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