BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,992,585	G→T	85.7%	S478I (AGC→ATC)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,992,585	0	G	T	85.7%	277.6 / 44.0	105	S478I (AGC→ATC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base G (9/6); new base T (42/47); total (52/53)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

Only 100 of 118 total aligned reads displayed.

ACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGA  >  NC_000913/3992439‑3992729

aCCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATcatcatcc                                                                                                                                                <  2:548707/149‑1 (MQ=255)

  cAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCGCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCTATcatcatccat                                                                                                                                              >  2:886775/1‑149 (MQ=255)

   aGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATcatcagccatc                                                                                                                                             <  3:282654/149‑1 (MQ=255)

    ggTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTAGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATcatcatccatca                                                                                                                                            <  3:691922/149‑1 (MQ=255)

     gTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATcatcatccatcag                                                                                                                                           >  2:714245/1‑149 (MQ=255)

      tAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATcatcatccatcagc                                                                                                                                          >  4:608342/1‑149 (MQ=255)

      tAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATcatcatccatcagc                                                                                                                                          <  3:880357/149‑1 (MQ=255)

          gTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCt                                                                                                                                      >  2:212278/1‑149 (MQ=255)

          gTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCt                                                                                                                                      >  4:722587/1‑149 (MQ=255)

           tGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTg                                                                                                                                     <  1:411099/149‑1 (MQ=255)

            gACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTgg                                                                                                                                    >  2:249657/1‑149 (MQ=255)

                   gTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATaa                                                                                                                             >  3:70754/1‑149 (MQ=255)

                     aaaCCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAAcc                                                                                                                           >  3:221822/1‑149 (MQ=255)

                     aaaCCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAAcc                                                                                                                           >  4:621158/1‑149 (MQ=255)

                      aaCCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTTGATCATCATCCATCAGCCGCTGGAATATAACCg                                                                                                                          >  3:652683/1‑149 (MQ=255)

                         ccGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATcatcatccatcagc                                                                                                                                          <  2:623779/130‑1 (MQ=255)

                         ccGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATcatcatccatcagc                                                                                                                                          >  1:623779/1‑130 (MQ=255)

                         ccGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTa                                                                                                                       >  3:409472/1‑149 (MQ=255)

                          cGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTAc                                                                                                                      >  1:389474/1‑149 (MQ=255)

                          cGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTAc                                                                                                                      >  4:387186/1‑149 (MQ=255)

                          cGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGACTATACCCGTTAc                                                                                                                      >  4:108316/1‑149 (MQ=255)

                           gCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTAcc                                                                                                                     <  4:235842/149‑1 (MQ=255)

                                    aCCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAAc                                                                                                            >  2:363849/1‑149 (MQ=255)

                                     cccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACt                                                                                                           >  3:274536/1‑149 (MQ=255)

                                      ccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATaa                                                                                                                             <  3:191720/130‑1 (MQ=255)

                                      ccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATaa                                                                                                                             >  4:191720/1‑130 (MQ=255)

                                      ccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGtt                                                                                                                        >  3:376516/1‑135 (MQ=255)

                                      ccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGtt                                                                                                                        <  4:376516/135‑1 (MQ=255)

                                      ccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtg                                                                                                          >  1:555298/1‑149 (MQ=255)

                                        gATCTCTCGGAACCGAATCTGATCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCTCGCGCCAAATATTGAGTCGATCATCATCAACCAGCCGCTGGAATATAACCGTTACCTGAATAAACtggt                                                                                                        <  4:762425/149‑1 (MQ=255)

                                        gATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtggt                                                                                                        <  2:387228/149‑1 (MQ=255)

                                           ctctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGc                                                                                                     <  1:429241/149‑1 (MQ=255)

                                            tctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCg                                                                                                    >  4:46228/1‑149 (MQ=255)

                                            tctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCAGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCg                                                                                                    >  2:391485/1‑149 (MQ=255)

                                             ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGt                                                                                                   <  1:212278/149‑1 (MQ=255)

                                             ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGt                                                                                                   <  2:389474/149‑1 (MQ=255)

                                              tcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTg                                                                                                  <  1:374356/149‑1 (MQ=255)

                                                ggAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTggg                                                                                                <  3:71396/149‑1 (MQ=255)

                                                 gAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGc                                                                                               >  2:517668/1‑149 (MQ=255)

                                                     cGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATgg                                                                                           <  3:608342/149‑1 (MQ=255)

                                                       aaTCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGtt                                                                                         <  2:108483/149‑1 (MQ=255)

                                                         tCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATg                                                                                            >  1:308550/1‑144 (MQ=255)

                                                         tCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATg                                                                                            <  2:308550/144‑1 (MQ=255)

                                                         tCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTa                                                                                       <  4:494320/149‑1 (MQ=255)

                                                          cTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTaa                                                                                      <  1:402087/149‑1 (MQ=255)

                                                           tGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAAc                                                                                     <  3:621158/149‑1 (MQ=255)

                                                           tGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAAc                                                                                     <  2:698566/149‑1 (MQ=255)

                                                            gACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACg                                                                                    <  1:4524/149‑1 (MQ=255)

                                                            gACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACg                                                                                    <  1:914385/149‑1 (MQ=255)

                                                              ccTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGc                                                                                  <  3:123175/149‑1 (MQ=255)

                                                               cTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGGCAAATATTGAGTCGATCATCATCCATCAGCCGGTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGcc                                                                                 >  4:882323/1‑149 (MQ=255)

                                                               cTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGcc                                                                                 >  3:4106/1‑149 (MQ=255)

                                                                tttatttaTGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCct                                                                                >  3:194329/1‑149 (MQ=255)

                                                                 ttatttatGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGGGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCctg                                                                               >  1:679834/1‑149 (MQ=255)

                                                                   atttatGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCctgct                                                                             <  2:858124/149‑1 (MQ=255)

                                                                               ccgggccgggCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCAc                                                                 <  4:4106/149‑1 (MQ=255)

                                                                                     cgggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGttt                                                           <  1:870907/149‑1 (MQ=255)

                                                                                     cgggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGttt                                                           <  3:858563/149‑1 (MQ=255)

                                                                                     cgggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGttt                                                           >  4:838893/1‑149 (MQ=255)

                                                                                        gcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGta                                                        <  3:224201/149‑1 (MQ=255)

                                                                                         cTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCATATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGaat                                                       <  4:70754/149‑1 (MQ=255)

                                                                                         cTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAACATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGGGAACCCGTTTGtat                                                       >  4:535075/1‑149 (MQ=255)

                                                                                               gTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTTGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                 <  4:652683/149‑1 (MQ=255)

                                                                                               gTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                 <  4:397079/149‑1 (MQ=255)

                                                                                               gTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                 <  3:332798/149‑1 (MQ=255)

                                                                                               gTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                 >  4:5410/1‑149 (MQ=255)

                                                                                                ttCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAg                                                <  3:698471/149‑1 (MQ=255)

                                                                                                ttCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAg                                                >  1:503280/1‑149 (MQ=255)

                                                                                                     gTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGt                                                                                          <  4:170347/102‑1 (MQ=255)

                                                                                                     gTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGt                                                                                          >  3:170347/1‑102 (MQ=255)

                                                                                                     gTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTaa                                                                                      >  4:230807/1‑106 (MQ=255)

                                                                                                      ttGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACg                                          >  4:752452/1‑149 (MQ=255)

                                                                                                       tGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACgg                                         <  1:262811/149‑1 (MQ=255)

                                                                                                           aTCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCAtt                                     >  2:303238/1‑149 (MQ=255)

                                                                                                           aTCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCAtt                                     >  2:582630/1‑149 (MQ=255)

                                                                                                                tataACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGa                                >  4:89103/1‑149 (MQ=255)

                                                                                                                tataACCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCAGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGa                                >  2:521517/1‑149 (MQ=255)

                                                                                                                tataACCGCGCGCCAAATATTGAGTCGATCATCACCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGa                                >  1:46353/1‑149 (MQ=255)

                                                                                                                   aaCCGCGCGCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGAttt                             <  1:363849/149‑1 (MQ=255)

                                                                                                                      cgcgcgCCAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGcc                          >  4:823167/1‑149 (MQ=255)

                                                                                                                             cAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGt                                   <  2:490352/133‑1 (MQ=255)

                                                                                                                             cAAATATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTg                   <  3:553644/149‑1 (MQ=255)

                                                                                                                             cAAATATTGAGACGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGt                                   >  1:490352/1‑133 (MQ=255)

                                                                                                                               aaTATTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCa                 <  3:89103/149‑1 (MQ=255)

                                                                                                                               aaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCa                 <  2:489464/149‑1 (MQ=255)

                                                                                                                                atatTGAGTCGATGATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTCACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAg                <  3:535075/149‑1 (MQ=255)

                                                                                                                                atatTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAg                <  3:901587/149‑1 (MQ=255)

                                                                                                                                atatTGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAg                <  3:126517/149‑1 (MQ=255)

                                                                                                                                    tGAGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGgaga            <  4:194329/149‑1 (MQ=255)

                                                                                                                                     gAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGAt           <  1:582630/149‑1 (MQ=255)

                                                                                                                                      aGTCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATg          <  4:415977/149‑1 (MQ=255)

                                                                                                                                        tCGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGt        >  4:830709/1‑149 (MQ=255)

                                                                                                                                        tCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGt        <  3:722587/149‑1 (MQ=255)

                                                                                                                                         cGATCATCATCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTc       <  1:517668/149‑1 (MQ=255)

                                                                                                                                           atcatcatcCATCAGCCGCTGGAATATAACCGTTACCTGAATGAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGc     >  3:544324/1‑149 (MQ=255)

                                                                                                                                           atcatcatcCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGc     >  4:595846/1‑149 (MQ=255)

                                                                                                                                           atcatcatcCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGc     >  1:648711/1‑149 (MQ=255)

                                                                                                                                           atcatcAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGc     >  3:883066/1‑149 (MQ=255)

                                                                                                                                            tcatcaTCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGcc    >  1:240180/1‑149 (MQ=255)

                                                                                                                                              atcatccatcagccGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGa  >  2:420117/1‑149 (MQ=255)

ACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGA  >  NC_000913/3992439‑3992729

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