BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,286,965	A→T	7.6%	intergenic (+439/+101)	narI → / ← rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,286,965	0	A	T	7.6%	740.0 / 3.0	286	intergenic (+439/+101)	narI/rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand): ref base A (97/144); new base T (6/17); total (113/173)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

CTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGT  >  NC_000913/1286816‑1287114

      cTTCCCCGCTTGGGCAGAATATTTGATTACGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCg                                                                                                                                                 >  1:1051445/1‑150 (MQ=255)

       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                                <  2:8330709/150‑1 (MQ=255)

         ccccGCTTGGGCAGAATTTTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                              <  2:3361842/150‑1 (MQ=255)

          cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGATATTTCCGGGGCTGTTGAGAGGGAGGGTGTTGGGGGAAGGATAAATCATCGCTGAGATCTTCGACCTTCGGGTCTGTGACTGCGGCAACGCTCTCTCGCTGTCGCTCGAggc                                                                                                                                             >  2:1455830/1‑148 (MQ=14)

          cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                                                             >  1:785878/1‑150 (MQ=255)

           ccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                                                            >  1:8979189/1‑150 (MQ=255)

           ccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                                                            >  1:2928252/1‑150 (MQ=255)

            cGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACACACGCAGACGCACCAGCCGa                                                                                                                                           >  1:6137202/1‑150 (MQ=14)

                    cAGAGTTTTTGATTGGGTCTGAGTTTGAGTATTCGGTGTATTATGTAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGt                                                                                                                                   <  2:5791462/150‑1 (MQ=14)

                    cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGt                                                                                                                                   >  2:2652040/1‑150 (MQ=12)

                    cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGt                                                                                                                                   <  1:4158308/150‑1 (MQ=255)

                    cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCACACTCGCTGGCGCGCGAGACGAACCATAGt                                                                                                                                   >  2:4060861/1‑150 (MQ=11)

                         tatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                              <  1:6216816/150‑1 (MQ=255)

                          atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGc                                                                                                                             <  1:5742018/150‑1 (MQ=255)

                               aTTGCGGATTCGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTtctc                                                                                                                        <  1:5254462/150‑1 (MQ=255)

                                ttGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                                                                                                       <  1:8160577/150‑1 (MQ=255)

                                ttGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCCCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCATAGTCGAATCTTCTCa                                                                                                                       >  2:5010531/1‑150 (MQ=255)

                                 tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                                                                                                      <  2:448471/150‑1 (MQ=255)

                                  gCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACACACCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCACTCACGCAGCCCCTCGAGCCGAACCTGAGTCGGAGCTTATAATc                                                                                                                     >  1:687117/1‑150 (MQ=14)

                                   gggATTAGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATcc                                                                                                                    <  2:6466769/149‑1 (MQ=255)

                                          gtgtGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCCGCCACGCACGCTCGCTGACCCGCGAGCCGAACCGGCGTCGAAGCATCGCCTCCTGCCCCg                                                                                                             >  2:7950997/4‑150 (MQ=14)

                                          gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCg                                                                                                             <  2:1101632/150‑1 (MQ=255)

                                          gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCACCAGCCGAACCTTACTCGAAGCTTCTCATCCTTCCCCg                                                                                                             >  2:6631492/1‑150 (MQ=255)

                                           tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGCTCTGCGCGCTTCGCCCTTCGGGTCGTTTCATCAGGCAACGCACTCTCGCTGACGCAAGAGCCGAACCTAAGTCGAAGCTTATCATACTTCCCCGc                                                                                                            >  2:3056803/1‑150 (MQ=14)

                                             tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                                          <  1:1861791/150‑1 (MQ=255)

                                             tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCACTCACCCAGACCCAAGAGCCGAAACTTAGTCGAAGATTATAATCCTTACCCGCaa                                                                                                          >  1:8249455/1‑149 (MQ=14)

                                              gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATg                                                                                                         >  2:9115677/1‑150 (MQ=255)

                                                gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGAGCATCGCCCGACTGGTAGTGACCGGAGGAAACGCGATCACGCTGACGCTCGAGACGCACCGTTGTCTAAGCTTCTAATCCTTCACCgaaatgg                                                                                                       >  2:227991/1‑144 (MQ=9)

                                                   ttCCGGGGCTTTTGAAAGTGATGGTGGTGGGTGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAg                                                                                                    <  1:2230579/150‑1 (MQ=11)

                                                     ccGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                                  <  2:4006926/150‑1 (MQ=40)

                                                       ggggCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                                                <  1:3249504/150‑1 (MQ=255)

                                                       ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCATCTCAACCTTCCCCGCAAGGGCAGAata                                                                                                >  1:5366606/1‑150 (MQ=255)

                                                          gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                             <  2:8979189/150‑1 (MQ=38)

                                                           cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAACCACACTAGCAGACGCTCGAGCCGAACCCTAGTAGAAGCTTATCATACATACCCCCAGGGGCAGTAGAtttt                                                                                            >  1:7323346/1‑149 (MQ=11)

                                                            ttttGCAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                                           <  1:1256471/150‑1 (MQ=255)

                                                            ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                                           <  1:1498097/150‑1 (MQ=38)

                                                                  aaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCgg                                                                                     <  1:8196703/150‑1 (MQ=39)

                                                                   agcgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                                                    <  2:2995749/147‑1 (MQ=255)

                                                                   agcgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                                                    <  2:6881189/147‑1 (MQ=255)

                                                                   aGTGATGGTGGTGGGGGAAGGATTACTCACCGCTGCGCGCATCGCCCTACGGGTCGTTGCAAGCGCCAAAGCTCTCACGCTGACGCAAGAGCCGAACCGGTGTCGAAGCTTCACATCCTTCCCCGCATGGGCACACTAGTTGAGTGCggg                                                                                    >  2:8786816/1‑149 (MQ=11)

                                                                      gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGAGCTTCGCCCTTCGGGTCGTTGCCTGCCGCAACGCACTCTCGCTGACCCTCGAGGCGAAACTTGGTAGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAATTGAATTCGGATTc                                                                                 >  2:4832866/1‑150 (MQ=11)

                                                                       aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGa                                                                                                                                                       >  2:3919505/1‑79 (MQ=31)

                                                                         ggttgtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCt                                                                              <  1:2734125/150‑1 (MQ=255)

                                                                          gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                             <  2:5414917/150‑1 (MQ=255)

                                                                          ggggtggGCGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                             <  1:8747482/148‑1 (MQ=255)

                                                                            ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                           <  1:8004636/150‑1 (MQ=255)

                                                                            ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                           <  1:5039362/150‑1 (MQ=255)

                                                                            ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                           <  1:6024295/150‑1 (MQ=255)

                                                                                   gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                    <  1:8809206/150‑1 (MQ=255)

                                                                                   gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                    <  1:4976544/150‑1 (MQ=255)

                                                                                    aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                   <  2:7730680/150‑1 (MQ=255)

                                                                                     aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                  <  1:9115677/150‑1 (MQ=255)

                                                                                     aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                  <  2:4356258/150‑1 (MQ=255)

                                                                                        aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                                >  2:6081596/1‑69 (MQ=35)

                                                                                             tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGAGAGTTGAGGGACTTTTg                                                          >  2:6092509/1‑150 (MQ=255)

                                                                                             tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTg                                                          >  2:7269949/1‑150 (MQ=255)

                                                                                               aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                             >  2:3156421/1‑113 (MQ=39)

                                                                                               aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCTCACACGCTGACGCTCGAGCCGAACCGTTGTCGAAGCGTCTCAACCTTACCCGCAGGGGCAGAAGATTTGATTGCGGAGTCGGGTGGCCGTTCAGGGACTTTTgta                                                        >  2:4155401/1‑148 (MQ=11)

                                                                                                gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                       <  1:4508754/150‑1 (MQ=255)

                                                                                                gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                       <  1:77562/150‑1 (MQ=255)

                                                                                                    tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                   <  2:8573933/150‑1 (MQ=255)

                                                                                                     gcgcgcTTCGCGCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAt                                                  <  1:439072/150‑1 (MQ=255)

                                                                                                      cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtg                                                 <  1:931900/150‑1 (MQ=255)

                                                                                                         gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtg                                              <  2:8745222/150‑1 (MQ=255)

                                                                                                            tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTATCATCCTTCCCCGCATGGGCAGAAAATTTGATTTCGGATTCGGTTGGGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                           >  2:3757603/1‑150 (MQ=14)

                                                                                                             cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                          <  2:2198114/150‑1 (MQ=255)

                                                                                                             cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                          <  1:7972702/150‑1 (MQ=255)

                                                                                                             cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                          <  1:2769105/150‑1 (MQ=255)

                                                                                                                   tcgggtcgggGCCTGCGGCAACGCTCTCGCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCGTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                    <  1:2505379/150‑1 (MQ=255)

                                                                                                                   tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                    <  1:71831/150‑1 (MQ=255)

                                                                                                                        tCGTTGCCTGCGGCTACGCTCTCTAGCTGACGCTCGAGCCGATCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                               <  1:2577383/150‑1 (MQ=255)

                                                                                                                        tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                               <  1:3074195/150‑1 (MQ=255)

                                                                                                                        tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                               >  2:2579108/1‑150 (MQ=255)

                                                                                                                         cGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGCTGGTGGTGGGGGACGGAttcg                              >  2:6420182/1‑148 (MQ=14)

                                                                                                                          gTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                             <  1:4050528/150‑1 (MQ=255)

                                                                                                                          gTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAttcga                             >  2:8429144/1‑147 (MQ=255)

                                                                                                                          gTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                             <  1:1948193/150‑1 (MQ=255)

                                                                                                                          gTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                             <  1:6631492/150‑1 (MQ=255)

                                                                                                                          gTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                             >  2:2921962/1‑150 (MQ=255)

                                                                                                                          gTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                             >  2:5965441/1‑150 (MQ=255)

                                                                                                                           ttGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAGTCGCTGGCGAGTGCAGGGCCTTTGGCAAGTGAGGGTGGTGGGGGACGGATTACTc                            >  2:4195848/1‑150 (MQ=14)

                                                                                                                           ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                            >  1:6911854/1‑150 (MQ=255)

                                                                                                                              ccTGCGGCAACGCTCACACGCTGACGCACGAGCCGAACCTTAGTCGAAGCTTCACAAACATCCCCCCAAGGGCAGAATATTTGATTTCGGATTAGCTTTAGAGTTAAGGGACTTTTGGAAGTTATTGTGGTGGGGGGAGGATTAATCAgc                         >  1:390436/1‑150 (MQ=255)

                                                                                                                               cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGTACATTAGTCGAAGCTTCTCTTCCGTCCCGGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcg                        <  1:5858973/150‑1 (MQ=14)

                                                                                                                               cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcg                        >  2:9063519/1‑150 (MQ=255)

                                                                                                                                  cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTg                     >  2:2618995/1‑150 (MQ=255)

                                                                                                                                  cGGCAACGCTCTCTAGCTGACGCACGAGCCGAACCTTAGTCGAAGCATCTCATCCATCCCCCCATGGGCAGAATATTTGATTGCGGATTCGCGTGAGAGTTGAGGGACTTTTGTAAGTGATGGTGGTGGGGGAAGGATTTCTCAGCGCTg                     >  1:7650555/1‑150 (MQ=17)

                                                                                                                                   ggCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgc                    <  2:8032333/150‑1 (MQ=255)

                                                                                                                                      aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTATTCGTCGCTtcgct                 >  2:5251179/1‑149 (MQ=255)

                                                                                                                                      aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcg                 >  2:1705083/1‑150 (MQ=255)

                                                                                                                                      aaCGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCTCAACCTACCCCGCAAGGGCCGAGTATTTGTTTGCGGTTTCGATTGAGAGTTCAGGGACTTTTGTAAGTGATGGTGGTGTGGGAAGGTTTACTCAgcgcggcgcg                 >  2:7779158/1‑150 (MQ=255)

                                                                                                                                        cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCt               >  1:7779568/1‑150 (MQ=255)

                                                                                                                                          ctctctCGCTGACTCTCGTTCAGATGCTTAGTCGACGGTTCTCATCCTTCCCCGCTTGGGCAGATTATTTGATTGCGGATTCGCTTGAGAGTTCTTTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc             <  1:1743351/150‑1 (MQ=14)

                                                                                                                                           tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg            <  1:6470146/150‑1 (MQ=255)

                                                                                                                                           tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg            <  1:3262036/150‑1 (MQ=255)

                                                                                                                                             tctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc          <  1:1435909/150‑1 (MQ=255)

                                                                                                                                             tctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc          <  1:5554213/150‑1 (MQ=255)

                                                                                                                                             tctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAACCTCTGCGACCATCGcc          >  2:3306329/1‑150 (MQ=17)

                                                                                                                                                 gCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc      >  1:3373148/1‑150 (MQ=255)

                                                                                                                                                   tGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCgg    >  2:7428142/1‑150 (MQ=255)

CTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGT  >  NC_000913/1286816‑1287114

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