BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	613,550	G→C	5.8%	A272P (GCC→CCC) ‡	fes →	enterochelin esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	613,550	0	G	C	5.8%	398.6 / 8.9	117	A272P (GCC→CCC) ‡	fes	enterochelin esterase
Reads supporting (aligned to +/- strand): ref base G (58/51); new base C (4/3); total (62/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

AGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTA  >  NC_000913/613401‑613699

aGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTg                                                                                                                                                       >  1:1330110/1‑150 (MQ=255)

  cAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGcc                                                                                                                                                     <  2:926802/150‑1 (MQ=255)

       tcctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCg                                                                                                                                                <  1:492467/150‑1 (MQ=255)

       tcctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAACGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCg                                                                                                                                                <  1:1516503/150‑1 (MQ=255)

        cctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGc                                                                                                                                               <  1:297845/150‑1 (MQ=255)

          tccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGGCACGATCGcta                                                                                                                                             <  2:1126547/150‑3 (MQ=255)

          tccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGGCACGATCGcta                                                                                                                                             <  1:1237112/150‑3 (MQ=255)

            ccgccgTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCg                                                                                                                                           >  2:286591/1‑150 (MQ=255)

             cgccgTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGt                                                                                                                                          <  1:741130/150‑1 (MQ=255)

               ccgTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTgg                                                                                                                                        <  1:723206/150‑1 (MQ=255)

                        tgtTGATCGACGCTATCGACACCACGCACCGCTCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCGCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGc                                                                                                                               <  2:1292258/150‑1 (MQ=255)

                        tgtTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGc                                                                                                                               <  1:1174605/150‑1 (MQ=255)

                            gATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCagag                                                                                                                           >  2:1341587/1‑150 (MQ=255)

                              tCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGtt                                                                                                                         >  2:1396549/1‑150 (MQ=255)

                                gACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGtttt                                                                                                                       >  1:809385/1‑150 (MQ=255)

                                gACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGtttt                                                                                                                       <  1:294713/150‑1 (MQ=255)

                                  cGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtgg                                                                                                                     <  1:374094/150‑1 (MQ=255)

                                    cTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggtg                                                                                                                   >  1:48083/1‑150 (MQ=255)

                                    cTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggtg                                                                                                                   >  1:135012/1‑150 (MQ=255)

                                      aTCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTggg                                                                                                                 >  2:51255/1‑150 (MQ=255)

                                        cGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGCTCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCt                                                                                                               >  2:665396/1‑150 (MQ=255)

                                        cGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCt                                                                                                               >  2:34172/1‑150 (MQ=255)

                                             ccccGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCcg                                                                                                          >  1:1392712/4‑150 (MQ=255)

                                             ccacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCcg                                                                                                          >  1:909802/1‑150 (MQ=255)

                                              cacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCcgc                                                                                                         >  1:414037/1‑150 (MQ=255)

                                                   aCCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtg                                                                                                                      <  1:419021/132‑1 (MQ=255)

                                                   aCCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtg                                                                                                                      >  2:419021/1‑132 (MQ=255)

                                                   aCCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtg                                                                                                                      <  1:424200/132‑1 (MQ=255)

                                                   aCCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGt                                                                                                    <  2:276051/150‑1 (MQ=255)

                                                      gcgcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATg                                                                                                 <  1:181341/150‑1 (MQ=255)

                                                        gcCCACGAACTGCCGTGTAATGCGGAGTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGcc                                                                                               <  1:805857/150‑1 (MQ=255)

                                                         cccACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCg                                                                                              >  1:84218/1‑150 (MQ=255)

                                                         cccACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCg                                                                                              >  1:957004/1‑150 (MQ=255)

                                                          ccACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCgg                                                                                             <  1:1180606/150‑1 (MQ=255)

                                                          ccACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCgg                                                                                             <  1:1214794/150‑1 (MQ=255)

                                                              gAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTg                                                                                         <  2:84218/150‑1 (MQ=255)

                                                                 cTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCAc                                                                                      <  1:949190/150‑1 (MQ=255)

                                                                      gtgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGcc                                                                                 >  1:869890/1‑150 (MQ=255)

                                                                           aTGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAAc                                                                            <  2:1303301/150‑1 (MQ=255)

                                                                            tGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACg                                                                           >  1:1316218/1‑150 (MQ=255)

                                                                              cGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCt                                                                         >  2:1065904/1‑150 (MQ=255)

                                                                                    tCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCt                                                                   >  2:91636/1‑150 (MQ=255)

                                                                                     cTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCtg                                                                  <  1:1283717/150‑1 (MQ=255)

                                                                                     atGGCACTCAGTACAGCAACAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCtg                                                                  >  2:713816/2‑150 (MQ=255)

                                                                                      tGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCtgt                                                                 >  2:866734/1‑150 (MQ=255)

                                                                                      tGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCtgt                                                                 >  1:192764/1‑150 (MQ=255)

                                                                                       ggCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCtgtg                                                                <  2:1350768/150‑1 (MQ=255)

                                                                                         cTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTa                                                              <  2:957004/150‑1 (MQ=255)

                                                                                            gCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTa                                                           <  1:1341587/150‑1 (MQ=255)

                                                                                             cAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTaa                                                          <  2:843331/150‑1 (MQ=255)

                                                                                             cAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTaa                                                          >  2:1450336/1‑150 (MQ=255)

                                                                                              aGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAg                                                         <  1:1396549/150‑1 (MQ=255)

                                                                                              aGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAg                                                         <  2:192764/150‑1 (MQ=255)

                                                                                                  cAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAg                                                     >  1:511062/1‑150 (MQ=255)

                                                                                                   aGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGt                                                    >  1:839092/1‑150 (MQ=255)

                                                                                                   aGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGt                                                    >  1:1377168/1‑150 (MQ=255)

                                                                                                    gCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTc                                                   >  1:344989/1‑150 (MQ=255)

                                                                                                    gCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTc                                                   >  1:488997/1‑150 (MQ=255)

                                                                                                     cAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCa                                                  <  2:203590/150‑1 (MQ=255)

                                                                                                      aaGAGTTATTACCCCTTGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCAGCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCTCTCTGTATGCCTGACTGCACTGGCCTGAACGCGTTGGCTGTGTATTAAGCCAGTCAg                                                 <  2:749737/150‑1 (MQ=255)

                                                                                                      aaGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCtg                                                                  <  2:363912/133‑1 (MQ=255)

                                                                                                      aaGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAg                                                 >  2:270448/1‑150 (MQ=255)

                                                                                                        gagTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGa                                               <  1:286591/150‑1 (MQ=255)

                                                                                                          gTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGGTGTGTATTAGGCCAGGCAGGagg                                             >  1:664906/1‑148 (MQ=255)

                                                                                                           ttattaCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCg                                            <  2:1178957/150‑1 (MQ=255)

                                                                                                            tattaCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGt                                           >  2:1249054/1‑150 (MQ=255)

                                                                                                               taCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACt                                        >  1:1523918/1‑150 (MQ=255)

                                                                                                                aCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtg                                       >  1:130511/1‑150 (MQ=255)

                                                                                                                 ccccTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtgg                                      >  1:1137692/1‑150 (MQ=255)

                                                                                                                 ccccTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtgg                                      <  1:1077309/150‑1 (MQ=255)

                                                                                                                  cccTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtggt                                     >  2:363005/1‑150 (MQ=255)

                                                                                                                      ggTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGcc                                 >  2:316890/1‑150 (MQ=255)

                                                                                                                       gTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCg                                <  1:1065904/150‑1 (MQ=255)

                                                                                                                         gAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCtg                                                                  >  1:84764/1‑114 (MQ=255)

                                                                                                                           aaGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATc                            >  1:1042975/1‑150 (MQ=255)

                                                                                                                             gCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGACTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCgg                          >  1:1289163/1‑150 (MQ=255)

                                                                                                                              cTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCggg                         <  2:436533/150‑1 (MQ=255)

                                                                                                                              cTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCggg                         <  1:1369452/150‑1 (MQ=255)

                                                                                                                               tATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCgggg                        >  2:256805/1‑150 (MQ=255)

                                                                                                                                  tGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcgg                     >  1:1301615/1‑150 (MQ=255)

                                                                                                                                   gcgccTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcggg                    <  2:1441574/147‑1 (MQ=255)

                                                                                                                                   gCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcggg                    >  1:623730/1‑150 (MQ=255)

                                                                                                                                   gCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcggg                    >  2:576731/1‑150 (MQ=255)

                                                                                                                                    ccccTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcgggc                   >  1:718458/1‑150 (MQ=255)

                                                                                                                                    ccccTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcgggc                   >  2:718588/1‑150 (MQ=255)

                                                                                                                                    ccccTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcgggc                   <  1:1319198/150‑1 (MQ=255)

                                                                                                                                    ccccTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCAGCATCGgggcgggc                   >  2:1493986/1‑150 (MQ=255)

                                                                                                                                        ttttAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGgcagca               <  1:226151/150‑1 (MQ=255)

                                                                                                                                        acggtgCGATCGGCACGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGgcagca               >  1:286752/6‑150 (MQ=255)

                                                                                                                                        acggtgCGATCGGCACGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGgcagca               >  1:1122881/6‑150 (MQ=255)

                                                                                                                                         tttAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCaa              <  1:676158/150‑1 (MQ=255)

                                                                                                                                           tAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAaga            <  2:414037/150‑1 (MQ=255)

                                                                                                                                           gtgCGATCGGCACGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAaga            >  1:433652/3‑150 (MQ=255)

                                                                                                                                            aGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAagag           >  2:291487/1‑150 (MQ=255)

                                                                                                                                              cGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAggg         <  2:1377168/150‑1 (MQ=255)

                                                                                                                                               gATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGc        <  2:1325676/150‑1 (MQ=255)

                                                                                                                                               gATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGc        >  2:1135975/1‑150 (MQ=255)

                                                                                                                                                   gTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCgtgt    <  2:514131/150‑1 (MQ=255)

                                                                                                                                                   gTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCgtgt    >  2:956827/1‑150 (MQ=255)

                                                                                                                                                     gCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTa  >  1:702197/1‑150 (MQ=255)

AGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTA  >  NC_000913/613401‑613699

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