BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,412,848	T→G	5.1%	G154G (GGT→GGG)	yjfM →	DUF1190 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,412,848	0	T	G	5.1%	625.6 / 3.8	203	G154G (GGT→GGG)	yjfM	DUF1190 family protein
Reads supporting (aligned to +/- strand): ref base T (62/126); new base G (13/2); total (75/128)
Fisher's exact test for biased strand distribution p-value = 5.98e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

AATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTC  >  NC_000913/4412699‑4412997

aaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTGCTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGt                                                                                                                                                       <  1:2628493/150‑1 (MQ=255)

aaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTTGGTTTTTTTTTAGGCCGGTTTTTGGTTAAGATAGCGGTGAGGCCGTTT‑TTTTTAAAAGCGGCGGt                                                                                                                                                       >  2:1140274/1‑150 (MQ=255)

aaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAATCCGGGTTATTCGTAAAGATAGCGAAGA‑GACGTTTTTTTATAAAAGCGGCGGt                                                                                                                                                       >  2:2180523/1‑150 (MQ=255)

     gACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCTCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCctc                                                                                                                                                  <  2:3371106/150‑1 (MQ=255)

     gACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTAGTAAATATGGCGATTA‑TCCGTTTTTTTTTAAAATCGGCGGttacaa                                                                                                                                                  >  2:3078907/1‑146 (MQ=255)

     gACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTGTTAAGCCGGGGTATTCGGAAAGAACGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTccac                                                                                                                                                  >  2:4338478/1‑148 (MQ=255)

     gACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCGGGCGTTTTGTTAAGCCGGGGGATGCGGAAAGATCGCGCTGA‑GCCCGTTGTGTAGAACACCCGCGGGTccac                                                                                                                                                  >  2:3945508/1‑148 (MQ=255)

       cTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCtt                                                                                                                                                <  2:1139332/150‑1 (MQ=255)

        tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCttt                                                                                                                                               <  2:2691602/150‑1 (MQ=255)

        tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCttt                                                                                                                                               <  2:2962559/150‑1 (MQ=255)

        tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGTTTTTTGTTTATACGGGGTATTGGTTAAGATAGCGGTGT‑GGCGTTTTTTTTTAAAAGCGGCGGGTGCGCttt                                                                                                                                               >  1:2169557/1‑150 (MQ=255)

            cagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTTGGTTTTTTTTTATCCGGGTTATTAGTAAAGATAGCGAGGA‑GCCGTTTTTTTTTAAAAGCGGCGGGTACGCTTTTggt                                                                                                                                           >  1:123234/1‑148 (MQ=255)

              gAACTGTCAGTCTAAGTCCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTc                                                                                                                                         <  1:1131515/150‑1 (MQ=255)

              gAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTc                                                                                                                                         <  1:2028375/150‑1 (MQ=255)

              gAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTc                                                                                                                                         <  1:3423502/150‑1 (MQ=255)

              gAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTc                                                                                                                                         <  1:1038355/150‑1 (MQ=255)

               aaCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTct                                                                                                                                        <  1:2377941/150‑2 (MQ=255)

                 cTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcgcg                                                                                                                                      <  2:4204228/150‑1 (MQ=255)

                  tGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCgcgc                                                                                                                                     <  2:433492/150‑1 (MQ=255)

                  tGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCgcgc                                                                                                                                     <  2:3239160/150‑1 (MQ=255)

                   gTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGcc                                                                                                                                    <  2:3904988/150‑1 (MQ=255)

                     cAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCa                                                                                                                                  <  2:1317375/150‑1 (MQ=255)

                     cAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCa                                                                                                                                  <  2:67817/150‑1 (MQ=255)

                        tCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGATTTTTTTTAATCCGGGTTATTAGTAAATATAGCGATTA‑GCCGGTTTTTTTTTAAATCGGCGGTTGCTCTTTTTCTTAGCGcccccga                                                                                                                               >  2:61808/1‑146 (MQ=255)

                             gTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgcg                                                                                                                          <  2:1525915/150‑1 (MQ=255)

                              tACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgcgc                                                                                                                         <  2:3742153/150‑1 (MQ=255)

                               aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTCAAGCTCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca                                                                                                                        <  2:3303665/150‑1 (MQ=255)

                                cGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTTTTAAACCCGGTTATTAGTAAAGATCGCGATGA‑GCCGGTTTTTTATAAAAACGGCGGGTAAGCATTTTCTTAGCGCCCACTAAGgcgccg                                                                                                                       >  2:3369337/1‑148 (MQ=255)

                                 gAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGGGGGTTTTTGTTAAGCCGGGTTATTCGTAAATGTCGCGCTGT‑TCCGGTTTTTTTTAAAAGCGGCGGTTGCTCTTTTTGGTCGCGCCCAGTCGGGGgccgc                                                                                                                      >  2:803739/1‑150 (MQ=255)

                                 gAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGGTTTTTTTTAAGCCGGGGTATTAGTAAATAAAGCGATGA‑GCCGTTTGTTTATAAAAGCGGCGGTTCCTCTTTTTCTTAGCGCCCCCTAGGGCgcagc                                                                                                                      >  2:2649528/1‑150 (MQ=255)

                                 gAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAATCCGGGTTATTAGGTAAGATCGCGATGA‑GCCGGTTGTTTATAAAAGCGGCGGTTCCTCTTTTGCGTCGCGCCCAGTCTGGGgccgc                                                                                                                      >  2:2425526/1‑150 (MQ=255)

                                 gAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTTTTAAGCCCGGGTATTCGGGAAGCTCGCGCTGG‑GCCCGGTGTTTATAACAGCCGCGGTTGCGCGTTTTGGTCGCGCCCCCTCGGgcgcccc                                                                                                                      >  2:1267051/1‑147 (MQ=255)

                                     aTTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACtt                                                                                                                  <  2:571988/150‑1 (MQ=255)

                                         ctactaCGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTgg                                                                                                              <  2:3323984/150‑1 (MQ=255)

                                            ctacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGa                                                                                                           <  1:3788663/150‑1 (MQ=255)

                                             tacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGAt                                                                                                          <  2:1987811/150‑1 (MQ=255)

                                                  cAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAATCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACATCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTGCTGGTGATTAATc                                                                                                     >  2:3992500/1‑150 (MQ=255)

                                                   aaTGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGAATTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑TCCGTTTTTTTATAACAGCGGCGGTTACTCTTTTGCTTAGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTcc                                                                                                    >  2:2868065/1‑150 (MQ=255)

                                                   aaTGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTGTTAATCCGGGTTATTAGTAAAGATCGCGATTA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCGTTTGCTTAGCGCCCAGTCAGGCGCAGCACTTCGGGTGATTaatac                                                                                                    >  2:129010/1‑146 (MQ=255)

                                                    aTGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCt                                                                                                   <  1:1730805/150‑1 (MQ=255)

                                                     tGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTg                                                                                                  <  1:366644/150‑1 (MQ=255)

                                                     tGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCACATTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTg                                                                                                  <  2:369635/150‑1 (MQ=255)

                                                      gTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTgg                                                                                                 <  2:2773216/150‑1 (MQ=255)

                                                         gAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGgcgc                                                                                              <  1:898429/150‑1 (MQ=255)

                                                          aaCAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAAACCGGGTTATTAGTAAACAACGCCAAGA‑GCCGTTTTTTTATAAAAGCGGCGGTTACACTTTTTATTAGCGCCCAGTCGGGCGCACCACAACTGGTGATTACTACAGgcgcg                                                                                             >  2:3793632/1‑149 (MQ=255)

                                                            cAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTcc                                                                                           <  2:3266452/150‑1 (MQ=255)

                                                                   ggATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                                                                                    <  2:2066251/150‑1 (MQ=255)

                                                                       cccGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaa                                                                                <  2:2842890/150‑1 (MQ=255)

                                                                             cGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga                                                                          <  1:991865/150‑1 (MQ=255)

                                                                             cGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga                                                                          <  2:3521847/150‑1 (MQ=255)

                                                                              gTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAagag                                                                         <  1:1066548/150‑1 (MQ=255)

                                                                              gTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAagag                                                                         <  1:62421/150‑1 (MQ=255)

                                                                              gTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAagag                                                                         >  1:1397419/1‑150 (MQ=255)

                                                                                ttCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTc                                                                       <  2:732133/150‑1 (MQ=255)

                                                                                   tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTa                                                                    >  1:3518049/1‑150 (MQ=255)

                                                                                    ggATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTAc                                                                   >  2:1742426/1‑150 (MQ=255)

                                                                                      aTTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActc                                                                 >  2:330932/1‑150 (MQ=255)

                                                                                      aTTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActc                                                                 >  2:2448879/1‑150 (MQ=255)

                                                                                        ttttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCtt                                                               <  2:1899694/150‑1 (MQ=255)

                                                                                          ttGTGAAGCCGGGTTAGTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTGTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTAAGACTCtttt                                                             <  1:3381477/150‑3 (MQ=255)

                                                                                           tGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCgg                                                            <  2:3890028/150‑1 (MQ=255)

                                                                                             ttAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggc                                                          >  1:3908267/1‑150 (MQ=255)

                                                                                             ttAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACATCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggc                                                          >  1:3027435/1‑150 (MQ=255)

                                                                                             ttAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTATTTTTCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTAGggc                                                          >  1:1683342/1‑150 (MQ=255)

                                                                                               aaGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcgg                                                        <  2:1355425/150‑1 (MQ=255)

                                                                                                 gCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCt                                                      <  1:1883613/150‑1 (MQ=255)

                                                                                                  ccGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCtt                                                     <  2:1071577/150‑1 (MQ=255)

                                                                                                  ccGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCtt                                                     <  2:411017/150‑1 (MQ=255)

                                                                                                     ggTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTac                                                  <  2:67175/150‑1 (MQ=255)

                                                                                                      gTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTacc                                                 <  2:3026589/150‑1 (MQ=255)

                                                                                                       ttATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTacca                                                >  2:2292097/1‑150 (MQ=255)

                                                                                                       ttATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCACTTTTGCTTCGCGCCCAGGCGGGCGCAGCACGTCGGGTGATTACTCCTGGCGCTCCGGTTCGGGCCAAAAAGATTCTTAATCATCGGGCGGCGTTacca                                                >  1:3461391/1‑150 (MQ=255)

                                                                                                        tATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTaccac                                               <  1:1553433/150‑1 (MQ=255)

                                                                                                         aTTCGTAAAGATCGCGATGA‑GCCGTTTTTTTATAAAAGCGGCGGTTCCACTTTTGCTTAGCGCCCAGTCTTGCGCAGCACTTCTGGTGAATAATCCTGGCGCGACGGTTCTGGCAAAAAAGAGACTTAATCGTCGGGCGGCTTTaccacc                                              >  1:1941512/1‑150 (MQ=255)

                                                                                                          ttCGTAAAGATCGCGCTGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCAGAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAATAGAGTCTTACTCTTCGGGCGGCTTTAccacca                                             <  1:3147674/150‑1 (MQ=255)

                                                                                                          ttCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCGTCGCGCCCAGTCTGGCGCACCACATCTGGTGATTACTCCAGGCGCTCCCGTTCTGGCAAAAAAGAGTCGTACTCGTCCGGCGGCGTTaccaccc                                             >  1:2724442/1‑149 (MQ=255)

                                                                                                          ttCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCACTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACACCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACACTTCGGGCGGCTTTAccacca                                             >  1:4265729/1‑150 (MQ=255)

                                                                                                             gTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCaaaa                                          >  2:3545908/1‑150 (MQ=255)

                                                                                                             gTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTACGCTTTTTCTTCGCGCCCAGTCTGGCGCAGCACATCAGGTGATTACACCAGGCGCTCCGGTTCTGGCCAAAAAGAGTCATACACTTAGGGCGGCTTTACCACCaaaa                                          >  1:3983107/1‑150 (MQ=255)

                                                                                                              tAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTTCTTCGCGCCCAGTCTGGCGCAGCACTTCAGGTGATTACTCCAGGCACTCCGGTTCTGACAAAAAATAGTCTTAAACGTCGGGCGGCTTTACCCCCaaaaa                                         >  1:2538816/1‑150 (MQ=255)

                                                                                                              tAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCaaaaa                                         >  2:3825363/1‑150 (MQ=255)

                                                                                                               aaaGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACGTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCaaaaaa                                        >  1:2850107/1‑150 (MQ=255)

                                                                                                                aaGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAg                                       >  2:4214838/1‑150 (MQ=255)

                                                                                                                 aGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGc                                      <  2:1824429/150‑1 (MQ=255)

                                                                                                                 aGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGc                                      <  2:3977929/150‑1 (MQ=255)

                                                                                                                 aGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGc                                      >  2:2574199/1‑150 (MQ=255)

                                                                                                                 aGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGc                                      >  2:3495264/1‑150 (MQ=255)

                                                                                                                      gcgTTGC‑GCCGTTTTTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCTGTCTGGCGCAGCACTTCTGGTGTTTACTCCTGGCGCTCCGGTTCTGGAAAAAATGATTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTa                                 <  1:4236335/150‑1 (MQ=255)

                                                                                                                        gATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCACCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCCGCTTTACCACCAAAAAAGCGACTAcc                               >  2:4045367/1‑150 (MQ=255)

                                                                                                                          tGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGt                             <  2:440237/150‑1 (MQ=255)

                                                                                                                              gCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTtc                          >  2:89501/1‑150 (MQ=255)

                                                                                                                                  tttgtttATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTcgc                      <  1:132668/150‑1 (MQ=255)

                                                                                                                                           aaCAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTAt             <  2:718507/150‑1 (MQ=255)

                                                                                                                                             cAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATgg           <  2:1065297/150‑1 (MQ=255)

                                                                                                                                               gcggcgGGTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTc         <  1:3378368/150‑1 (MQ=255)

                                                                                                                                                cggcggTTCCTCTTTTTCTTCGCGCCCAGTCTGTCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCTTTTCTGTCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCg        <  1:549588/150‑1 (MQ=255)

                                                                                                                                                  gcggTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGtt      <  2:222007/150‑1 (MQ=255)

                                                                                                                                                    ggTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGttct    <  1:2284376/150‑1 (MQ=255)

                                                                                                                                                    ggTTCCTCTTTTGCTTCACGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGttct    <  1:73283/150‑1 (MQ=255)

                                                                                                                                                      ttCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGttcttc  <  2:3197254/150‑1 (MQ=255)

AATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGA‑GCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTC  >  NC_000913/4412699‑4412997

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