BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	329,526	G→A	8.9%	A22T (GCC→ACC)	betT →	choline:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	329,526	0	G	A	8.9%	340.8 / 15.6	112	A22T (GCC→ACC)	betT	choline:H(+) symporter
Reads supporting (aligned to +/- strand): ref base G (53/49); new base A (4/6); total (57/55)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

GTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTG  >  NC_000913/329377‑329674

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

                tCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTttgtt                                                                                                                                      <  1:331883/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)

                     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‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAACAACAAAATCAGTCCGGCGGAGGTGTAAAACACCACCGGATTGATTTTGTCCTTTTccc                                                                                                                              <  1:786173/150‑1 (MQ=18)

                        aaaTGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAACAACAAAATCAGTCCGGCGGAGGTGTAAAACACCACCGGATTGATTTTGTCCTTTTccc                                                                                                                              <  1:176219/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)

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

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

                                    aTTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATcc                                                                                                                  <  1:826505/150‑1 (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                                                                                                                                        >  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)

                                                           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)

                                                                           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)

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

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

                                                                                                   tCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTc                                                    <  2:980474/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)

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

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

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

                                                                                                             gggAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACt                                          >  2:791036/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)

                                                                                                                aaaaGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATc                                       <  1:1146824/150‑1 (MQ=255)

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

                                                                                                                  aaGGACAAAATCAATCCGGTGGTGTTTTACAACTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGACGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATctg                                     >  2:1369338/1‑150 (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)

                                                                                                                                 ccGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGctc                      <  2:569108/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)

                                                                                                                                        tggtttACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTg               <  1:1400376/147‑1 (MQ=255)

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

                                                                                                                                          ttttACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTc             <  2:1381460/150‑1 (MQ=255)

                                                                                                                                             tACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCttt          >  1:1314212/1‑150 (MQ=255)

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

                                                                                                                                                acCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTtgtg       <  1:930424/150‑1 (MQ=255)

                                                                                                                                                 ccTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTgg      >  1:934178/1‑150 (MQ=255)

                                                                                                                                                  cTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGt     <  2:1017170/150‑1 (MQ=255)

                                                                                                                                                   tCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTc    >  2:1067133/1‑150 (MQ=255)

                                                                                                                                                     cgccgGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTg  <  2:1032506/150‑1 (MQ=255)

GTTATTGATTGGACGTTCAATATAAAATGTGTCTTAATTGTTACGAATTTGATTTTAAATAGTAACAATAACAGTGGGGAT‑ACTGGATGACAGACCTTTCACACAGCAGGGAAAAGGACAAAATCAATCCGGTGGTGTTTTACACCTCCGCCGGACTGATTTTGTTGTTTTCCCTGACAACGATCCTGTTTCGCGACTTCTCGGCCCTGTGGATTGGCCGCACGCTGGACTGGGTTTCTAAAACCTTCGGTTGGTACTATCTGCTGGCGGCAACGCTCTATATTGTCTTTGTGGTCTG  >  NC_000913/329377‑329674

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