BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,869,780	G→A	100%	intergenic (‑92/‑160)	cysD ← / → iap	sulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,869,780	0	G	A	97.2%	704.4 / 7.0	216	intergenic (‑92/‑160)	cysD/iap	sulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand): ref base G (5/1); new base A (124/86); total (129/87)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 216 total aligned reads displayed.

TGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTA  >  CP009273/2869636‑2869925

tgGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAagta                                                                                                                                               >  2:761315/1‑149 (MQ=255)

 ggATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAagtag                                                                                                                                              <  2:1254664/149‑1 (MQ=255)

 ggATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAagtag                                                                                                                                              <  1:993628/149‑1 (MQ=255)

  gATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAagtagt                                                                                                                                             >  2:282009/1‑149 (MQ=255)

    tGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTc                                                                                                                                           >  2:2838056/1‑149 (MQ=255)

         tCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTaa                                                                                                                                      <  1:3622246/149‑1 (MQ=255)

          ccGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAAc                                                                                                                                     <  1:622457/149‑1 (MQ=255)

               tCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAAt                                                                                                                                >  2:701373/1‑149 (MQ=255)

                 cAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATaa                                                                                                                              >  2:207526/1‑149 (MQ=255)

                    ttGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCATTAACGGAATAACGa                                                                                                                           >  2:1192487/1‑149 (MQ=255)

                    tagCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGa                                                                                                                           <  2:1428630/147‑1 (MQ=255)

                      gCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTACTTTGCAATACCGCTATTTTCTTTCCATCAGATGTTTCGACTATCTGGAGCGTAAGAGAACGAATGAAATTAACAATTAGAATGAGTAGTTCCTTAACGGAATaaacatt                                                                                                                         >  1:1410090/1‑144 (MQ=255)

                        cGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTg                                                                                                                       <  1:2123817/149‑1 (MQ=255)

                         gCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAAGACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTgg                                                                                                                      <  2:3697786/149‑1 (MQ=255)

                          cAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGc                                                                                                                     <  1:1033763/149‑1 (MQ=255)

                          cAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGc                                                                                                                     <  2:2499289/149‑1 (MQ=255)

                           aGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCa                                                                                                                    >  1:1400867/1‑149 (MQ=255)

                                   aaGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAat                                                                                                            >  1:2357324/1‑149 (MQ=255)

                                    aGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAata                                                                                                           <  2:3926592/149‑1 (MQ=255)

                                      tCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATc                                                                                                         >  2:3544713/1‑149 (MQ=255)

                                       cGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCa                                                                                                        >  2:3443564/1‑149 (MQ=255)

                                       cGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCa                                                                                                        <  2:536392/149‑1 (MQ=255)

                                        gTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaa                                                                                                       >  2:3255476/1‑149 (MQ=255)

                                           tttGATCCATAACCGTTCCTTTGCAATACCGCTATTTTGGGGCCGGCCGCGGGGGCGCCGCGCGGGCGCGGCAGAGCCCGCATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaaaa                                                                                                    >  2:1878089/1‑149 (MQ=255)

                                           tttGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaaaa                                                                                                    >  1:762881/1‑149 (MQ=255)

                                            ttGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAg                                                                                                   >  2:407534/1‑149 (MQ=255)

                                              gATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTg                                                                                                 <  1:3940708/149‑1 (MQ=255)

                                                tCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCt                                                                                               <  1:207526/149‑1 (MQ=255)

                                                   aTAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTaa                                                                                            <  1:2907825/149‑1 (MQ=255)

                                                    tAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAg                                                                                           >  2:4006977/1‑149 (MQ=255)

                                                     aaCCGTTCCTTTGCAATACCGCTATTGTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGAAGTTCCTTAACGGAATAACGATATGTCAAAGCTAATAAAAAAAAGTGCATaatg                                                                                          >  1:2539204/1‑147 (MQ=255)

                                                       ccGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCa                                                                                        >  1:2719872/1‑149 (MQ=255)

                                                         gTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAcc                                                                                      >  2:1425610/1‑149 (MQ=255)

                                                          ttCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCg                                                                                     >  1:817746/1‑149 (MQ=255)

                                                           tCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCgg                                                                                    >  2:2239388/1‑149 (MQ=255)

                                                            ccTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGa                                                                                   >  1:310427/1‑149 (MQ=255)

                                                            ccTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGa                                                                                   >  1:2180749/1‑149 (MQ=255)

                                                             cTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAt                                                                                  <  2:3216682/149‑1 (MQ=255)

                                                             cTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAt                                                                                  <  1:904235/149‑1 (MQ=255)

                                                             cTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAt                                                                                  <  1:1085564/149‑1 (MQ=255)

                                                              tttGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCANTTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAtt                                                                                 >  1:2082668/1‑149 (MQ=255)

                                                               ttGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAttt                                                                                <  1:1068675/149‑1 (MQ=255)

                                                                    aTACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGc                                                                           >  1:4024326/1‑149 (MQ=255)

                                                                    aTACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGc                                                                           <  2:3965813/149‑1 (MQ=255)

                                                                    aTACCGCTATTTTCGTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGc                                                                           <  2:3966144/149‑1 (MQ=255)

                                                                           tATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAgg                                                                    >  2:272195/1‑149 (MQ=255)

                                                                              tttCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAg                                                                 >  1:3293902/1‑149 (MQ=255)

                                                                               ttCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGa                                                                <  1:282009/149‑1 (MQ=255)

                                                                                 cTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAtt                                                              <  1:3520217/149‑1 (MQ=255)

                                                                                  ttGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAttt                                                             <  1:1940807/149‑1 (MQ=255)

                                                                                   tGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTg                                                            <  2:85632/149‑1 (MQ=255)

                                                                                       aTCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt                                                        <  1:3126026/149‑1 (MQ=255)

                                                                                       aTCAGATGTTTAGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt                                                        <  1:3126141/149‑1 (MQ=255)

                                                                                        tCAGATGTTTCGACTATAGGGAGCGAAAGAGAACGAATTAAATTACCAATTATAATAAGTAGTACATTAAAGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTTCTTAAGGCACCGGATTTCTGGCGATTATGAAGATTTGAAAtt                                                       >  1:1512966/1‑149 (MQ=255)

                                                                                          aGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGt                                                     >  2:3029042/1‑149 (MQ=255)

                                                                                          aGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACAGGATTTCGGGCGTTTATGAAGATTTGAAATTGt                                                     >  2:2622054/1‑149 (MQ=255)

                                                                                           gATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGtt                                                    <  2:277971/149‑1 (MQ=255)

                                                                                           gATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGtt                                                    <  2:3567650/149‑1 (MQ=255)

                                                                                           gATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGtt                                                    >  1:1209589/1‑149 (MQ=255)

                                                                                           gATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGtt                                                    >  1:1795153/1‑149 (MQ=255)

                                                                                                 tCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAgcg                                              >  2:1793761/1‑149 (MQ=255)

                                                                                                    aCTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAg                                           >  2:1336580/1‑149 (MQ=255)

                                                                                                    aCTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTGGAAATTGTTTTAGCGCAg                                           >  1:2693347/1‑149 (MQ=255)

                                                                                                      tataGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCg                                         <  1:735479/149‑1 (MQ=255)

                                                                                                        taGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGc                                       <  2:168121/149‑1 (MQ=255)

                                                                                                         aGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCa                                      >  1:1571917/1‑149 (MQ=255)

                                                                                                          gggAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGccg                                     >  2:3630661/1‑147 (MQ=255)

                                                                                                          gggAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAg                                     <  2:1620650/149‑1 (MQ=255)

                                                                                                          gggAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAg                                     >  2:1584248/1‑149 (MQ=255)

                                                                                                          gggAGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAg                                     <  2:762881/149‑1 (MQ=255)

                                                                                                             aGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGttt                                  >  1:1837666/1‑149 (MQ=255)

                                                                                                             aGCGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGttt                                  >  1:1837987/1‑149 (MQ=255)

                                                                                                               cGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCa                                <  2:871857/149‑1 (MQ=255)

                                                                                                               cGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCa                                <  1:2211195/149‑1 (MQ=255)

                                                                                                                 tAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATa                              >  2:1365594/1‑149 (MQ=255)

                                                                                                                 tAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATa                              <  1:2626019/149‑1 (MQ=255)

                                                                                                                  aaGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATAc                             >  1:3867868/1‑149 (MQ=255)

                                                                                                                   agagAACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATATCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACt                            >  1:2527827/1‑149 (MQ=255)

                                                                                                                    gagaACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTa                           <  2:2072969/149‑1 (MQ=255)

                                                                                                                     agaACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTAt                          <  1:1736361/149‑1 (MQ=255)

                                                                                                                      gaACGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATg                         >  1:1390552/1‑149 (MQ=255)

                                                                                                                      gaACGAATGAAATTACAAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATg                         <  2:2508565/149‑1 (MQ=255)

                                                                                                                       aaCGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATgg                        >  2:2670712/1‑149 (MQ=255)

                                                                                                                         cGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCg                      >  1:1002286/1‑149 (MQ=255)

                                                                                                                         cGAATGAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCg                      <  2:2106583/149‑1 (MQ=255)

                                                                                                                              gAAATTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaa                 >  2:499860/1‑149 (MQ=255)

                                                                                                                                aaTTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaaaa               >  2:521871/1‑149 (MQ=255)

                                                                                                                                 aTTACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaaaaa              >  2:1141090/1‑149 (MQ=255)

                                                                                                                                  ttACCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaaaaaa             >  2:944895/1‑149 (MQ=255)

                                                                                                                                    aCCAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAtt           >  1:1874037/1‑149 (MQ=255)

                                                                                                                                      cAATTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTg         >  2:30122/1‑149 (MQ=255)

                                                                                                                                       aaTTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGc        >  2:1545539/1‑149 (MQ=255)

                                                                                                                                       aaTTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGc        <  1:1475083/149‑1 (MQ=255)

                                                                                                                                        aTTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa       >  1:2242235/1‑149 (MQ=255)

                                                                                                                                        aTTAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa       <  1:1195646/149‑1 (MQ=255)

                                                                                                                                         ttAGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCAt      <  1:3084152/149‑1 (MQ=255)

                                                                                                                                           aGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATgg    >  2:543240/1‑149 (MQ=255)

                                                                                                                                           aGAATAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGCTTCATACTATGGCGGTAAAAAAATTTGCATgg    >  1:1244581/1‑149 (MQ=255)

                                                                                                                                             aaTAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTa  >  2:45015/1‑149 (MQ=255)

                                                                                                                                             aaTAAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTa  <  1:3544713/149‑1 (MQ=255)

TGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTA  >  CP009273/2869636‑2869925

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