BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	329,493	G→A	6.7%	D11N (GAC→AAC)	betT →	choline:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	329,493	0	G	A	6.7%	332.7 / 8.7	104	D11N (GAC→AAC)	betT	choline:H(+) symporter
Reads supporting (aligned to +/- strand): ref base G (49/47); new base A (2/5); total (51/53)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

Only 100 of 108 total aligned reads displayed.

TCCATTCATCAGCGGTGTTTATCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTG  >  NC_000913/329344‑329642

tCCATTCATCAGCGGTGTTTATCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAgg                                                                                                                                                       <  2:1537741/150‑1 (MQ=255)

tCCATTCATCAGCGGTGTTTATCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAgg                                                                                                                                                       <  2:788471/150‑1 (MQ=255)

  cATTCATCAGCGGTGTTTATCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGTTTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGAc                                                                                                                                                     >  1:501759/1‑150 (MQ=255)

        tcaGCGGTGTTTATCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAaatc                                                                                                                                               <  1:997619/150‑1 (MQ=255)

        tcaGCGGTGTTTATCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAaatc                                                                                                                                               <  2:1198300/150‑1 (MQ=255)

            cGGTGTTTATCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAaatcaatc                                                                                                                                           <  2:912185/150‑1 (MQ=255)

                    aTCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGgtggtg                                                                                                                                   >  2:826505/1‑150 (MQ=255)

                     tCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGgtgt                                                                                                                                  >  1:1301801/1‑150 (MQ=255)

                       tATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACaa                                                                                                                                                   <  1:73755/131‑1 (MQ=255)

                       tATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGttt                                                                                                                                >  1:1478397/1‑150 (MQ=255)

                        aTTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGtttt                                                                                                                               >  2:974933/1‑150 (MQ=255)

                        aTTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGtttt                                                                                                                               >  2:1135043/1‑150 (MQ=255)

                          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                                                                                                                           >  2:571226/1‑150 (MQ=255)

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

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

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

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

                             aGCGGTTATTGATTGGACGTTCAATATAAAATGTGCCTTAATTGTTACGAATTTGATTTAAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACAcc                                                                                                                          <  1:801534/150‑1 (MQ=255)

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

                                   tATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTccgcc                                                                                                                    >  2:938031/1‑150 (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:928779/1‑150 (MQ=255)

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

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

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

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

                                                 tCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttgtt                                                                                                      <  1:331883/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                                                                                                  <  2:102096/150‑1 (MQ=255)

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

                                                      ataAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGtttt                                                                                                 >  2:284746/1‑150 (MQ=255)

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

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

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

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

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

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

                                                          aaTGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGTCCTTTCACACAGCAGGGAAAACAACAAAATCAGTCCGGCGGAGGTGTAAAACACCACCGGTTTGATTTTGTCCTTTTCCCt                                                                                             <  2:1267/150‑1 (MQ=18)

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

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

                                                             gtgtCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGAc                                                                                          <  1:735760/150‑1 (MQ=255)

                                                                   tAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGAt                                                                                    >  2:356209/1‑150 (MQ=255)

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

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

                                                                             gAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttg                                                                                                        >  1:195903/1‑120 (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)

                                                                                tttGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGAc                                                                       >  1:1532710/1‑150 (MQ=255)

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

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

                                                                                         aaaTAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGcc                                                              >  2:1096831/1‑150 (MQ=255)

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

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

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

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

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

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

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

                                                                                                                    cTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTg                                    <  2:290309/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:1017170/1‑150 (MQ=255)

                                                                                                                                  tttCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCt                      >  1:1032506/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:621536/150‑1 (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)

                                                                                                                                              gggAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACt          >  2:791036/1‑150 (MQ=255)

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

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

                                                                                                                                                 aaaaGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATc       <  1:1146824/150‑1 (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)

TCCATTCATCAGCGGTGTTTATCTATTAAAGCGGTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTG  >  NC_000913/329344‑329642

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