BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,498,612	C→A	11.2%	R213R (CGC→CGA)	insQ →	IS609 transposase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,498,612	0	C	A	11.2%	660.0 / 15.5	215	R213R (CGC→CGA)	insQ	IS609 transposase B
Reads supporting (aligned to +/- strand): ref base C (122/61); new base A (2/30); total (124/91)
Fisher's exact test for biased strand distribution p-value = 5.43e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

TGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCA  >  CP009273/1498465‑1498760

tGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCa                                                                                                                                                     >  1:4308745/1‑149 (MQ=255)

    gTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAAt                                                                                                                                                 >  1:1274631/1‑149 (MQ=255)

    gTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCATGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAAt                                                                                                                                                 <  2:871908/149‑1 (MQ=255)

      ggCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATAc                                                                                                                                               >  2:4308014/1‑149 (MQ=255)

      ggCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATAc                                                                                                                                               >  2:488352/1‑149 (MQ=255)

         tAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTTTAAACAGTTTTAATAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGc                                                                                                                                            <  1:418292/149‑1 (MQ=255)

         tAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGc                                                                                                                                            <  1:4396037/149‑1 (MQ=255)

           aaCTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTTTAAACAGTTTTCAGAAAAACCTGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGa                                                                                                                                          <  1:1249443/149‑1 (MQ=255)

           aaCTCGCCACGCTGTCAGATGGCACAGTATTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGa                                                                                                                                          <  1:1952445/149‑1 (MQ=255)

              tCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTg                                                                                                                                       >  2:3781750/1‑149 (MQ=255)

                   aCGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTc                                                                                                                                  >  1:1367817/1‑149 (MQ=255)

                   aCGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTc                                                                                                                                  >  1:1367880/1‑149 (MQ=255)

                     gCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCt                                                                                                                                >  1:210785/1‑149 (MQ=255)

                          cAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATc                                                                                                                           >  1:3378894/1‑149 (MQ=255)

                               ggCACTGTCTTTGCGCCTTTAAACTGTTTTAATAAAAACAAGAACACGCGGGCGATACTTCAGCGACAGTTAAGCCGAAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACCGCGACTGCATTCCTGTATCGCaaa                                                                                                                      <  1:1466934/149‑1 (MQ=255)

                               ggCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCaaa                                                                                                                      >  2:244210/1‑149 (MQ=255)

                               ggCACAGTCTTTGAGACTGTAAACAGTTTTAAGAAAAACCATAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGAAAGTACAAATTCAGCAACAACTGGAAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCaaa                                                                                                                      <  2:1765795/149‑1 (MQ=255)

                                   caGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAAAAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATc                                                                                                                  <  2:2255828/149‑1 (MQ=255)

                                       cTTTGAGCCTTTAAACTGTTTTCAGAAAAACCATAAGAAGCTGGCTAGACTTCAGCGACATTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCa                                                                                                              <  2:1817523/149‑1 (MQ=255)

                                           gAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGa                                                                                                          >  2:2382902/1‑149 (MQ=255)

                                           gAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGa                                                                                                          >  2:1671378/1‑149 (MQ=255)

                                           gAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGa                                                                                                          >  1:3789071/1‑149 (MQ=255)

                                             gCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACt                                                                                                        >  1:1051772/1‑149 (MQ=255)

                                              ccTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTa                                                                                                       >  2:1279338/1‑149 (MQ=255)

                                                tGTAAACTGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCATGGTCAACTTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTAcc                                                                                                     <  2:700060/149‑1 (MQ=255)

                                                   aaaCAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTc                                                                                                  >  2:644733/1‑149 (MQ=255)

                                                   aaaCAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTc                                                                                                  >  1:497347/1‑149 (MQ=255)

                                                     aCAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTcac                                                                                                >  1:95829/1‑149 (MQ=255)

                                                         ttttCATAAAAACCAGAAGACGCTGGCTAGAATTCAGCGACAGTTAAGCCGCAAGGTAAAATTCAGAAACAACTGGCAAAAGCAGAAACGCAAAATACCGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAg                                                                                            <  2:3290118/149‑1 (MQ=255)

                                                         ttttCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAg                                                                                            >  1:86766/1‑149 (MQ=255)

                                                         ttttCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAg                                                                                            >  1:4473646/1‑149 (MQ=255)

                                                         ttttCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAg                                                                                            >  1:1102378/1‑149 (MQ=255)

                                                           ttCAGAAAAACCAGAAGAAGCTGGCGAGACGTCAGCGACAGTTAAGCCGAAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTc                                                                                          <  2:1058148/149‑1 (MQ=255)

                                                             cAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTcac                                                                                        <  2:1472149/149‑1 (MQ=255)

                                                               taaaaaaCACAAGAAGCTGGCTAGACTTAAGCTAATTTAAAGCCGAAAGGTAAAATTCACAAACAACTGGAAAAAGCAGAAACGAAAAATACCGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACaa                                                                                      <  1:4518745/148‑1 (MQ=255)

                                                               taaaaaCCAGAATAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGAAAGGTAAAATTCAGCAACAACTGGCAAAAGCATAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACaa                                                                                      <  1:251933/148‑1 (MQ=255)

                                                                aaaaaCAATAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGAAAAAACTGGCAAAAGATGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAAc                                                                                     <  1:381950/149‑1 (MQ=255)

                                                                 aaaaCCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACg                                                                                    >  1:645488/1‑149 (MQ=255)

                                                                 aaaaCCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACg                                                                                    >  1:3717661/1‑149 (MQ=255)

                                                                  aaaCCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACCCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACgg                                                                                   >  2:3038666/1‑149 (MQ=255)

                                                                  aaaCCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGAAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACgg                                                                                   <  2:176295/149‑1 (MQ=255)

                                                                    aCCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGcc                                                                                 >  1:2816117/1‑149 (MQ=255)

                                                                     ccAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCg                                                                                >  2:1957350/1‑149 (MQ=255)

                                                                        gaatatgCTGGCGATACTTAAGCGACATTTACGCCGAAAGTTAAAATTCAGAAAAAACTGTAAAAAGCATAAACTAAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCa                                                                             <  2:901979/143‑1 (MQ=255)

                                                                        gaagaaGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCa                                                                             >  1:3302957/1‑149 (MQ=255)

                                                                          ataagATGGCGAGAATTCAGAGACAGTAAAGACGAAAGGTCAAATTCAGCAACAAATTGAAAAAGAATAAAAGAAAAAAAAAGCGACTGCATTACTGTATAGAAAATATCATAATGTAATACCTTCACAAAGTCACAACGACATTCAGc                                                                           >  2:3661128/3‑149 (MQ=255)

                                                                          agaagCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGc                                                                           >  1:1161642/1‑149 (MQ=255)

                                                                          agaagCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGc                                                                           >  2:2297988/1‑149 (MQ=255)

                                                                           gaagCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCa                                                                          >  2:3443902/1‑149 (MQ=255)

                                                                            aagCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCaa                                                                         >  1:3121025/1‑149 (MQ=255)

                                                                               cTGGCGAGAATTCAGCTAATTTTATGCCGAATGTTAAATTTCAGAAACACCTGCAAAAACCAGAAACCAAAAATACCGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCaaaaa                                                                      <  1:1697758/149‑1 (MQ=255)

                                                                                 ggCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAAcc                                                                    >  2:3214078/1‑149 (MQ=255)

                                                                                 ggCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAAcc                                                                    >  2:3597897/1‑149 (MQ=255)

                                                                                  gCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCa                                                                   >  2:2544343/1‑149 (MQ=255)

                                                                                  gCGAGACTTCAGCGACAGTTAAGCCGAATGTTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGAAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCa                                                                   <  1:3417246/149‑1 (MQ=255)

                                                                                     agaCTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGc                                                                >  2:984256/1‑149 (MQ=255)

                                                                                      gaCTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCa                                                               >  1:1131996/1‑149 (MQ=255)

                                                                                      gaCTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCa                                                               >  1:233544/1‑149 (MQ=255)

                                                                                       agttCAGCGACTGTTAAGCGGAATGTTCAATTTCAGCAACAACTGGAAAATGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCaa                                                              <  2:821370/147‑1 (MQ=255)

                                                                                          tCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGa                                                           >  1:4202367/1‑149 (MQ=255)

                                                                                                 cAGTTAAGCCGAAAGGTAAAATTAAGAAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCAt                                                    <  2:3688938/149‑1 (MQ=255)

                                                                                                    ttAAGCCTCATGGTAAAATTCCGAAAAAACTGTAAAAAGAATAAACGAAAAATCCCGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGa                                                 <  1:3215210/149‑1 (MQ=255)

                                                                                                    ttAAGCCGCAAGGTCAAATTCAGCAACAACGGGCAAAAGCAGAAACGCAAAAGACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGa                                                 <  2:220003/149‑1 (MQ=255)

                                                                                                      aaGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAgg                                               <  2:1367880/149‑1 (MQ=255)

                                                                                                      aaGCCGCAAGGTCAAATTCAGCAACAACTGGAAAACGCAGAAACGAAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAgg                                               <  2:1367817/149‑1 (MQ=255)

                                                                                                        gCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGAt                                             >  2:1353146/1‑149 (MQ=255)

                                                                                                        gCCGCAAGGACAAATTCAGAAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGAt                                             <  2:2664382/149‑1 (MQ=255)

                                                                                                         ccGAAAGGAAAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGAtt                                            <  2:940623/149‑1 (MQ=255)

                                                                                                          ctcAAGGTAAAATTCAGAAACAACTGGCAAAAGCAGAAACGAAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGAttt                                           <  2:965789/147‑1 (MQ=255)

                                                                                                          cGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACCGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGAttt                                           <  2:2002241/149‑1 (MQ=255)

                                                                                                           gCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTg                                          >  1:1663172/1‑149 (MQ=255)

                                                                                                           gCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTg                                          >  2:460351/1‑149 (MQ=255)

                                                                                                             atgttaaaaTTCAGCAACAACTGCAAAAAGCAAAAACGAAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGaa                                        <  1:4501927/143‑1 (MQ=255)

                                                                                                              aGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAg                                       >  2:272474/1‑149 (MQ=255)

                                                                                                               ggTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAgg                                      >  2:753596/1‑149 (MQ=255)

                                                                                                                  cAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCa                                   <  2:3052033/149‑1 (MQ=255)

                                                                                                                    aaTTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATCGTCATTGAGGATTTGAAGGTCAGc                                 >  2:4071831/1‑149 (MQ=255)

                                                                                                                    aaTTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGc                                 >  2:1227181/1‑149 (MQ=255)

                                                                                                                    aaTTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGc                                 >  2:2709032/1‑149 (MQ=255)

                                                                                                                    aaTTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGc                                 >  2:4095162/1‑149 (MQ=255)

                                                                                                                     aTTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCa                                >  1:3494328/1‑149 (MQ=255)

                                                                                                                          gCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCACTGATAGTCATTGAGGATTTGAAGGTCAGCAACATg                           >  2:3003134/1‑149 (MQ=255)

                                                                                                                            aacaacTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTc                         >  2:3847729/1‑149 (MQ=255)

                                                                                                                             acaacTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCa                        >  1:3740963/1‑149 (MQ=255)

                                                                                                                             acaacTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCCTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCa                        >  2:3504972/1‑149 (MQ=255)

                                                                                                                                 cTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGt                    >  2:4504163/1‑149 (MQ=255)

                                                                                                                                  tGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTc                   >  2:4144646/1‑149 (MQ=255)

                                                                                                                                  tGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTc                   >  2:4315684/1‑149 (MQ=255)

                                                                                                                                      aaaaGCATAAACCAAAAAACCAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTcagca               <  1:838579/149‑1 (MQ=255)

                                                                                                                                      aaaaGCAGAAACGAAAAATACCGCGACTGCATTCCTGTATCGCAAATAGCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTcagca               <  1:2858081/149‑1 (MQ=255)

                                                                                                                                      aaaaGCACAACCCAAAAAACCCGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTcagca               <  1:1279338/149‑1 (MQ=255)

                                                                                                                                      aaaaCCACAAACCAAAAAAACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTcagca               <  1:1295949/149‑1 (MQ=255)

                                                                                                                                       aaaGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATCGTCCTTGAGGATTTGAAGGTCAGCACCATGTCAAAGTcagcag              >  2:3934137/1‑149 (MQ=255)

                                                                                                                                        aaGCAGAAACGAAAAATACCGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCagcagc             <  1:4108673/149‑1 (MQ=255)

                                                                                                                                          gCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCgg           >  1:2748307/1‑149 (MQ=255)

                                                                                                                                          gCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACAGGTCAAAGTCAGCAGCgg           >  1:454412/1‑149 (MQ=255)

                                                                                                                                            aGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGt         >  2:3983320/1‑149 (MQ=255)

                                                                                                                                             gAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTa        >  2:1874316/1‑149 (MQ=255)

                                                                                                                                                   cAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTca  >  1:1218291/1‑149 (MQ=255)

                                                                                                                                                   cAAAAAACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTca  <  2:505827/149‑1 (MQ=255)

TGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGAAGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCA  >  CP009273/1498465‑1498760

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