BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	329,519	C→A	8.7%	Y19* (TAC→TAA)	betT →	choline:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	329,519	0	C	A	8.7%	376.9 / 16.7	114	Y19* (TAC→TAA)	betT	choline:H(+) symporter
Reads supporting (aligned to +/- strand): ref base C (54/50); new base A (4/6); total (58/56)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 121 total aligned reads displayed.

TAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGT  >  NC_000913/329370‑329668

tAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTac                                                                                                                                                       >  2:621536/1‑150 (MQ=255)

tAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTac                                                                                                                                                       >  2:38627/1‑150 (MQ=255)

tAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTac                                                                                                                                                       <  1:745639/150‑1 (MQ=255)

  aaGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTacac                                                                                                                                                     <  1:782228/150‑1 (MQ=255)

   aGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGTCAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACAcc                                                                                                                                                    <  2:781287/150‑1 (MQ=255)

   aGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACAcc                                                                                                                                                    >  2:98664/1‑150 (MQ=255)

   aGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACAcc                                                                                                                                                    >  2:1146824/1‑150 (MQ=255)

       gTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTAAAACACCAccg                                                                                                                                                <  1:840719/150‑1 (MQ=255)

            tgattgGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGa                                                                                                                                           >  1:290031/1‑150 (MQ=255)

                 ggACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGAt                                                                                                                                      >  2:177002/1‑150 (MQ=255)

                  gACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGAtt                                                                                                                                     <  1:1514871/150‑1 (MQ=255)

                   aCGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTccgccg                                                                                                                                             >  2:1467441/1‑141 (MQ=255)

                   aCGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTccgccg                                                                                                                                             <  1:1467441/141‑1 (MQ=255)

                   aCGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGAttt                                                                                                                                    <  1:1084310/150‑1 (MQ=255)

                     gTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttg                                                                                                                                  >  1:569108/1‑150 (MQ=255)

                      ttCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttgt                                                                                                                                 <  1:1363469/150‑1 (MQ=255)

                      ttCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttgt                                                                                                                                 <  2:928779/150‑1 (MQ=255)

                      ttCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttgt                                                                                                                                 <  1:1530886/150‑1 (MQ=255)

                       tCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttgtt                                                                                                                                <  1:422235/150‑1 (MQ=255)

                        cAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGTGGCT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttgttg                                                                                                                               <  1:1488680/150‑1 (MQ=255)

                        cAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGTCAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttgttg                                                                                                                               <  1:1000697/150‑1 (MQ=255)

                        cAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttgttg                                                                                                                               <  2:462239/150‑1 (MQ=255)

                         aaTATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTtgttgt                                                                                                                              >  2:381880/1‑150 (MQ=255)

                          atatAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTgttgtt                                                                                                                             >  1:306325/1‑150 (MQ=255)

                           tataAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGttt                                                                                                                            <  1:119049/150‑1 (MQ=255)

                           tataAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGttt                                                                                                                            <  2:102096/150‑1 (MQ=255)

                             taAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAAACCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTc                                                                                                                          >  2:948913/1‑150 (MQ=255)

                             taAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAACAACAAAATCAGTCCGGCGGAGGTGTAAAACACCACCGGATTGATTTTGTCCTTTTc                                                                                                                          <  1:762182/150‑1 (MQ=18)

                               aaaTGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTccc                                                                                                                        <  1:1353947/150‑1 (MQ=255)

                               aaaTGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAACAACAAAATCAGTCCGGCGGAGGTGTAAAACACCACCGGATTGATTTTGTCCTTTTccc                                                                                                                        <  1:786173/150‑1 (MQ=18)

                               aaaTGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAACAACAAAATCAGTCCGGCGGAGGTGTAAAACACCACCGGATTGATTTTGTCCTTTTccc                                                                                                                        <  1:176219/150‑1 (MQ=18)

                                aaTGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAACAACAAAATCAGTCCGGCGGAGGTGTAAAACACCACCGGATTGATTTTGTCCTTTTCCCt                                                                                                                       <  1:837311/150‑1 (MQ=18)

                                  tgtgTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGa                                                                                                                     <  1:1275589/150‑1 (MQ=255)

                                          aaTTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATc                                                                                                             <  2:392063/150‑1 (MQ=255)

                                           aTTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATcc                                                                                                            <  1:826505/150‑1 (MQ=255)

                                            ttGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCt                                                                                                           >  1:282093/1‑150 (MQ=255)

                                              gTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGt                                                                                                         >  1:575178/1‑150 (MQ=255)

                                                  cGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTcg                                                                                                     >  2:289121/1‑150 (MQ=255)

                                                   taaTTTTTTTTAAAATAGAAAAAATAACA‑TGGGGCTAACTGTATGCCATACCTTTCACACAGCCGGTAAACGTACAAATTAAAT‑CGGTGGGTGTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTcgcg                                                                                                   <  1:638187/149‑1 (MQ=255)

                                                   gAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttg                                                                                                                                  <  2:195903/120‑1 (MQ=255)

                                                   gAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTcgc                                                                                                    >  2:125289/1‑150 (MQ=255)

                                                     aTTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGa                                                                                                  >  2:1149250/1‑150 (MQ=255)

                                                     aTTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATACTGTTTCGCGa                                                                                                  >  2:1361508/1‑150 (MQ=255)

                                                        tGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACtt                                                                                               >  1:354006/1‑150 (MQ=255)

                                                          aTTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCTGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTtct                                                                                             >  1:362691/1‑150 (MQ=255)

                                                                  tAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCtg                                                                                     <  2:1395458/150‑1 (MQ=255)

                                                                    gTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCtgtg                                                                                   <  1:188421/150‑1 (MQ=255)

                                                                         aaTAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTg                                                                              <  1:766373/150‑1 (MQ=255)

                                                                         aaTAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTg                                                                              >  1:1083741/1‑150 (MQ=255)

                                                                          aTAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTgg                                                                             <  2:1083741/150‑1 (MQ=255)

                                                                                  ggggAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGc                                                                     >  1:1381460/1‑150 (MQ=255)

                                                                                   gggAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCt                                                                    >  1:16407/1‑150 (MQ=255)

                                                                                          cTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTg                                                              <  2:290309/150‑1 (MQ=255)

                                                                                          cTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTg                                                              <  2:700083/150‑1 (MQ=255)

                                                                                            ggATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTggg                                                            >  1:1354692/1‑150 (MQ=255)

                                                                                             gATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGt                                                           <  1:540747/150‑1 (MQ=255)

                                                                                                gACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTc                                                        <  2:354006/150‑1 (MQ=255)

                                                                                                 aCAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCt                                                       <  2:306325/150‑1 (MQ=255)

                                                                                                  cAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTa                                                      <  1:571226/150‑1 (MQ=255)

                                                                                                    gACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTaaa                                                    >  2:1002096/1‑150 (MQ=255)

                                                                                                        tttCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCt                                                >  1:1032506/1‑150 (MQ=255)

                                                                                                        tttCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCt                                                >  1:1017170/1‑150 (MQ=255)

                                                                                                          tCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTc                                              <  2:980474/150‑1 (MQ=255)

                                                                                                          tCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTc                                              <  2:1355093/150‑1 (MQ=255)

                                                                                                             cacaGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTAGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGt                                           >  2:567389/1‑150 (MQ=255)

                                                                                                              ccagcagGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACCCCGCCGCCGGACTGATTTTGTTGTTTTCCCTGACCACGCTCCTGTGTCGCGCCTTCTCGGCCCGGTGGATTGGCCGCACGCTGGACTGGGGGTCTAAAACCTTCgggt                                          >  1:1401809/2‑148 (MQ=255)

                                                                                                              acaGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGtt                                          <  1:356209/150‑1 (MQ=255)

                                                                                                               tgtcagGGAAAACAACAAAATCAGTCCGGCGGAGGTGTAAAACACCAGCGGATTGATTTTGTCCTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTg                                         >  1:975569/4‑150 (MQ=17)

                                                                                                                  cagGGAAAACAACAAAATCAGTCCGGCGGAGGTGTAAAACACCACCGGATTGATTTTGTCCTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTa                                      >  1:36213/1‑150 (MQ=18)

                                                                                                                   agGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTAc                                     >  1:693161/1‑150 (MQ=255)

                                                                                                                    gggAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACt                                    <  2:501759/150‑1 (MQ=255)

                                                                                                                     ggAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTa                                   >  1:924570/1‑150 (MQ=255)

                                                                                                                      gAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTAt                                  >  1:1273645/1‑150 (MQ=255)

                                                                                                                         aaGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATctg                               >  1:645512/1‑150 (MQ=255)

                                                                                                                         aaGGACAAAATCAATCCGGTGGTGTTTTACAACTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGACGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATctg                               >  2:1369338/1‑150 (MQ=255)

                                                                                                                           ggACAAAATCAATCCCGTGGTGTTTTACACATCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATctgat                             >  2:849853/1‑148 (MQ=255)

                                                                                                                            gACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATctgctg                            >  1:281965/1‑150 (MQ=255)

                                                                                                                            aaCAAAATCAGTCCGGCGGAGGTGTAAAACACCACCGGATTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATctgctg                            >  2:52913/2‑150 (MQ=25)

                                                                                                                             aCAAAATCAGTCCGGCGGAGGTGTAAAACACCAGCGGATTGATTTTGTCCTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTgg                           <  2:975569/150‑1 (MQ=21)

                                                                                                                              cAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTggc                          >  2:1379398/1‑150 (MQ=255)

                                                                                                                              cAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTggc                          >  2:1335997/1‑150 (MQ=255)

                                                                                                                                aaaTCAGTCCGGCGGAGGTGTAAAACACCACCGGATTGATTTTGTCCTTTTCACTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTggcgg                        >  2:689260/1‑150 (MQ=25)

                                                                                                                                aaaTCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTggcgg                        >  2:1210063/1‑150 (MQ=255)

                                                                                                                                 aatcaatcCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTggcggc                       >  1:695484/1‑150 (MQ=255)

                                                                                                                                  atcaatcCGGTTGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCa                      <  2:819941/150‑1 (MQ=255)

                                                                                                                                  atcaatcCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCa                      >  1:160064/1‑150 (MQ=255)

                                                                                                                                  atcaatcCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGATTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGTCa                      <  2:825197/150‑1 (MQ=255)

                                                                                                                                   tcaatcCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCaa                     >  1:1331653/1‑150 (MQ=255)

                                                                                                                                     aatcCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACg                   >  2:172526/1‑150 (MQ=255)

                                                                                                                                      atcCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCATGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTCGCGGCAACGc                  >  2:295128/1‑150 (MQ=255)

                                                                                                                                        ccGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGctc                <  2:569108/150‑1 (MQ=255)

                                                                                                                                         cGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGctct               <  1:1149250/150‑1 (MQ=255)

                                                                                                                                         cGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGctct               >  2:1063133/1‑150 (MQ=255)

                                                                                                                                          ggtggtGTTTTACACCTCCGCCTGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCta              <  2:362691/150‑1 (MQ=255)

                                                                                                                                          ggtggtGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCTGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCta              <  1:125289/150‑1 (MQ=255)

                                                                                                                                          ggtggtGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCta              <  1:265680/150‑1 (MQ=255)

                                                                                                                                               tgtTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTg         >  1:3525/1‑150 (MQ=255)

                                                                                                                                                 ttttACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTc       >  1:891195/1‑150 (MQ=255)

                                                                                                                                                      cacCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTtgt  >  2:840473/1‑150 (MQ=255)

                                                                                                                                                      cacCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTtgt  <  2:59571/150‑1 (MQ=255)

TAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGT  >  NC_000913/329370‑329668

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