BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	5,537,065	0	T	G	27.6%	324.9 / 15.2	197	intergenic (+42/+61)	nadE CDS (NH(3)‑dependent NAD(+) synthetase)/PP_4870 CDS (Azurin)	–/–
Reads supporting (aligned to +/- strand): ref base T (59/83); new base G (1/53); total (60/137)
Fisher's exact test for biased strand distribution p-value = 2.18e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 208 total aligned reads displayed.

TGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCG  >  AE015451_phaM1/5536929‑5537201

tGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                             <  1:6884714/137‑1 (MQ=255)

tGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                             <  1:3292352/137‑1 (MQ=255)

tGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                             <  1:18118682/137‑1 (MQ=255)

tGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                             <  2:3200123/137‑1 (MQ=255)

tGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                             <  2:7513196/137‑1 (MQ=255)

tGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCGTGATAAGTGGCCTGCAGGAGCGGGTTTACCCGCGAATACGAGGGt                                                                                                                                             <  2:14258301/137‑1 (MQ=255)

    tGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAc                                                                                                                                         <  2:9398037/137‑1 (MQ=255)

     gggCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGTAAGAGCGGGTTTCCCCGCGAATACGAGGGGGGAcc                                                                                                                                        >  1:4441367/1‑137 (MQ=255)

         aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAcc                                                                                                                                    <  1:6943099/137‑1 (MQ=255)

              gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGc                                                                                                                               <  1:3562625/137‑1 (MQ=255)

               ttGACCAGGCAGCGTTCGACATCATCGTCGTCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCa                                                                                                                              <  2:5914985/137‑1 (MQ=255)

               ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCa                                                                                                                              <  1:3086951/137‑1 (MQ=255)

               ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCa                                                                                                                              <  1:10027303/137‑1 (MQ=255)

               ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCa                                                                                                                              <  1:3750375/137‑1 (MQ=255)

                tGACCAGGCAGCTTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAATTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCAt                                                                                                                             <  2:13794111/137‑1 (MQ=255)

                tGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCAt                                                                                                                             <  2:9712351/137‑1 (MQ=255)

                tGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCAt                                                                                                                             <  2:14372957/137‑1 (MQ=255)

                tGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCAt                                                                                                                             <  1:12032125/137‑1 (MQ=255)

                tGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCAt                                                                                                                             <  2:9410989/137‑1 (MQ=255)

                tGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCACGAACTGCTGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGATGGTGAACCCACCGACGCAt                                                                                                                             <  2:15432701/137‑1 (MQ=255)

                   ccAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcg                                                                                                                          <  1:9569427/137‑1 (MQ=255)

                     aGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgcg                                                                                                                        <  2:15185497/137‑1 (MQ=255)

                                  catcttCTTCTCCACTCACCGCAAAACCCAGCAAAAGCGCAAACTCCCTTTCGCTCCTGGATAATTGGCCTGCAGGAGCGGGTTCACCCGCAAAGACGGGGGTGAACCCACCGACGCATTCGCGGGGGACCCCGCTc                                                                                                           <  2:4413968/137‑1 (MQ=21)

                                  catcatCGTCACCACCTCCCCAAAGACCCAACACAATCGCGAACTGCTGTGCGCTTCTTGAGAAGTGGCCTGCGGGAGCGTGTTCACCTGGAAGTAGGGGGGGGAACCCACCGACGCATTCGTGGGGGAACCCGCTc                                                                                                           <  2:12955624/137‑1 (MQ=255)

                                               acttACCGAAAGACCCAACACAAGCGTGAATGGCCGTTCGCTCCGGGAGAGGTGGCCGGAGGGAGCGGGTTCACCCGCGAATACGGGGGGGAACCCACCGACGCCTTCGCGGGTGACCCCGCTCCTGCAGGGTTCGc                                                                                              <  2:16834714/134‑1 (MQ=255)

                                                cttACCACAAGACCCAGCACAAGCGCGAATTGCTGTTCGCTCCTTGATAGGTGGCCGGCAGGAGAGGGTTCCCCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGcc                                                                                             <  2:12656583/135‑1 (MQ=37)

                                                       cAAGACCCGGCACAAGCGCGAACTGCTGTTCGCTCCTGGATAAGTGGCCTGAAGGGGCGGGTTCACCCGCGAATACGGGGGGGAACCCACCGACGCGTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGcc                                                                                      <  2:8848660/137‑1 (MQ=37)

                                                         acaCCCAGCACAATCGCGAACTGCCTTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTAAACCCACCGACGCATTCGCGGGTGAGCCCGCTCCTGCAGGGTTCGCCGTCAGCCtt                                                                                    <  2:2665579/135‑1 (MQ=38)

                                                          gACCCAACAAAAGCGCGAACTGCCGTTTTCTCTTGGATAGGTGCCCTGCGGGAGCGGGTTAACCCGAGAATACGGGGGGGAACCAACCGACGCGTTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCttt                                                                                   <  2:1350451/137‑1 (MQ=255)

                                                             ccAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAATTGGCCTGAAGGAGCGGGTTCACCCGCGAATACGAGGGGGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGAt                                                                                <  2:9078717/137‑1 (MQ=255)

                                                             ccAGCACAAGCGCGAACTGCCGTTCGCTCCTTGAGAAGTGGCCTGCAGGAGCGGTTTCACCCGCGAATACGAGGGTAAACCCACCGACGCGTTCGTGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGAt                                                                                <  2:17727043/137‑1 (MQ=255)

                                                                  acaAGCGCGAACTGCCTTTCGCTCCTGGAGAGGTGTCCTGGAGGAGTGGTTTCACCCGCGAAGACGAGGGTAAACCCACCGACGCGTTTGTGGGGGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                           <  2:16690826/137‑1 (MQ=255)

                                                                   caatagcgAACGTCCGTTCGCTCCTTGATAAGTGGCCTGCAGGGGTGGGTTCACCCGCGAGTACGAGGGTGAACCCACTGACGCGTTCGGGGGGGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTg                                                                          <  2:7888648/132‑1 (MQ=255)

                                                                   caAGCGCAAATTGCTGTTCGCTCCTTGATAATTGGCCTGCAGGAGTGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTTTCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTg                                                                          <  2:8204658/137‑1 (MQ=255)

                                                                    aatcgcgATCTGCCTTGTTCTCTTTGGTAGGTGGTCTGCAGGAGTGGGTTCCCATACGAGGACGGGGTTGAACCCACAGCACCATTTGCGGGTGAACCCTCTCCTGCAGGGTTCGCCGTCAGCCTGTGATTACTTGa                                                                         <  2:11952931/134‑1 (MQ=255)

                                                                     agtgcgAACTGCCGTTCGCTCTTTGATAATTTGCCTGAAGGAGGGGTTTCACCCGCAAGTACGAGGGGGAACCCACCGACGCGTTCGTGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGAc                                                                        <  2:16937672/134‑1 (MQ=255)

                                                                          gAACTGCCTTTCGCCCCTGGATAAGTGTCCCGCAGAGGTGGGTTCACCCGAGAATATGGGGGGGACCCCACCGACGCGTCTGCGGGGGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGa                                                                   <  2:10293934/137‑1 (MQ=255)

                                                                            cttGCTGTTCGCTCTTTGAGAGGTGTCCTGCAGGGGTGGTTTCACCCGCGAATACGGGGGTGAACCCACCAACGCGTTTGTGGGGGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGAcc                                                                 <  2:407340/135‑1 (MQ=255)

                                                                                 tgTTCGCTCCTGGATGAGTGGCTGGCAGGAGTGGGTTCACCCGCGAAGATGAGGGTAAACCCACCGACGCATTTGCGGGGGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGcc                                                            <  2:7182945/136‑1 (MQ=255)

                                                                                   ttCGCTCTTGGAGAATTGGCCTGCAGGAGGGGGTTCACCGGAAAAGACGAGGGGGAACCCACCGACGCGTTCGCGGGGAAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCtt                                                          <  2:12560659/137‑1 (MQ=255)

                                                                                       ccccTTGATAAGTTGCCCGCGTAAGCGGGTTCACCCGCGACTACGAGGGTGAACCCAACGACGCATTCACGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑tcat                                                   <  2:17337502/135‑1 (MQ=255)

                                                                                             gATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACTGACGCGTTCGTGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCga                                             <  2:11853286/137‑1 (MQ=255)

                                                                                             gATAAGGGTCCTGCAGGAGCGGGTTCACCAGAGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCga                                             <  2:2399491/137‑1 (MQ=255)

                                                                                              aTAGGTGGCCTGGAGGAGTGGGTTCACCCGTGGATATGAGGGGAAACCCACCGACGCGTTTGTGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCgag                                            <  1:4085190/137‑1 (MQ=255)

                                                                                                aaGTTGCCTGCGGGAGCGGTTTCACCCGCGAATACGGGGGTGAACCCACCGACGTGTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCTTCAGCCTTTGATTACTTTACGATTACCCTGCCTT‑‑‑TCACCATCgatcc                                          <  2:2653262/137‑4 (MQ=25)

                                                                                                        tgagggAGCGGGTTCACCCGCGGATACGAGGTTGAACCCACCGACGCGTTTGTGGTGGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCg                                  <  1:8758431/133‑1 (MQ=255)

                                                                                                        tgaaGGAGAGGATCCACCGGAAAATAAGGGGGTAAACCCACAGCCGTTTTCGCGGGGGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCg                                  <  1:2863353/134‑1 (MQ=255)

                                                                                                        tacgggAGCGGGTCCACCCGCGAATATGAGGGTGAACCCACCGACGCATTCGTGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCg                                  <  2:15642535/133‑1 (MQ=255)

                                                                                                        tacAGGAGTGGTTTCACCCGCGAATATGAGGGTGAACCCACCGACGCCTTCGCGGGGAACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCg                                  <  1:8411465/135‑1 (MQ=255)

                                                                                                        tGCAGGGGCGGGTTCACCCGAGAATACGGGGGTGAACCAACCGACGCATTCGCGGGGGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCg                                  <  2:4457875/137‑1 (MQ=255)

                                                                                                        tGCAGGAGTGGTTTCACCCGCGGAGAGGGGGGGGAACCCACCGATGCGTTTGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCg                                  <  2:1095501/137‑1 (MQ=255)

                                                                                                        ggaaggaggGGTTCCACCCGAAAAAAAGAGGGGAACCCCACCGACGCGTTCGGGGGTAACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCg                                  <  1:1446009/134‑1 (MQ=255)

                                                                                                        ggCAGGGGTGGTTTCACCCGGGAATACGGGGGTGAACCAACCGACGCGTTTGTGGGGAACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCg                                  <  1:17593875/136‑1 (MQ=255)

                                                                                                        agaaGGAGCGGGTCCACCCGCGAATATGAGGGTAAACCCACCGACGCGTTCGTGGGGAAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCg                                  <  1:5832744/134‑1 (MQ=255)

                                                                                                           aGGAGCGGGTTCACCCGCGAATACGATGGGGAACCCACCGACGCGTTCTCGGGTGAACCCGCTCCTGCAGGGTTCGCCTTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCATCGAGATGTGGCCCg                                  <  2:18496431/137‑1 (MQ=37)

                                                                                                             gAGCGGGTCAACCAGAGAAACCGAGGGGGAACCCACCGATGCATTTGTGGTCGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAg                             <  2:15316188/137‑1 (MQ=255)

                                                                                                              ggtgggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAgg                            >  1:4631344/4‑137 (MQ=255)

                                                                                                              aGCGGGTTCACCCGCGAATACGAGGGGGAACCCACCGACGCATTCGCGGGTAAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAgg                            <  2:1336226/137‑1 (MQ=255)

                                                                                                               gtgggTTCACTAGAAAATAGGGGGGGGAACCCCCCGACGCGTTATTGGTTACACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATTATCGAGATGTGGCCCGGGAAGGt                           <  1:7943596/135‑1 (MQ=255)

                                                                                                                 gggTTCACCCGCGAAGAGGGGGGTGAACCCACCGACGCATTCGTGGGGGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGc                         <  2:3023486/137‑1 (MQ=255)

                                                                                                                 gggTTCACCCGAGAAGATGGGGGGGAACCAACCGACGCGTTTGTGGGGAAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGc                         <  1:13038105/137‑1 (MQ=255)

                                                                                                                  ggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGcc                        >  1:1376721/1‑136 (MQ=255)

                                                                                                                   tttCACCCGAGAAGACGGGGGTGAACCCACCGAAGAATTCGCGGGGAAACCCGCTCCTGCAGGGTTCGCCGTCAGTCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCag                       <  1:2287061/136‑1 (MQ=255)

                                                                                                                   ttccACCCGAGAAGAAGAGGGTGAACCCACAGACGCATTCGCGGTTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCag                       <  1:2394115/134‑1 (MQ=255)

                                                                                                                   gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTTCCTT‑‑‑TCCTCATCGAGATGTGGCCCGGGAAGGTGCag                       >  1:14707603/1‑137 (MQ=255)

                                                                                                                   gTTCACCCGCGAAGACGAGGGGGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCag                       <  2:18195367/137‑1 (MQ=255)

                                                                                                                   gTTCACCCGAAAATACGAGGGGAAACCCACCGACGCGTTCGCGGGTAAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCag                       <  1:13610011/137‑1 (MQ=39)

                                                                                                                      cACCCTAGAAGACGAGGGGAAACCCACCGACGCATTCGCGGGGGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCagaag                    <  2:2896379/137‑1 (MQ=255)

                                                                                                                      cACCCGCGAATATGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCagaag                    <  2:12015461/137‑1 (MQ=255)

                                                                                                                      cACCCGCGAATACGGGGGGGAACCCACCGACGTATTCGAGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCagaag                    <  2:11447675/137‑1 (MQ=255)

                                                                                                                      cACCCGCGAATACGAGGGTAAACCCACCGACGTATTCGTGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCagaag                    <  2:9209172/137‑1 (MQ=255)

                                                                                                                      cACCCGCGAAGATGGGGGGGAACCCACCGACGCGTTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCagaag                    <  2:5422582/137‑1 (MQ=255)

                                                                                                                      cACCCGAGAAGAGGAGGGGGAACCCACAGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCagaag                    <  1:15568599/137‑1 (MQ=255)

                                                                                                                      cACCCCCGAAGACGGGGTTGAACCCACCGACGAATTCGAGGGGGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTG‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCagaag                    <  2:5613075/137‑1 (MQ=255)

                                                                                                                      cACCCAAGAATACGAGGGGAACCCCACCAACTCCTTCGGGGGTAACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCagaag                    <  1:17701407/137‑1 (MQ=255)

                                                                                                                       aCCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCagaaga                   >  2:15871662/1‑137 (MQ=255)

                                                                                                                        cccgaaaaTACGAGGGTGAACCCACCGACGCATTCGTGGGGGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAgaagaa                  <  2:11155345/131‑1 (MQ=255)

                                                                                                                         ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCACGGGAAGGTGCAGaagaag                 >  2:6927153/1‑137 (MQ=255)

                                                                                                                         ccGCAACTACGAGGGTGAAACCACCGCCTCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTGTGATTACTTGACGATGACCTTACCTT‑‑‑CCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                 <  2:5159923/137‑1 (MQ=255)

                                                                                                                          cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCCTCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGc                >  1:10113953/1‑137 (MQ=255)

                                                                                                                          cgcgAATACGAGGGTGAACCCACCGAAGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGc                >  1:532827/1‑137 (MQ=255)

                                                                                                                          cgcgAATACAAGGGTGAACCCACCGACGCACTCGCGGGTGAACCCGCTCCTGAAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGTCTT‑‑‑TCCTCATCAAGATGTGGCACGGGAAGGTGCAGAAGAAGc                >  2:1862335/1‑137 (MQ=255)

                                                                                                                            cgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg              >  1:12788871/1‑137 (MQ=255)

                                                                                                                            cgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg              >  1:8327106/1‑137 (MQ=255)

                                                                                                                             gAATACGAGGGGGAACCCACCGAGGCATTCGCGGGGGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGt             <  2:10639679/137‑1 (MQ=255)

                                                                                                                               aTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAg           >  1:906407/1‑137 (MQ=255)

                                                                                                                               aTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAg           >  1:1216299/1‑137 (MQ=255)

                                                                                                                               aTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAg           >  2:8578301/1‑137 (MQ=255)

                                                                                                                               aTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAg           >  1:6649963/1‑137 (MQ=255)

                                                                                                                               aTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAg           >  1:8540007/1‑137 (MQ=255)

                                                                                                                                tACGAGGGTGAACCCACCGACGCATTCGCGGGTGCACCCGCTCCTGCAGAGTTCGCCGTCGGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCACCGAGACGTAGCCCGGGAAAGCGCAGAAGAAGCTGTAGt          >  2:16185446/1‑137 (MQ=255)

                                                                                                                                tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt          >  1:697638/1‑137 (MQ=255)

                                                                                                                                tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt          >  1:5010971/1‑137 (MQ=255)

                                                                                                                                tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt          >  1:17545315/1‑137 (MQ=255)

                                                                                                                                 aCGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc         >  1:16880177/1‑137 (MQ=255)

                                                                                                                                 aCGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCAGGGAAGGTGCAGAAGAAGCTGTAGTc         >  1:18386758/1‑137 (MQ=255)

                                                                                                                                 aCGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCAAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc         >  2:18336428/1‑137 (MQ=255)

                                                                                                                                     gggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca     >  1:14864243/1‑137 (MQ=255)

                                                                                                                                      ggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac    >  1:11073701/1‑137 (MQ=255)

                                                                                                                                        tGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCg  >  1:946864/1‑137 (MQ=255)

TGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTT‑‑‑TCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCG  >  AE015451_phaM1/5536929‑5537201

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 9 ≤` `ATCG/ATCG` `< 10 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`