BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,860,653	C→A	89.3%	E44* (GAA→TAA)	ibpA ←	heat shock chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,860,653	0	C	A	89.3%	485.5 / 54.1	193	E44* (GAA→TAA)	ibpA	heat shock chaperone
Reads supporting (aligned to +/- strand): ref base C (10/11); new base A (81/91); total (91/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 193 total aligned reads displayed.

TCAGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTG  >  CP009273/3860508‑3860797

tCAGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGt                                                                                                                                                <  2:1563802/148‑1 (MQ=255)

  aGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCt                                                                                                                                              <  1:2609134/148‑1 (MQ=255)

   gCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTa                                                                                                                                             <  2:2810598/148‑1 (MQ=255)

    cGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTAc                                                                                                                                            >  1:1496215/1‑148 (MQ=255)

     gATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTAcc                                                                                                                                           <  1:819637/148‑1 (MQ=255)

           cTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCa                                                                                                                                     <  1:2963873/148‑1 (MQ=255)

            tGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCaa                                                                                                                                    <  2:3247991/148‑1 (MQ=255)

              gTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACg                                                                                                                                  <  1:2320948/148‑1 (MQ=255)

                aCAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGtt                                                                                                                                <  1:668663/148‑1 (MQ=255)

                aCAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGtt                                                                                                                                <  2:3122918/148‑1 (MQ=255)

                aCAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGtt                                                                                                                                <  1:245506/148‑1 (MQ=255)

                   gATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTtata                                                                                                                             <  2:2741123/148‑1 (MQ=255)

                    aTAGGTGCGCTCTTTTTGTTCTTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATAc                                                                                                                            >  2:983067/1‑148 (MQ=255)

                    aTAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATAc                                                                                                                            >  1:192236/1‑148 (MQ=255)

                    aTAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATAc                                                                                                                            >  1:1384539/1‑148 (MQ=255)

                      aGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACgg                                                                                                                          <  1:2432884/148‑1 (MQ=255)

                       ggTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGa                                                                                                                         <  1:804494/148‑1 (MQ=255)

                        gTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAg                                                                                                                        <  1:1654180/148‑1 (MQ=255)

                        gTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAg                                                                                                                        >  2:2881464/1‑148 (MQ=255)

                          gcgcTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAggg                                                                                                                      <  2:1223575/148‑1 (MQ=255)

                          gcgcTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAggg                                                                                                                      >  2:1373359/1‑148 (MQ=255)

                             ctctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAg                                                                                                                   >  1:842942/1‑148 (MQ=255)

                              tctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgc                                                                                                                  >  1:1009969/1‑148 (MQ=255)

                              tctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgc                                                                                                                  <  1:2083803/148‑1 (MQ=255)

                              tctTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgc                                                                                                                  <  1:262180/148‑1 (MQ=255)

                                 ttttGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgccgc                                                                                                               >  2:611167/1‑148 (MQ=255)

                                  tttGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAgccgcc                                                                                                              >  1:305113/1‑148 (MQ=255)

                                     gTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCAtt                                                                                                           <  2:912472/148‑1 (MQ=255)

                                        cgtcgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTAct                                                                                                        <  2:1012590/148‑1 (MQ=255)

                                         gtcgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActc                                                                                                       <  2:1100833/148‑1 (MQ=255)

                                         gtcgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActc                                                                                                       >  2:2680359/1‑148 (MQ=255)

                                          tcgGCGTGAGCACCTGTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActct                                                                                                      <  2:2452541/148‑1 (MQ=255)

                                           cgGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctg                                                                                                     <  1:2219106/148‑1 (MQ=255)

                                            ggCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctgg                                                                                                    >  1:1940594/1‑148 (MQ=255)

                                            ggCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctgg                                                                                                    >  1:1940619/1‑148 (MQ=255)

                                                  aGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctggctctgg                                                                                              <  2:599138/148‑1 (MQ=255)

                                                  aGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctggctctgg                                                                                              <  1:2880447/148‑1 (MQ=255)

                                                  aGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctggctctgg                                                                                              >  2:2086375/1‑148 (MQ=255)

                                                  aGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTActctggctctgg                                                                                              <  2:1883939/148‑1 (MQ=255)

                                                   gCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgt                                                                                             <  1:2921751/148‑1 (MQ=255)

                                                      ccTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGgttgt                                                                                          <  1:1805782/148‑1 (MQ=255)

                                                        tttCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGttt                                                                                        >  1:2889081/1‑148 (MQ=255)

                                                         ttCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGtttt                                                                                       <  2:2643284/148‑1 (MQ=255)

                                                             ccaccaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTaa                                                                                   >  1:942872/1‑148 (MQ=255)

                                                              caccaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAg                                                                                  <  1:1013266/148‑1 (MQ=255)

                                                               accaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGt                                                                                 >  1:1295425/1‑148 (MQ=255)

                                                               accaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGt                                                                                 >  1:2366426/1‑148 (MQ=255)

                                                               accaGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGt                                                                                 >  2:1340999/1‑148 (MQ=255)

                                                                  agcagATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGt                                                                              <  1:2542411/148‑1 (MQ=255)

                                                                      gATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaa                                                                          >  2:2954832/1‑148 (MQ=255)

                                                                      gATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaa                                                                          >  1:968327/1‑148 (MQ=255)

                                                                        ttATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaaca                                                                        >  1:1345602/1‑148 (MQ=255)

                                                                        ttATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaaca                                                                        >  2:2447333/1‑148 (MQ=255)

                                                                        ttATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaaca                                                                        >  2:2021553/1‑148 (MQ=255)

                                                                         tATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGAAAAACCATCCACAGCGATAGCAATGCGGTAATGTTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaacaa                                                                       <  1:555219/148‑1 (MQ=255)

                                                                         gaTCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaaacaa                                                                       <  2:2344219/147‑1 (MQ=255)

                                                                           tCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTaa                                                                           <  2:2362440/142‑1 (MQ=255)

                                                                            ccTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACg                                                                    >  2:3082721/1‑148 (MQ=255)

                                                                            ccTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACg                                                                    >  2:3082819/1‑148 (MQ=255)

                                                                            ccTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACg                                                                    >  2:3160353/1‑148 (MQ=255)

                                                                             cTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGTTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACgg                                                                   <  1:2225710/148‑1 (MQ=255)

                                                                               gggCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAATCAAACGGTc                                                                 >  2:528236/1‑148 (MQ=255)

                                                                                    gTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACACACGGTCAAATc                                                            >  2:1306778/1‑148 (MQ=255)

                                                                                    gTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATc                                                            >  2:1618594/1‑148 (MQ=255)

                                                                                    gTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATc                                                            >  2:747940/1‑148 (MQ=255)

                                                                                      aaTTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCa                                                          >  2:2433563/1‑148 (MQ=255)

                                                                                        tttCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAAt                                                        <  1:2946328/148‑1 (MQ=255)

                                                                                        tttCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAAt                                                        <  2:3002982/148‑1 (MQ=255)

                                                                                          tCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGATCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAg                                                      >  2:219685/1‑148 (MQ=255)

                                                                                                tCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAAc                                                >  2:3109762/1‑148 (MQ=255)

                                                                                                tCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAAc                                                >  1:1227085/1‑148 (MQ=255)

                                                                                                 cGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACg                                               <  1:3007143/148‑1 (MQ=255)

                                                                                                  gCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACgg                                              >  1:502433/1‑148 (MQ=255)

                                                                                                   ctctCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGt                                             >  2:1589789/1‑148 (MQ=255)

                                                                                                   ctctCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGt                                             >  1:1544709/1‑148 (MQ=255)

                                                                                                         gCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCg                                       <  2:148818/148‑1 (MQ=255)

                                                                                                          cAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCgg                                      <  1:371800/148‑1 (MQ=255)

                                                                                                          cAAAACCAGCCACAGCGATAGCAATGAGGTAATGTTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATAACTCTGCCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGAAAAGCgg                                      >  1:2829960/1‑148 (MQ=255)

                                                                                                            aaaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCgggg                                    <  2:2018178/148‑1 (MQ=255)

                                                                                                             aaCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGa                                   >  2:1147747/1‑148 (MQ=255)

                                                                                                              aCCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGAt                                  >  1:318295/1‑148 (MQ=255)

                                                                                                                  gCCACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAAt                              <  1:2183945/148‑1 (MQ=255)

                                                                                                                   ccACAGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATc                             >  2:450648/1‑148 (MQ=255)

                                                                                                                    cacaGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCa                            <  2:1696065/148‑1 (MQ=255)

                                                                                                                       aGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAg                         <  2:1945293/148‑1 (MQ=255)

                                                                                                                       aGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAg                         <  2:795922/148‑1 (MQ=255)

                                                                                                                       aGCGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAg                         <  1:443193/148‑1 (MQ=255)

                                                                                                                         cGATAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGtt                       <  2:634412/148‑1 (MQ=255)

                                                                                                                         cGATAGCAATGCGGTAATGGTTTTAGTCTACCAGATCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGtt                       >  2:479897/1‑148 (MQ=255)

                                                                                                                            tAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACg                    >  2:967158/1‑148 (MQ=255)

                                                                                                                            tAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACg                    >  2:943715/1‑148 (MQ=255)

                                                                                                                            tAGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACg                    >  2:3004030/1‑148 (MQ=255)

                                                                                                                             aGCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGc                   <  1:1152497/148‑1 (MQ=255)

                                                                                                                              gCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCa                  <  2:2178316/148‑1 (MQ=255)

                                                                                                                              gCAATGCGGTAATGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCa                  <  1:2806653/148‑1 (MQ=255)

                                                                                                                              gCAATGCGGTAATGGTTTTAGTCTACCAGCTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCa                  >  1:2090600/1‑148 (MQ=255)

                                                                                                                                          tGGTTTTAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCt      >  1:762867/1‑148 (MQ=255)

                                                                                                                                           ggTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGAGGGGATAAATCAAAGTTACGCATAATCAATAGCTc     <  1:156509/148‑1 (MQ=255)

                                                                                                                                           ggTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGAGGGGATAAATCAAAGTTACGCATAATCAATAGCTc     <  1:156451/148‑1 (MQ=255)

                                                                                                                                              tttAGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTg  >  2:2231324/1‑148 (MQ=255)

TCAGCGATGCCCTGGTACAGATAGGTGCGCTCTTTTTGTTCGTCGGCGTGAGCACCTTTCACCACCAGCAGATTATCCTGGGCGGTAATTTCCAGTTCGCTCTCAGCAAAACCAGCCACAGCGATAGCAATGCGGTAATGGTTTTCGTCTACCAGTTCAACGTTATACGGAGGGTAGCCGCCATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTG  >  CP009273/3860508‑3860797

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