BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,070,301	C→G	7.3%	G325A (GGT→GCT) ‡	PP_4465 ←	porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,070,301	0	C	G	7.3%	652.3 / 28.4	199	G325A (GGT→GCT) ‡	PP_4465	porin
Reads supporting (aligned to +/- strand): ref base C (86/98); new base G (4/11); total (90/109)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

CCCCGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCATT  >  NC_002947/5070153‑5070449

ccccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacc                                                                                                                                                     >  2:1447704/1‑150 (MQ=255)

 cccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacct                                                                                                                                                    >  1:6053125/1‑150 (MQ=255)

 cccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacct                                                                                                                                                    >  2:3950924/1‑150 (MQ=255)

 cccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacct                                                                                                                                                    >  1:8837494/1‑150 (MQ=255)

  ccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacctc                                                                                                                                                   >  1:4064817/1‑150 (MQ=255)

  ccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacctc                                                                                                                                                   >  2:4024089/1‑150 (MQ=255)

   cGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacctcg                                                                                                                                                  >  1:4633329/1‑150 (MQ=255)

    gCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacctcga                                                                                                                                                 <  1:7837505/150‑1 (MQ=255)

    gCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacctcga                                                                                                                                                 >  1:337525/1‑150 (MQ=255)

       tttAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCa                                                                                                                                              <  2:4673264/150‑1 (MQ=255)

          aGGCCTGGAACACCTACGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTCTTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaa                                                                                                                                           <  2:3006654/150‑4 (MQ=255)

          aGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaa                                                                                                                                           <  2:4955563/150‑4 (MQ=255)

          aGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaa                                                                                                                                           <  2:7588005/150‑4 (MQ=255)

          aGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATAGAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaa                                                                                                                                           <  1:8711157/150‑4 (MQ=255)

           ggCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaag                                                                                                                                          <  1:2697930/150‑5 (MQ=255)

           ggCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaag                                                                                                                                          <  1:406516/150‑5 (MQ=255)

             ccTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAAcc                                                                                                                                        <  1:4006878/150‑1 (MQ=255)

                 gAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGt                                                                                                                                    <  2:2708236/150‑1 (MQ=255)

                 gAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGt                                                                                                                                    >  1:4394952/1‑150 (MQ=255)

                  aaCACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGtt                                                                                                                                   <  1:1874608/150‑1 (MQ=255)

                  aaCACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGtt                                                                                                                                   <  2:1036599/150‑1 (MQ=255)

                   acacCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTg                                                                                                                                  >  2:519208/1‑150 (MQ=255)

                      ccTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGcta                                                                                                                               <  2:1812947/150‑1 (MQ=255)

                       cTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGctac                                                                                                                              >  2:6584535/1‑150 (MQ=255)

                       cTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGctac                                                                                                                              >  1:52308/1‑150 (MQ=255)

                        tAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGctact                                                                                                                             >  1:9132232/1‑150 (MQ=255)

                                 gAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGt                                                                                                                    <  1:1606454/150‑1 (MQ=255)

                                  aaaTCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGGTATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGtt                                                                                                                   <  2:2827324/150‑1 (MQ=255)

                                  aaaTCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGtt                                                                                                                   <  2:4781542/150‑1 (MQ=255)

                                    aTCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGc                                                                                                                 <  1:2914945/150‑1 (MQ=255)

                                       gTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCaaa                                                                                                              >  2:2946508/1‑150 (MQ=255)

                                       gTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCaaa                                                                                                              <  2:8873419/150‑1 (MQ=255)

                                        tAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAg                                                                                                             <  2:2524362/150‑1 (MQ=255)

                                             gTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGa                                                                                                        <  2:6265684/150‑1 (MQ=255)

                                             gTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGa                                                                                                        <  2:9016500/150‑1 (MQ=255)

                                                 aGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTa                                                                                                    <  1:6796027/150‑1 (MQ=255)

                                                     cAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCAc                                                                                                >  2:7193214/1‑150 (MQ=255)

                                                      aaaaGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACt                                                                                               <  1:645842/150‑1 (MQ=255)

                                                           tAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACtt                                                                                          <  2:5670686/150‑1 (MQ=255)

                                                           tAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACtt                                                                                          >  1:2130026/1‑150 (MQ=255)

                                                             gAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACAGACttgt                                                                                        >  2:3317571/1‑150 (MQ=255)

                                                              aaTTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACttgtt                                                                                       <  2:4062394/150‑1 (MQ=255)

                                                                   tCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGc                                                                                  <  1:244433/150‑1 (MQ=255)

                                                                         gcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaaa                                                                            <  2:1987007/150‑1 (MQ=255)

                                                                         gcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaaa                                                                            <  2:1340002/150‑1 (MQ=255)

                                                                         gcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaaa                                                                            >  1:1007289/1‑150 (MQ=255)

                                                                         gcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaaa                                                                            <  1:6695235/150‑1 (MQ=255)

                                                                         gcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaaa                                                                            <  2:846294/150‑1 (MQ=255)

                                                                         gcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaaa                                                                            <  1:2570160/150‑1 (MQ=255)

                                                                         gcgcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTaaa                                                                            <  1:2052956/150‑1 (MQ=255)

                                                                           gcgcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGt                                                                          >  1:3948113/1‑150 (MQ=255)

                                                                             gcgTGAAGCTATTCACCATCACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCa                                                                        <  2:7323291/150‑1 (MQ=255)

                                                                             gcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCa                                                                        <  1:6363776/150‑1 (MQ=255)

                                                                             gcgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCa                                                                        <  2:8971965/150‑1 (MQ=255)

                                                                              cgTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAg                                                                       <  1:6902750/150‑1 (MQ=255)

                                                                               gTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAgc                                                                      <  1:9103982/150‑1 (MQ=255)

                                                                                 gAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAgcgc                                                                    >  2:2754308/1‑150 (MQ=255)

                                                                                   aGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCAt                                                                  <  1:5410367/150‑1 (MQ=255)

                                                                                    gCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATg                                                                 >  1:8958584/1‑150 (MQ=255)

                                                                                            ccATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCg                                                         >  2:7979799/1‑150 (MQ=255)

                                                                                            ccATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCg                                                         >  1:1263259/1‑150 (MQ=255)

                                                                                               tAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCtt                                                      <  2:8482927/150‑1 (MQ=255)

                                                                                                aaCGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCCGCGCATGACCACCCGCttt                                                     >  2:8771302/1‑150 (MQ=255)

                                                                                                 aCGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTta                                                    <  1:314106/150‑1 (MQ=255)

                                                                                                   gTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTtata                                                  <  1:8599656/150‑1 (MQ=255)

                                                                                                   gTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTtata                                                  <  1:4782610/150‑1 (MQ=255)

                                                                                                      ggTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGt                                               <  2:6053125/150‑1 (MQ=255)

                                                                                                       gTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTg                                              <  1:8595805/150‑1 (MQ=255)

                                                                                                           atagatatagGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAg                                          <  1:76888/150‑1 (MQ=255)

                                                                                                            tagatatagGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGc                                         <  1:806320/150‑1 (MQ=255)

                                                                                                             agatatagGGATTGGACGCCCCCTTCACCTTCTGGTCGACGTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGTCTTATCACTGTGTTGTTGTTGCCGTATAGTCAGCGCATGACCACCCGCTTTATAGGTGTATGcc                                        >  2:4630761/1‑150 (MQ=255)

                                                                                                                 atagGGATTGGACGCCCACTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAg                                    >  1:8139956/1‑150 (MQ=255)

                                                                                                                    gggATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGAcc                                 >  1:2938479/1‑150 (MQ=255)

                                                                                                                    gggATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGAcc                                 >  1:754587/1‑150 (MQ=255)

                                                                                                                       aTTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACAATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACAACCCGCTTTATAGGTGAAAGCCAAAGACCAgg                              >  1:145751/1‑150 (MQ=255)

                                                                                                                        tttgACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAggt                             <  2:6893088/147‑1 (MQ=255)

                                                                                                                        ttGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAggt                             <  1:5046534/150‑1 (MQ=255)

                                                                                                                         tGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAggtg                            <  1:3590974/150‑1 (MQ=255)

                                                                                                                              gCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtatt                       >  2:2597267/1‑150 (MQ=255)

                                                                                                                              gCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtatt                       <  2:4245675/150‑1 (MQ=255)

                                                                                                                              gCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtatt                       >  2:1156461/1‑150 (MQ=255)

                                                                                                                              gCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACGTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtatt                       >  2:2384042/1‑150 (MQ=255)

                                                                                                                                ccccTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtattat                     <  1:5980418/150‑1 (MQ=255)

                                                                                                                                ccccTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtattat                     >  2:7197595/1‑150 (MQ=255)

                                                                                                                                  ccTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCt                   >  2:3804580/1‑150 (MQ=255)

                                                                                                                                   cTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTa                  <  2:2921122/150‑1 (MQ=255)

                                                                                                                                    ttCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTAt                 >  1:5658181/1‑150 (MQ=255)

                                                                                                                                       ggttgatGGTCGAGGTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTAtttt              >  2:5762117/7‑150 (MQ=255)

                                                                                                                                       ggttgatGGTCGAGGTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTAtttt              >  2:2949135/7‑150 (MQ=255)

                                                                                                                                        ccTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTc             >  2:6996983/1‑150 (MQ=255)

                                                                                                                                         cTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTcc            >  1:7440976/1‑150 (MQ=255)

                                                                                                                                          ttCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTccc           <  1:4868638/150‑1 (MQ=255)

                                                                                                                                               gTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCATTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGGAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCAtttt      >  2:1554752/1‑150 (MQ=255)

                                                                                                                                               gTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCAtttt      >  1:5456698/1‑150 (MQ=255)

                                                                                                                                               gTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCAtttt      >  2:3241874/1‑150 (MQ=255)

                                                                                                                                               gTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCAtttt      >  2:913492/1‑150 (MQ=255)

                                                                                                                                                tcgacctcgaccATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTc     <  2:8282203/150‑1 (MQ=255)

                                                                                                                                                  gacctcgaccATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCAt   >  2:5452418/1‑150 (MQ=255)

                                                                                                                                                   acctcgaccATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCAtt  <  2:7788346/150‑1 (MQ=255)

                                                                                                                                                   acctcgaccATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCAtt  <  2:337525/150‑1 (MQ=255)

CCCCGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCATT  >  NC_002947/5070153‑5070449

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