BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,618,536	C→A	100%	R267R (CGG→AGG)	yhhI →	putative transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,618,536	0	C	A	97.6%	664.5 / 8.2	246	R267R (CGG→AGG)	yhhI	putative transposase
Reads supporting (aligned to +/- strand): ref base C (2/4); new base A (130/110); total (132/114)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 246 total aligned reads displayed.

AGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACG  >  CP009273/3618393‑3618680

    tGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCata                                                                                                                                          >  1:1321263/1‑148 (MQ=17)

    tGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCata                                                                                                                                          >  2:1157747/1‑148 (MQ=17)

    tGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTGTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCata                                                                                                                                          <  2:2866033/148‑1 (MQ=17)

     gACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCataa                                                                                                                                         <  2:539077/148‑1 (MQ=17)

       cAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCataata                                                                                                                                       <  1:1055436/148‑1 (MQ=17)

          ttACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCa                                                                                                                                    >  1:632337/1‑148 (MQ=17)

           tACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAg                                                                                                                                   <  2:3301096/148‑1 (MQ=17)

            aCGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGa                                                                                                                                  <  1:994145/148‑1 (MQ=17)

             cGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGaa                                                                                                                                 >  1:3225860/1‑148 (MQ=17)

             cGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGaa                                                                                                                                 >  1:1759204/1‑148 (MQ=17)

               cAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACa                                                                                                                               <  2:1509324/148‑1 (MQ=17)

                  tCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaaaa                                                                                                                            >  2:463421/1‑148 (MQ=14)

                   cAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAaaag                                                                                                                           <  1:1000907/148‑1 (MQ=11)

                           agagTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGcc                                                                                                                   >  2:885457/1‑148 (MQ=11)

                               tCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGaa                                                                                                               <  2:323590/148‑1 (MQ=18)

                               tCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGaa                                                                                                               <  2:1293732/148‑1 (MQ=18)

                                cACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGaaa                                                                                                              >  2:2437939/1‑148 (MQ=21)

                                cACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGaaa                                                                                                              >  1:3208960/1‑148 (MQ=21)

                                 aCGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAAt                                                                                                             >  2:2738063/1‑148 (MQ=21)

                                   ggCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGa                                                                                                           <  2:2420982/148‑1 (MQ=21)

                                      agagAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACgg                                                                                                        <  1:2073658/148‑1 (MQ=21)

                                      agagAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACgg                                                                                                        >  1:1945468/1‑148 (MQ=21)

                                      agagAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACgg                                                                                                        <  1:3324950/148‑1 (MQ=21)

                                       gagaAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAATAGCCAGAAATGACGGt                                                                                                       <  2:1420910/148‑1 (MQ=21)

                                       gagaAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGt                                                                                                       <  1:2302259/148‑1 (MQ=21)

                                       gagaAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGt                                                                                                       >  1:291609/1‑148 (MQ=21)

                                       gagaAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGt                                                                                                       >  1:1300296/1‑148 (MQ=21)

                                        agaagaAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTc                                                                                                      <  2:700746/148‑1 (MQ=21)

                                        agaagaAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTc                                                                                                      <  2:367581/148‑1 (MQ=21)

                                         gaagaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCa                                                                                                     >  1:165310/1‑148 (MQ=21)

                                            gaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGat                                                                                                  <  2:2194686/148‑1 (MQ=21)

                                             aaaTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGata                                                                                                 >  2:2314021/1‑148 (MQ=21)

                                             aaaTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGata                                                                                                 <  1:1468911/148‑1 (MQ=21)

                                             aaaTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGata                                                                                                 >  2:3361946/1‑148 (MQ=21)

                                               aTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                                               >  2:1720745/1‑148 (MQ=21)

                                               aTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCAAGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                                               <  2:3427856/148‑1 (MQ=21)

                                                tCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatta                                                                                              <  1:1815042/148‑1 (MQ=21)

                                                 ccGTCTTCATATTGTTTGCGATTTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtattat                                                                                             >  1:1118329/1‑148 (MQ=21)

                                                 ccGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtattat                                                                                             >  1:3317934/1‑148 (MQ=21)

                                                  cGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATtata                                                                                            >  1:2818401/1‑148 (MQ=34)

                                                  cGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATtata                                                                                            >  1:3362319/1‑148 (MQ=21)

                                                   gTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTatat                                                                                           <  2:801754/148‑1 (MQ=21)

                                                   gTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTatat                                                                                           >  1:702918/1‑148 (MQ=21)

                                                    tCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATc                                                                                          <  2:3250541/148‑1 (MQ=21)

                                                       tCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGt                                                                                       <  2:1934297/148‑1 (MQ=21)

                                                       tCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGt                                                                                       <  2:491089/148‑1 (MQ=21)

                                                       tCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGt                                                                                       >  1:2866604/1‑148 (MQ=21)

                                                        cATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGtt                                                                                      <  2:3236694/148‑1 (MQ=21)

                                                          tatTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTct                                                                                    >  2:1341207/1‑148 (MQ=21)

                                                            ttGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctgc                                                                                  <  2:899296/148‑1 (MQ=21)

                                                             tGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctgct                                                                                 >  1:868017/1‑148 (MQ=21)

                                                                    gATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTaa                                                                          <  2:2370005/148‑1 (MQ=21)

                                                                     aTGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGACAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAAc                                                                         >  1:1975844/1‑148 (MQ=21)

                                                                       gTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCg                                                                       <  1:3312896/148‑1 (MQ=21)

                                                                       gTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCg                                                                       <  1:3353909/148‑1 (MQ=21)

                                                                          ccTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTg                                                                    >  2:1114558/1‑148 (MQ=25)

                                                                             gatgaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTgaga                                                                 <  2:2626324/148‑1 (MQ=37)

                                                                              atgaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGaa                                                                <  2:601611/148‑1 (MQ=37)

                                                                                  aCTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTc                                                            >  2:1055463/1‑148 (MQ=255)

                                                                                   cTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCg                                                           >  1:1467503/1‑148 (MQ=255)

                                                                                   cTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCg                                                           >  1:1803549/1‑148 (MQ=255)

                                                                                   cTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCg                                                           >  1:222811/1‑148 (MQ=255)

                                                                                   cTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCg                                                           >  1:239488/1‑148 (MQ=255)

                                                                                    ttATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGc                                                          <  1:463421/148‑1 (MQ=255)

                                                                                      aTTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCca                                                        <  1:2525984/148‑1 (MQ=255)

                                                                                       ttGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCcac                                                       <  1:896603/148‑1 (MQ=255)

                                                                                         gATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAg                                                     <  1:2146739/148‑1 (MQ=255)

                                                                                           tttCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCg                                                   >  1:3245899/1‑148 (MQ=255)

                                                                                             tCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAATAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGAt                                                 >  2:1133001/1‑148 (MQ=255)

                                                                                             tCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGAt                                                 <  2:3363507/148‑1 (MQ=255)

                                                                                                cGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTCGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCg                                              >  1:895009/1‑148 (MQ=255)

                                                                                                     gAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATc                                         <  2:2151144/148‑1 (MQ=255)

                                                                                                     gAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATc                                         <  2:2321934/148‑1 (MQ=255)

                                                                                                       aTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCAc                                       >  1:1023229/1‑148 (MQ=255)

                                                                                                       aTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCAc                                       >  1:3342234/1‑148 (MQ=255)

                                                                                                        tGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACt                                      <  1:1648280/148‑1 (MQ=255)

                                                                                                          gAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTgg                                    >  1:2960477/1‑148 (MQ=255)

                                                                                                           aaaGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGc                                   >  1:1688008/1‑148 (MQ=255)

                                                                                                             aGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACAGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCAc                                 >  2:3039241/1‑148 (MQ=255)

                                                                                                              gggCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACg                                >  1:2305195/1‑148 (MQ=255)

                                                                                                               ggCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGt                               >  2:1555794/1‑148 (MQ=255)

                                                                                                                gCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTg                              <  1:227902/148‑1 (MQ=255)

                                                                                                                 cTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTgg                             >  1:1088966/1‑148 (MQ=255)

                                                                                                                 cTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTgg                             >  1:2026164/1‑148 (MQ=255)

                                                                                                                  tGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGga                            <  1:506301/148‑1 (MQ=255)

                                                                                                                  tGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGga                            >  1:2114717/1‑148 (MQ=255)

                                                                                                                   gaagaaATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgag                           <  2:2818401/148‑1 (MQ=255)

                                                                                                                    aagaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgaga                          <  2:1834046/148‑1 (MQ=255)

                                                                                                                    aagaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgaga                          <  2:3473661/148‑1 (MQ=255)

                                                                                                                        aaTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATaa                      <  2:2031277/148‑1 (MQ=255)

                                                                                                                         aTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAg                     <  2:1849012/148‑1 (MQ=255)

                                                                                                                          ttATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGc                    <  1:2676346/148‑1 (MQ=255)

                                                                                                                              gCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCa                >  2:3007376/1‑148 (MQ=255)

                                                                                                                               cGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCAc               >  1:3143319/1‑148 (MQ=255)

                                                                                                                                 tGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTg             <  2:2254779/148‑1 (MQ=255)

                                                                                                                                 tGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTg             >  2:2680194/1‑148 (MQ=255)

                                                                                                                                  ggCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTgg            >  1:299211/1‑148 (MQ=255)

                                                                                                                                    cAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCg          >  1:436581/1‑148 (MQ=255)

                                                                                                                                        ctcCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg      >  1:547770/1‑148 (MQ=255)

                                                                                                                                            tttAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACg  >  1:1575539/1‑148 (MQ=255)

AGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACG  >  CP009273/3618393‑3618680

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