BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,927,597	G→T	9.6%	R415S (CGC→AGC)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,927,597	0	G	T	9.6%	250.8 / 4.1	101	R415S (CGC→AGC)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base G (38/51); new base T (2/10); total (40/61)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

ACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGG  >  NC_000913/1927448‑1927746

aCTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCg                                                                                                                                                       <  1:802964/150‑1 (MQ=255)

 cTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAACCCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCgg                                                                                                                                                      >  1:541647/1‑150 (MQ=255)

  tGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGt                                                                                                                                                     <  1:1643971/150‑1 (MQ=255)

       tCGGCATCAATACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTc                                                                                                                                                >  1:947595/1‑150 (MQ=255)

           cATCAATACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTcgccg                                                                                                                                            <  2:1461322/150‑1 (MQ=255)

            aTCAATACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACTACATCACGGGCGACTATCCAAAGGTGTTCACTGCGGTAATTcgccgc                                                                                                                                           >  2:1915241/1‑150 (MQ=255)

             tCAATACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCa                                                                                                                                          >  2:1093281/1‑150 (MQ=255)

              cAATACTTGCGCCGTAAGATCCTTAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCAt                                                                                                                                         <  1:107819/150‑1 (MQ=255)

                  aCTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATaaaa                                                                                                                                     <  2:953008/150‑1 (MQ=255)

                   cTTGCGCCTTAAGAACCATAGCCATACAACTGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCAATCACTGGCGACTATACACAGGTGTTAACTGATGTAAATAGCAGCATaaaaa                                                                                                                                    >  1:368952/1‑150 (MQ=255)

                   cTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATaaaaa                                                                                                                                    >  2:1510319/1‑150 (MQ=255)

                   cTTGCGCCGTAAGAACCATAGCCATACACAAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGGGTTCACTGCGGTAATTCGCCGCAaaaaaa                                                                                                                                    >  1:363430/1‑150 (MQ=255)

                    ttGCGCCGTAAGAACCATAGCCTTACACCCGCAACGGGTTGCGGCCTGTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGCACTTCGACGCCATCAGGGGCGACTATCCACAGGTGTTCACTGGGTTAATTCGCCGCATAAAAAc                                                                                                                                   <  2:1263567/150‑1 (MQ=255)

                    ttGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAAc                                                                                                                                   <  1:689817/150‑1 (MQ=255)

                     tGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAAcc                                                                                                                                  >  1:2008937/1‑150 (MQ=255)

                     tGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAAcc                                                                                                                                  <  1:844234/150‑1 (MQ=255)

                     tGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGCTAATTCGCCGCATAAAAAcc                                                                                                                                  >  2:1068850/1‑150 (MQ=255)

                      gcgcCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTCTCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCa                                                                                                                                 >  2:350723/1‑150 (MQ=255)

                       cgcAGTAAGACCCAAAGCCATACACCACGAACGGGTGGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAg                                                                                                                                <  1:170347/150‑1 (MQ=255)

                        gcCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAgg                                                                                                                               <  2:792407/150‑1 (MQ=255)

                                atccATAGCGATACACCAGCAACGGGTTGTGTCCTTTGCGAATATGTTTGCGATGGTAGACCAACGATACCGGTACTTCGACGCGATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCg                                                                                                                       <  1:883507/148‑1 (MQ=255)

                                  ccATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCgtt                                                                                                                     <  1:1466370/150‑1 (MQ=255)

                                  ccATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCgtt                                                                                                                     >  1:1286564/1‑150 (MQ=255)

                                   cATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCgttt                                                                                                                    <  1:2043928/150‑1 (MQ=255)

                                   cATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATAACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCgttt                                                                                                                    >  1:1173860/1‑150 (MQ=255)

                                       gCCATACACCAGCAACGGGTTGTGACCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCCCTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCgtttgttt                                                                                                                >  2:377707/1‑150 (MQ=255)

                                         cATACACCCGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCCTCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTa                                                                                                              <  2:858144/150‑1 (MQ=255)

                                         cATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTa                                                                                                              <  2:907653/150‑1 (MQ=255)

                                         cATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACGTCCGTTTGTTTTa                                                                                                              >  2:1550281/1‑150 (MQ=255)

                                          atccacCAGCAACTGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACGAACGAATCCGGAACTTCGAAGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTaa                                                                                                             <  2:1503231/147‑1 (MQ=255)

                                           tACTCCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAAt                                                                                                            <  2:1865648/150‑1 (MQ=255)

                                           tACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAAt                                                                                                            <  1:613981/150‑1 (MQ=255)

                                            acacCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATa                                                                                                           <  2:502511/150‑1 (MQ=255)

                                            acacCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATa                                                                                                           >  1:1010483/1‑150 (MQ=255)

                                            acacCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATa                                                                                                           <  2:680164/150‑1 (MQ=255)

                                             cacCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATac                                                                                                          >  2:1169800/1‑150 (MQ=255)

                                               ccAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATacac                                                                                                        >  1:1221730/1‑150 (MQ=255)

                                                  gCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGCAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACacgac                                                                                                     >  2:1805733/1‑150 (MQ=255)

                                                  gCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGCTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACacgac                                                                                                     >  1:1018240/1‑150 (MQ=255)

                                                    aaCGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCAGTTTGTTTTAATACACGACGc                                                                                                   >  1:78239/1‑150 (MQ=255)

                                                         gTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACATGTGTTAAATTCTGTAATTGGCCTCATATAAACAATGAACTTCAGTTTGTTTTAATACACGACGCTCAcc                                                                                              <  2:2003066/150‑1 (MQ=255)

                                                           tgggtCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCAAATTTTTTTAATGCTTTAATTGTCCTCATATATACAATGATCTTTCGTTTGTTTTAATACACGACGCTCACCgg                                                                                            <  2:460505/147‑1 (MQ=255)

                                                                ccTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGtatcc                                                                                       <  1:1608016/150‑1 (MQ=255)

                                                                   ttGCGAAAATGTTTGCTATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACATGTTTTAAATGCGTTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGtatccata                                                                                    <  2:487298/150‑1 (MQ=255)

                                                                     gCGAAACTGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGAGACTATCCTTATTTTTTTTTTTGTGTAATTGTCCTCATAAATATAAGGAACTTTAGTTTGTTTTAATACACGACGCTCACCGGtatccatatc                                                                                  <  2:1446766/150‑1 (MQ=255)

                                                                      cGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGtatccatatcc                                                                                 >  1:743235/1‑150 (MQ=255)

                                                                      cGAAAATCTTTGCGATGGTACACCAACGAAACCGGAACTTCGACGCCATCACGGGCGAGTATCCATATTTTTTTTATTGTTTATTTGTCCTCATATATATAATGAACTTTATATTGTTTAAATACACGACGCTCACCGGtatccatatcc                                                                                 <  2:269347/150‑1 (MQ=255)

                                                                           aTGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTc                                                                            <  2:69596/150‑1 (MQ=255)

                                                                            tGTTTGCGATGGTAGACCAACCAAACCGGAACATCGACGCCATCACGGGCGACTATCCAAATTTTTTTTATTGTTTATTTGTCCTCATATATATAATGATCTTTAGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTca                                                                           <  2:1191538/150‑1 (MQ=255)

                                                                               ttGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTcaaac                                                                        <  1:1698255/150‑1 (MQ=255)

                                                                               ttGCGATGGGAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCCCCGGTATCCATATCCAGTTcaaac                                                                        <  1:1525842/150‑1 (MQ=255)

                                                                                  cGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTcaaacaaa                                                                     <  1:499702/150‑1 (MQ=255)

                                                                                   gATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACTACGCTCACCGGTATCCATATCCAGTTCAAACAAAg                                                                    >  2:544587/1‑150 (MQ=255)

                                                                                   gATGGAACACCACCGACACCGGAACTTCGACGCCATCACGGGCGACTATCCAAATTTTTTTTATTGTTTAATAGTCCTCATATATATAATGATCTAAAGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAg                                                                    <  2:1176071/150‑1 (MQ=255)

                                                                                     tGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAgtg                                                                  >  1:1612041/1‑150 (MQ=255)

                                                                                        tAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTTTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCt                                                               <  2:1999942/150‑1 (MQ=255)

                                                                                          gACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCATTTTTTTTTTATTGTTTATGTGTCCTCATATATATAATGATCTTTAGATTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTgg                                                             <  2:145980/150‑1 (MQ=255)

                                                                                           aCCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGgt                                                            <  1:141931/150‑1 (MQ=255)

                                                                                           aCCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGgt                                                            >  2:1677845/1‑150 (MQ=255)

                                                                                            ccAACGAAACCTGAACTTCGACGCCATCACGGGCGACTATCCAAATTTGTTCACTGCTGTAATTCGCCGCATAAAAACCAGGAACTTACGTTTTTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGgtg                                                           <  2:911606/150‑1 (MQ=255)

                                                                                            ccAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGgtg                                                           <  1:1095433/150‑1 (MQ=255)

                                                                                            ccAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGgtg                                                           >  2:308476/1‑150 (MQ=255)

                                                                                             cAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGgtgt                                                          >  2:521390/1‑150 (MQ=255)

                                                                                             cAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCAATGCGTTATTTGGCCGCATAAATACCATGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGgtgt                                                          <  2:639117/150‑1 (MQ=255)

                                                                                             aaaCGAAACCGGAACTTCGACGCCATCACGGGCGAATATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGgtgt                                                          <  1:1726909/149‑1 (MQ=255)

                                                                                                cGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGt                                                       >  1:1701050/1‑150 (MQ=255)

                                                                                                cGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGt                                                       >  1:886405/1‑150 (MQ=255)

                                                                                                     ccGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTAACTTCGGTAATTCGCCGCATAAAAACCAGGAACTTTCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATgg                                                  <  2:1009027/150‑1 (MQ=255)

                                                                                                       ggAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGgaa                                                <  2:960691/150‑1 (MQ=255)

                                                                                                        gcaCTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCTGTAATTCTCCGCATAAAAACCAGGATCTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGgaag                                               <  2:1461501/148‑1 (MQ=255)

                                                                                                        gAACTTCGACGCCATCACGGGCGACTATCCACAGGTTTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGgaag                                               <  2:1707451/150‑1 (MQ=255)

                                                                                                            ttCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATa                                           <  1:1416647/150‑1 (MQ=255)

                                                                                                            ttCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATa                                           <  1:1653466/150‑1 (MQ=255)

                                                                                                             tCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGCCGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATaa                                          >  1:560618/1‑150 (MQ=255)

                                                                                                              cGACGCCATCACGGGCGACTATCTAGAGGTGTTCACTGCGGTACTTCGCCGCATAAATATCATGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAAc                                         <  2:113011/150‑1 (MQ=255)

                                                                                                                 cGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCaa                                      >  2:808028/1‑149 (MQ=255)

                                                                                                                    cATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACCCGACGCTCACCGGTATCCATATCCAGTTCAATCAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAAc                                   >  1:1573776/1‑150 (MQ=255)

                                                                                                                          gggCGACTATCCACAGGTGTTCACTGCGGTAATTCTCCGCATAAAAACCAGGATCTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATc                             <  2:1075946/150‑1 (MQ=255)

                                                                                                                           ggCGACTATCCACATGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTAACCGGTATCCATATCCAGATCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATcg                            <  1:377405/150‑1 (MQ=255)

                                                                                                                           ggCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATcg                            >  1:590644/1‑150 (MQ=255)

                                                                                                                               aCTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCAAATCCAGATCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATcgcgcg                        <  1:1241507/150‑1 (MQ=255)

                                                                                                                                cTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCgg                       <  1:1510320/150‑1 (MQ=255)

                                                                                                                                    ccACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGttta                   >  1:462917/1‑149 (MQ=255)

                                                                                                                                    ccACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTc                   >  1:1085810/1‑150 (MQ=255)

                                                                                                                                    ccACAGGTGTTCACTGCGGTAATTCGCAGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTc                   >  2:341483/1‑150 (MQ=255)

                                                                                                                                       caGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAgg                <  1:258766/150‑1 (MQ=255)

                                                                                                                                       caGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAgg                <  1:1770031/150‑1 (MQ=255)

                                                                                                                                        aGGAGTTCACTGCGGTAAATCGCCGCATATAAACCAGGTACTTCCGTTTGTTTTAATACACGACGCTCACCTGTTTCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGt               <  1:377707/150‑1 (MQ=255)

                                                                                                                                         ggTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGGCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGtt              >  1:349076/1‑150 (MQ=255)

                                                                                                                                            tttCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTg           <  2:680693/149‑1 (MQ=255)

                                                                                                                                            gtTCACTGCTTTAATTCTCCTCATAAATACCAGGATCTTTAGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTg           <  2:1173859/150‑1 (MQ=255)

                                                                                                                                              tCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGa         <  2:2008936/150‑1 (MQ=255)

                                                                                                                                                  tGCGGTAATTCTCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGt     <  2:572386/150‑1 (MQ=255)

                                                                                                                                                   gCGGTAATTCGCCGCATAAAAACCAGGATCTTCAGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTa    <  2:875756/150‑1 (MQ=255)

                                                                                                                                                   gCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTa    >  1:1745619/1‑150 (MQ=255)

                                                                                                                                                    cGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAg   >  1:1994306/1‑150 (MQ=255)

                                                                                                                                                    cGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAg   <  2:924817/150‑1 (MQ=255)

                                                                                                                                                    cGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGGATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAg   >  1:1084048/1‑150 (MQ=255)

                                                                                                                                                     ggTAATTCGCCTCATAAAAACCATGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAgg  <  2:557318/150‑1 (MQ=255)

                                                                                                                                                     ggTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAgg  <  1:1740869/150‑1 (MQ=255)

ACTGAAATCGGCATCAATACTTGCGCCGTAAGAACCATAGCCATACACCAGCAACGGGTTGTGTCCTTTGCGAAAATGTTTGCGATGGTAGACCAACGAAACCGGAACTTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGG  >  NC_000913/1927448‑1927746

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