BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,123,426	C→A	5.1%	noncoding (380/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,426	0	C	A	5.1%	426.1 / 8.1	128	noncoding (380/398 nt)	PP_mr13	cspA regulatory region
Reads supporting (aligned to +/- strand): ref base C (56/65); new base A (5/2); total (61/67)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

GCACAGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAAGT  >  NC_002947/1123277‑1123574

gCACAGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCAc                                                                                                                                                      >  2:4453563/1‑150 (MQ=255)

 cacaGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACt                                                                                                                                                     >  1:676318/1‑150 (MQ=255)

    aGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTgg                                                                                                                                                  >  2:5040626/1‑150 (MQ=255)

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

          aaGTCATAGACCATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAgggg                                                                                                                                            <  2:2428922/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:373272/150‑1 (MQ=255)

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

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

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

                   ttcATCTCTGCGAGGCGCCGGC‑CCAGGGCGACCGGTGGAAAATCCCCGTGGGGAGTGAGTAGAGGGTTGGACATGCAGGACGCTCCAAAAGTCGAGCCCACTATTCTACGGGCAATATTTTGCCCGTGCACTTGGCAGGGGCTTTTCaaa                                                                                                                                   >  2:4207687/3‑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)

                      aTCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTc                                                                                                                                >  2:816626/1‑150 (MQ=255)

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

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

                         gCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATCCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTgg                                                                                                                             <  1:1598166/150‑1 (MQ=255)

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

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

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

                             ccAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCt                                                                                                                         <  1:2422702/150‑1 (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)

                                           cAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCACAGTCTGGCGCTGGATCAGACAGGt                                                                                                            <  2:2878503/150‑1 (MQ=255)

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

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

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

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

                                                         gAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGc                                                                                              >  2:2640363/1‑150 (MQ=255)

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

                                                          aaaaTCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGcc                                                                                             <  2:2172127/150‑1 (MQ=255)

                                                           aaaTCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAACAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCt                                                                                            <  2:3810654/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                   ggCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTg                                                                    <  1:2254519/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                                       tCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAATAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTgg                                                >  2:4753567/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:1753342/1‑150 (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)

                                                                                                                       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)

                                                                                                                                 gggCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                      <  2:1497749/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)

                                                                                                                                   actttgaaaaGCCCCTGCCAAGTGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGc                    >  2:3353991/7‑150 (MQ=39)

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

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

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

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

                                                                                                                                     aaaaaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCt                  <  2:4546382/150‑1 (MQ=255)

                                                                                                                                      aaaaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCtt                 <  2:3093711/150‑1 (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:340847/1‑150 (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:1043454/150‑1 (MQ=255)

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

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

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

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

GCACAGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGC‑CCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAAGT  >  NC_002947/1123277‑1123574

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