BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,286,965	A→T	9.1%	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	9.1%	862.2 / 7.3	338	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 (129/156); new base T (8/26); total (148/190)
Fisher's exact test for biased strand distribution p-value = 1.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 340 total aligned reads displayed.

CTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG  >  NC_000913/1286816‑1287112

ctcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGa                                                                                                                                                     <  1:8970832/150‑1 (MQ=255)

   aTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGACGc                                                                                                                                                  >  1:3353061/1‑150 (MQ=255)

     ccTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTc                                                                                                                                                >  1:8602860/1‑150 (MQ=255)

     ccTTCCCCGCTGGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTGCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTc                                                                                                                                                <  1:975600/150‑1 (MQ=255)

       ttCCCCGCTTGGGCAGAATTCTTGATTGCGGATTCGTTTGAGCATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                              <  2:1705704/150‑1 (MQ=255)

       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                              >  1:6519383/1‑150 (MQ=255)

       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCACTCACGCTGACGCTCGa                                                                                                                                              >  1:10445843/1‑150 (MQ=255)

        tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                                                                                                             <  2:9036993/150‑1 (MQ=255)

         ccacGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                            <  2:552237/147‑1 (MQ=255)

          cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                                                           <  1:1699239/150‑1 (MQ=255)

           ccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTTCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                                                          >  1:5971348/1‑150 (MQ=255)

               ttGGGCAGATTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAcc                                                                                                                                      <  1:8230888/150‑1 (MQ=255)

                   gCAGAATATTTGATTGCGGATTCGTTTGAGAATTACGGGGGTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTTCGAGCTTCGCCCTTCGGGTCGTTGGCTGCGGCAACGCTCTCTCGCTGGCGCTCGTGGGGTACCTTAg                                                                                                                                  >  1:1290862/1‑150 (MQ=11)

                       actctttGATTGCGGATTCGTCTGAGCATTCCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                              <  2:6510430/146‑1 (MQ=255)

                       aaTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGACGAACCATTGTCGa                                                                                                                              >  2:417974/1‑150 (MQ=255)

                        atctttGATTGCGGAGTCGTTTGAGCATTCCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGaa                                                                                                                             <  2:10925032/147‑1 (MQ=255)

                         tagttGATTGCGGAGTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                            <  2:9247034/147‑1 (MQ=255)

                           tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAAACCACACAAGCTGACGCTCGAGACGAACCATAGTCGAAGCt                                                                                                                          >  2:4074481/1‑150 (MQ=14)

                            ttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGAGCTTCGCCCTTCGGGTCGATGCCGGCCGCAAAGCACTCTCGCTGACCCTAGAGCCGAACCTTAGTCGAAGCtt                                                                                                                         >  2:10816175/1‑150 (MQ=14)

                            ttGATTGCGGATTCGATTGAGAATTAAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTATCGCTGACGCTCGAGACGCAACGTAGTAGAAGcgt                                                                                                                         >  1:437897/1‑148 (MQ=11)

                              gATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtct                                                                                                                       <  2:5662559/150‑1 (MQ=255)

                                 tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGACGAACCTTGGTCGCAGCATCTCaa                                                                                                                    >  2:7010657/1‑149 (MQ=11)

                                     gATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCtt                                                                                                                <  1:7487391/150‑1 (MQ=255)

                                         cGTTTGAGAATTCCGTGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                            <  1:6384732/150‑1 (MQ=255)

                                            ttGAGTATTCGGTGGCTTTAGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                                         <  2:7385352/150‑1 (MQ=14)

                                             tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCTCTATCGCTGACGCGCGAGCCGAAACTTTGTCGAAGCTTATCATACGTACCCGCag                                                                                                        >  1:6239143/1‑149 (MQ=14)

                                                gtaTTCCGGGTCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                     <  1:10558447/148‑1 (MQ=255)

                                                gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACACTCGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCTAATACATCCCCGCCTggg                                                                                                     >  1:1574770/1‑150 (MQ=14)

                                                 aaTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGATGCGCGCTTAGGCCTTCGGGTCGTTTGCTGGGTCAAAGCTAAAAGGGTGTCGCTCGAGAAGAAACTTGGTCGAAGCTTGTCATCCTTCCCCGCATGGGc                                                                                                    >  2:6713531/1‑150 (MQ=2)

                                                         ggCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                            <  2:10614923/150‑1 (MQ=38)

                                                         ggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                            <  1:2087480/150‑1 (MQ=38)

                                                         ggCTTATGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTATCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                            <  2:695782/150‑1 (MQ=255)

                                                          gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCGGCCGCAACGCGCACTAGCGGACCCGCCAGCCCGACCATAGGAGAAGCTTCTCAACCAACACCGCAAGGGCACaagaatt                                                                                           >  1:3199759/1‑145 (MQ=2)

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

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

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

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

                                                            ttttGAAAGTGATGGTGGTGGGGGAAGGATTACAAAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGCCAACGCTCACAAGCTTACGCTCCACCCGAAACATAGTCCAAGCGTCACATCCATCCCCCCAAGGGCAGGATATTTgt                                                                                         >  2:2359325/1‑149 (MQ=2)

                                                              ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCAAAGTCGAATCTTCACAACCTTCCCCGCATGGGCAGAATATTTTAtt                                                                                       >  2:4278374/1‑150 (MQ=11)

                                                                  aaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGATGCCGGACGCAACGCTCTCTCGCTGACGCAACAGGCGAACCTTAGTCGAAGCTTCTAATACGTCCCCGCATGGGCAGAAGATATTATTGCgg                                                                                   >  2:1150510/1‑150 (MQ=11)

                                                                      gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                                               <  1:3154771/150‑1 (MQ=39)

                                                                        tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGc                                                                             <  2:7813578/150‑1 (MQ=40)

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

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

                                                                          gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCACGCTGACGCTCGAGCCCAACCTTAGTCGAAGCTTCTCATCCATCCCCGCAGGGGCAGAAAATTTGATTGCGGATTCGCtt                                                                           >  1:10004996/1‑150 (MQ=14)

                                                                           tggtggGGGTAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTg                                                                          <  1:4150845/150‑1 (MQ=255)

                                                                              tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                                       <  2:10825558/150‑1 (MQ=255)

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

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

                                                                                       gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGa                                                              <  1:3451710/150‑1 (MQ=255)

                                                                                              cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGa                                                       <  1:3496141/150‑1 (MQ=255)

                                                                                               aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCGGGGGAGTTCAGGGACTTTTGaa                                                      >  2:10095309/1‑150 (MQ=255)

                                                                                                gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                     <  1:8891456/150‑1 (MQ=255)

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

                                                                                                     gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAt                                                <  2:1620494/150‑1 (MQ=255)

                                                                                                      cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAg                                                                                                  >  2:6359192/1‑99 (MQ=18)

                                                                                                       gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCGTGGGAGTTCAGGGACTTTTGAAAGTGAtgg                                              >  2:6272866/1‑150 (MQ=255)

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

                                                                                                          cTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGGAGCTTCTCATCCTTCCCCGCAGGGGCCGAATATGTGATTTCGGATTCGCTGGGGCGTTCCGGGGCTGTTGTTAGTGGtggggg                                           >  2:10433791/1‑150 (MQ=14)

                                                                                                           ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                          <  1:11029956/150‑1 (MQ=255)

                                                                                                           ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                          <  2:8862155/150‑1 (MQ=255)

                                                                                                              gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                     <  1:3617544/88‑1 (MQ=255)

                                                                                                              gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                     >  2:3617544/1‑88 (MQ=255)

                                                                                                              gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                       <  2:11150018/150‑1 (MQ=255)

                                                                                                               cccTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGCACCTTAGTCGAAGCTTCTCATACTTCCCCGCATGGGCAGAATATTTTATTGCCGATTCGCTTTTGGGTTGAGGGGGGTTTGGAGGTGTTTGTGGTgggg                                      >  2:7516913/1‑150 (MQ=14)

                                                                                                                  ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                   <  2:4900181/150‑1 (MQ=255)

                                                                                                                  ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                   <  1:6285073/150‑1 (MQ=255)

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

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

                                                                                                                     agagcGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                >  1:9092579/5‑150 (MQ=255)

                                                                                                                       agcGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGGAGGAtt                              >  2:4388789/3‑150 (MQ=255)

                                                                                                                        tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGGGAGTTCAGGGACTGTTGAAAGTGATGGTGGTGGGGGAAGGATTa                             >  2:11083961/1‑150 (MQ=255)

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

                                                                                                                        tCGTGGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGTATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                             <  1:10709330/150‑1 (MQ=255)

                                                                                                                        gcggtggcTGCGGCATGGCTCTCTCGCTGTCGCTCGAGCGGACGGTTGGTCGTAGCGTCTCATCCTTCCCCGGTTGGGCAGATGTTTTGATTGCGGAGTCGCTTGAGTGTTCTGGGACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTa                             <  1:6171828/143‑1 (MQ=2)

                                                                                                                        gcGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGTTGGTGGTGGGGGAAGGATTa                             >  2:8128089/2‑150 (MQ=255)

                                                                                                                         cGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAc                            <  1:7328449/150‑1 (MQ=255)

                                                                                                                          gTTGACTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                           <  2:8054790/150‑1 (MQ=255)

                                                                                                                           ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGGTTGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTc                          >  2:2673655/1‑150 (MQ=255)

                                                                                                                           ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGTGAGTTCAGGGACTTTTGATAGTGATGGTGGTGGGGGTAGGATTACta                          >  2:3413288/1‑149 (MQ=255)

                                                                                                                           ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGTGAGTTCAGGGACTTTTGAAATTGATGGTGGTGGGGGACGGATTACTc                          >  1:6745120/1‑150 (MQ=255)

                                                                                                                            ggCCCGCGGCCCCGCTCTCTCGCGGTCGCGCGAGCCGCCCCTTCGTCGCAGCCTCCCAGCCGCCCCCGCCTGGGCAGCCTCTGCGCTTGCGGCGGCGCTTGAGGGTTCCGGGACTTTTGCAAGCGATGGTGGTGGGGGAAGGATTACTcc                         <  2:3763178/149‑2 (MQ=9)

                                                                                                                             gCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAg                        <  1:7774322/150‑1 (MQ=255)

                                                                                                                             gCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAg                        <  2:3122611/150‑1 (MQ=255)

                                                                                                                              ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTACCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgc                       >  1:6084590/1‑150 (MQ=255)

                                                                                                                                tGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTGAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGATTTGAGGGACGTTTGTAAGTGATGGGGGTGGGGGAAGGATTACTCAtcgc                     >  1:2568787/1‑150 (MQ=17)

                                                                                                                                 gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTACCCCGCATGGGCAGAATATGTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGACAGTGAGGGGGGTGGGGGAAGGAGTACTCAGCGCt                    >  2:9558221/1‑150 (MQ=17)

                                                                                                                                  cGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTTGTCGAAGCTTCTCTTCCTTCGCCGCATGGGGAGAATATTTGATTGCGGATTCGTTTGAGTGTTCTTGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTg                   <  1:3811547/150‑1 (MQ=11)

                                                                                                                                    gCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAATCTTCTCCTCCTTCCCCGCTTGGTCAGACTTATTGATTGCGTATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcg                 <  1:8077080/150‑1 (MQ=14)

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

                                                                                                                                        cGCTCTCTCGCCGGCGCTCGCGCCGAACCGTTGTCGAAGCGTCTCATCCTTCCCCGCTTGGGCAGATTATTTGATTGCGGATTCGATTGAGATTTCATGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCt             <  1:2595944/150‑1 (MQ=11)

                                                                                                                                         gCTCTCTCTCTGCAGCTCGAGTCGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCTTGGGCAGATTCTTTTATTGCGGATGAGCTTGAGTGTTCAGGGAATTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCtt            <  1:6579155/150‑1 (MQ=14)

                                                                                                                                         gCTCTCTCGCTGTCGCTCTAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCtt            <  1:10842885/150‑1 (MQ=255)

                                                                                                                                          ctctctCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc           <  1:692866/150‑1 (MQ=255)

                                                                                                                                           tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtgg                                              <  2:9865824/114‑1 (MQ=255)

                                                                                                                                            ctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGc         <  1:1848545/150‑1 (MQ=255)

                                                                                                                                            ctctcGCTGACGCGCTAGCCGACCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGc         <  1:4892766/150‑1 (MQ=255)

                                                                                                                                             tcgcGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc        <  2:5862986/147‑1 (MQ=255)

                                                                                                                                               tcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTTCGGATTAGCTTGTGAGTTCAGGGAATTTTGTAAGTTTTGGTGGTGGGGGAAGGATTACTAATCGCTGCGCGCTTCGCCCt      >  1:10139390/1‑150 (MQ=17)

                                                                                                                                                cGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCtt     <  2:8691733/150‑1 (MQ=255)

CTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGG  >  NC_000913/1286816‑1287112

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