BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,859,994	A→T	6.8%	L88L (CTT→CTA)	ibpB ←	heat shock chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,859,994	0	A	T	6.8%	1291.8 / 61.0	367	L88L (CTT→CTA)	ibpB	heat shock chaperone
Reads supporting (aligned to +/- strand): ref base A (165/176); new base T (16/9); total (181/186)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 371 total aligned reads displayed.

GGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTT  >  CP009273/3859845‑3860143

 gACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATaa                                                                                                                                                      <  1:582104/150‑1 (MQ=255)

  aCGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAg                                                                                                                                                     <  2:1344732/150‑1 (MQ=255)

   cggtCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGc                                                                                                                                                    <  2:4466382/147‑1 (MQ=255)

    gTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATTAGcc                                                                                                                                                   >  1:2433128/1‑150 (MQ=255)

     ttCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGccc                                                                                                                                                  <  1:1122704/150‑1 (MQ=255)

       cGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGATTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCtt                                                                                                                                                <  1:4135201/150‑1 (MQ=255)

          tGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGAt                                                                                                                                             <  2:4212407/150‑1 (MQ=255)

          tGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATGACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGAt                                                                                                                                             <  2:2475140/150‑1 (MQ=255)

             tAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCa                                                                                                                                          >  1:4888857/1‑150 (MQ=255)

               gCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGc                                                                                                                                        <  1:5179400/150‑1 (MQ=255)

                 gATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATTAGCCCTTGATGCAGCCa                                                                                                                                      <  2:3080660/150‑1 (MQ=255)

                  aTACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCAt                                                                                                                                     <  2:1063690/150‑1 (MQ=255)

                    aCGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCAttt                                                                                                                                   <  1:1214623/150‑1 (MQ=255)

                     cgctgcgctgCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCAtttt                                                                                                                                  >  2:288091/1‑150 (MQ=255)

                     cgctgcgctgCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCAtttt                                                                                                                                  >  2:288245/1‑150 (MQ=255)

                       ctgcgctgcgATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCAtttttt                                                                                                                                <  2:3323251/150‑1 (MQ=255)

                        tgcgctgcgATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTtc                                                                                                                               >  1:3393750/1‑150 (MQ=255)

                        tgcgctgcgATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCATTGATGCAGCCATTTTTtc                                                                                                                               >  2:3703777/1‑150 (MQ=255)

                           gctgcgATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTctct                                                                                                                            >  1:651459/1‑150 (MQ=255)

                            ctgcgATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCtctt                                                                                                                           <  1:2024044/150‑1 (MQ=255)

                              gcgATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCtcttct                                                                                                                         <  2:827273/150‑1 (MQ=255)

                               cgATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTcttctt                                                                                                                        <  1:3822159/150‑1 (MQ=255)

                                gATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCttt                                                                                                                       >  1:4594009/1‑150 (MQ=255)

                                 aTGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCtttt                                                                                                                      >  2:728554/1‑150 (MQ=255)

                                    ggTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGgc                                                                                                                   <  1:617697/150‑1 (MQ=255)

                                    ggTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGgc                                                                                                                   >  1:4490554/1‑150 (MQ=255)

                                      ttCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGgctg                                                                                                                 >  2:4411862/1‑150 (MQ=255)

                                        cAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGgctgct                                                                                                               >  2:32281/1‑150 (MQ=255)

                                        cAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGgctgct                                                                                                               >  2:873380/1‑150 (MQ=255)

                                         aGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTc                                                                                                              <  1:3083346/150‑1 (MQ=255)

                                         aGGCTCATTACGAATTAAATCAATACGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTc                                                                                                              <  2:2148556/150‑1 (MQ=255)

                                              cATTACGAATTAAATAAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATACATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCg                                                                                                         >  2:2695578/1‑150 (MQ=255)

                                               aTTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGt                                                                                                        >  1:972298/1‑150 (MQ=255)

                                                    gAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCtt                                                                                                   >  2:4291642/1‑150 (MQ=255)

                                                     aaTTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCttt                                                                                                  >  2:3521363/1‑150 (MQ=255)

                                                     aaTTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCttt                                                                                                  >  2:3871553/1‑150 (MQ=255)

                                                       ttAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTa                                                                                                <  2:1375830/150‑1 (MQ=255)

                                                             cAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCa                                                                                          >  1:5085811/1‑150 (MQ=255)

                                                              aaTATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAg                                                                                         <  1:2826813/150‑1 (MQ=255)

                                                                tatGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGgc                                                                                       >  1:2171772/1‑150 (MQ=255)

                                                                tatGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGgc                                                                                       <  1:1036060/150‑1 (MQ=255)

                                                                           ccGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATATGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTcc                                                                            >  1:1718220/1‑150 (MQ=255)

                                                                             gTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAg                                                                          <  1:1085394/150‑1 (MQ=255)

                                                                               ttACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGtt                                                                        >  1:2313369/1‑150 (MQ=255)

                                                                                 aCGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCGTTTGGCTGCTCCGGCGTGCCTTTTAAGCTCAGGCGCGTACCTTCCAGTTGa                                                                      >  2:3814789/1‑150 (MQ=255)

                                                                                   gAAGGTTGCTCCAGAGCATTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAATTGGCTGATTCATAAGCCCTTGTTGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAAt                                                                    <  1:2569602/150‑1 (MQ=255)

                                                                                    aaGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATTAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGACTTTTACGCTCAGGCGCGTACCTTCCAGTTGAAtc                                                                   >  2:117447/1‑150 (MQ=255)

                                                                                    aaGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATTAGCCCTTGATGCAGCCATTTTTTCTCGTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAAtc                                                                   >  2:4858649/1‑150 (MQ=255)

                                                                                    aaGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAAtc                                                                   <  1:3586446/150‑1 (MQ=255)

                                                                                     aGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTACCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAAtct                                                                  >  1:1657602/1‑150 (MQ=255)

                                                                                      ggTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTTCCTTCCAGTTGAAtctc                                                                 >  2:3197579/1‑150 (MQ=255)

                                                                                       gTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATctct                                                                >  1:3662423/1‑150 (MQ=255)

                                                                                          gcgcCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTaaa                                                             <  2:3307381/150‑1 (MQ=255)

                                                                                           cgcCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTACAt                                                            >  2:2306689/1‑150 (MQ=255)

                                                                                             ccAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCt                                                          >  2:3999495/1‑150 (MQ=255)

                                                                                              cAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCtt                                                         >  2:3790339/1‑150 (MQ=255)

                                                                                              cAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCtt                                                         >  2:4340430/1‑150 (MQ=255)

                                                                                               agagACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATTAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTc                                                        >  1:3971617/1‑150 (MQ=255)

                                                                                                 agaCTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCt                                                      <  1:2728103/150‑1 (MQ=255)

                                                                                                    cTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGAc                                                   >  2:4119450/1‑150 (MQ=255)

                                                                                                    cTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGAc                                                   <  1:1891009/150‑1 (MQ=255)

                                                                                                     ttCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACg                                                  >  1:3734818/1‑150 (MQ=255)

                                                                                                     ttCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACg                                                  >  2:518253/1‑150 (MQ=255)

                                                                                                     ttCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTGTTGCTCTACTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAGTCTCTAAATCTTCCTGACg                                                  >  2:2807278/1‑150 (MQ=255)

                                                                                                        cATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGTCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATATTCCTGACGGaa                                               <  1:3309586/150‑1 (MQ=255)

                                                                                                        cATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGaa                                               <  1:427617/150‑1 (MQ=255)

                                                                                                         atatTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGaaa                                              <  1:1541225/150‑1 (MQ=255)

                                                                                                           atTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAAcc                                            >  1:3551061/1‑150 (MQ=255)

                                                                                                            ttttCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCt                                           <  1:3460052/150‑1 (MQ=255)

                                                                                                              ttCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAATTGAATCTCTAAATCTTCCTGACGGAAACCTGc                                         <  1:4791060/150‑1 (MQ=255)

                                                                                                                cagccagcGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCa                                       >  1:787337/1‑150 (MQ=255)

                                                                                                                  gccagcGTAAAGCTCAGGCTAAATGTCTGATTCATAAGCCCTTGATGCAGCAATTTTTTCTCTTCTTTTGGAAGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCAAGTTGAATCTCTAAATCTACCAGACTGAAACCTGCCAgc                                     >  1:4341474/1‑150 (MQ=255)

                                                                                                                  gccagcGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAgc                                     >  1:747544/1‑150 (MQ=255)

                                                                                                                  gccagcGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAgc                                     >  1:156827/1‑150 (MQ=255)

                                                                                                                  gccagcGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCGTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAAGCTTCCTGACGGAAACCTGCCAgc                                     >  2:4580687/1‑150 (MQ=255)

                                                                                                                    cagcGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAgcgc                                   >  2:1883937/1‑150 (MQ=255)

                                                                                                                     agcGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCa                                  <  1:4198917/150‑1 (MQ=255)

                                                                                                                        gTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAgg                               <  1:2253824/150‑1 (MQ=255)

                                                                                                                         tAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGTTGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAggg                              <  2:4098553/150‑1 (MQ=255)

                                                                                                                         tAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAggg                              >  2:1250561/1‑150 (MQ=255)

                                                                                                                            aGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTaa                           >  1:1564060/1‑150 (MQ=255)

                                                                                                                              cTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATg                         <  1:3194490/150‑1 (MQ=255)

                                                                                                                              cTCAGGCTAAATGGCTGAGTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATg                         >  2:1978077/1‑150 (MQ=255)

                                                                                                                                  ggCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGt                     <  1:2788723/150‑1 (MQ=255)

                                                                                                                                  ggCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTATTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGt                     <  2:2376285/150‑1 (MQ=255)

                                                                                                                                    cTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAg                   <  2:1443495/150‑1 (MQ=255)

                                                                                                                                     tAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGt                  <  1:1386832/150‑1 (MQ=255)

                                                                                                                                        aTGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCACGGGTAATGCGGTAGTGGt               >  2:4553546/1‑150 (MQ=255)

                                                                                                                                        aTGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGt               <  1:1214394/150‑1 (MQ=255)

                                                                                                                                        aTGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGt               <  2:3181793/150‑1 (MQ=255)

                                                                                                                                         tGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGtt              <  1:2929302/150‑1 (MQ=255)

                                                                                                                                          ggCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCGTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTa             <  2:603285/150‑1 (MQ=255)

                                                                                                                                            cTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTAtc           <  1:375575/150‑1 (MQ=255)

                                                                                                                                            cTGAGTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTAtc           <  2:5047306/150‑1 (MQ=255)

                                                                                                                                             tGATTCATTAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCCTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTAtcg          <  1:4307638/150‑1 (MQ=255)

                                                                                                                                              gATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTAtcgt         >  1:2635393/1‑150 (MQ=255)

                                                                                                                                              gATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTAtcgt         >  1:91779/1‑150 (MQ=255)

                                                                                                                                               aTTCATTAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTAtcgtc        <  2:1704151/150‑1 (MQ=255)

                                                                                                                                               aTTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTAtcgtc        <  2:1525941/150‑1 (MQ=255)

                                                                                                                                                ttCATTAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTAtcgtcg       <  2:1770167/150‑1 (MQ=255)

GGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAGGCTAAATGGCTGATTCATAAGCCCTTGATGCAGCCATTTTTTCTCTTCTTTTGGCTGCTCCGGCGTGCCTTTTACGCTCAGGCGCGTACCTTCCAGTTGAATCTCTAAATCTTCCTGACGGAAACCTGCCAGCGCAAGGGTAATGCGGTAGTGGTTATCGTCGCTTTT  >  CP009273/3859845‑3860143

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