BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,400,330	C→G	6.9%	L15V (CTG→GTG)	hfq →	RNA‑binding protein Hfq

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,400,330	0	C	G	6.9%	367.2 / 12.8	113	L15V (CTG→GTG)	hfq	RNA‑binding protein Hfq
Reads supporting (aligned to +/- strand): ref base C (58/47); new base G (6/2); total (64/49)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 115 total aligned reads displayed.

TGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCC  >  NC_000913/4400181‑4400477

tgttgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCAc                                                                                                                                                      <  2:350597/150‑1 (MQ=255)

tgttgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCAc                                                                                                                                                      <  1:137829/150‑1 (MQ=255)

  ttgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAATGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTg                                                                                                                                                    <  2:1440169/150‑1 (MQ=255)

  ttgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTg                                                                                                                                                    <  1:199238/150‑1 (MQ=255)

      tGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑Tc                                                                                                                                               >  2:92367/1‑150 (MQ=255)

       gCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCg                                                                                                                                              >  1:301122/1‑150 (MQ=255)

          aTCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGa                                                                                                                                           >  2:380443/1‑150 (MQ=255)

             gCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACg                                                                                                                                        <  2:943942/150‑1 (MQ=255)

             gCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACg                                                                                                                                        <  2:1012384/150‑1 (MQ=255)

             gCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCCG‑ACGCAGTGCGTTCAGGAACg                                                                                                                                        <  1:620907/150‑1 (MQ=255)

                ggCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCCG‑ACGCAGTGCGTTCAGGAACggat                                                                                                                                     <  1:128550/150‑4 (MQ=255)

                   tGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACCTGTTcc                                                                                                                                  >  2:806274/1‑150 (MQ=255)

                    gAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCa                                                                                                                                 >  1:104645/1‑150 (MQ=255)

                         tgtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGtttg                                                                                                                            >  2:113502/1‑149 (MQ=255)

                         tgtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTc                                                                                                                            <  1:437108/150‑1 (MQ=255)

                          gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCggg                                                                                                                                            <  1:1418365/133‑1 (MQ=255)

                          gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCggg                                                                                                                                            >  2:1418365/1‑133 (MQ=255)

                          gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCt                                                                                                                           >  2:738802/1‑150 (MQ=255)

                          gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCt                                                                                                                           <  2:902738/150‑1 (MQ=255)

                          gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGCGTTCCAGTTTCt                                                                                                                           <  2:771765/150‑1 (MQ=255)

                             cAATTGAGACGTATCGTGCGCAATTTTTTCAGATTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCtatt                                                                                                                        <  1:320298/150‑1 (MQ=255)

                               aTTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCtattta                                                                                                                      <  1:48028/150‑1 (MQ=255)

                                 tGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCtatttatt                                                                                                                    <  1:686575/150‑1 (MQ=255)

                                  gagaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTatttattt                                                                                                                   >  2:992667/1‑150 (MQ=255)

                                    gaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTgg                                                                                                                 <  2:1171554/150‑1 (MQ=255)

                                     aCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGt                                                                                                                <  2:1187161/150‑1 (MQ=255)

                                      cGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTg                                                                                                               <  2:1512561/150‑1 (MQ=255)

                                       gTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGa                                                                                                              <  1:286302/150‑1 (MQ=255)

                                         aTCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAAt                                                                                                            >  2:1288804/1‑150 (MQ=255)

                                             tGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTa                                                                                                        <  1:380443/150‑1 (MQ=255)

                                             tGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTa                                                                                                        <  1:315106/150‑1 (MQ=255)

                                              gcgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTAt                                                                                                       <  2:359593/150‑1 (MQ=255)

                                                gcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTa                                                                                                     >  2:776469/1‑150 (MQ=255)

                                                   aTTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGc                                                                                                  >  2:648667/1‑150 (MQ=255)

                                                    ttttttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCt                                                                                                 <  2:91470/150‑1 (MQ=255)

                                                      ttttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGc                                                                                               >  1:886000/1‑150 (MQ=255)

                                                       tttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCa                                                                                              >  1:729177/1‑150 (MQ=255)

                                                        ttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCaa                                                                                             >  1:1068211/1‑150 (MQ=255)

                                                         tCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAg                                                                                            <  2:289557/150‑1 (MQ=255)

                                                          cAGAATCCAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAAAGTGTTACAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAgg                                                                                           >  2:603108/1‑150 (MQ=255)

                                                           aGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAggg                                                                                          >  2:914633/1‑150 (MQ=255)

                                                           aGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAggg                                                                                          <  1:738802/150‑1 (MQ=255)

                                                            gAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGc                                                                                         <  1:309084/150‑1 (MQ=255)

                                                              aTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCaa                                                                                       >  1:594444/1‑150 (MQ=255)

                                                                cGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAAt                                                                                     <  2:675327/150‑1 (MQ=255)

                                                                      gTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTc                                                                               <  2:829581/150‑1 (MQ=255)

                                                                         cAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTACTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCttt                                                                            >  1:809496/1‑150 (MQ=255)

                                                                          aaaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCtttt                                                                           >  1:316346/1‑150 (MQ=255)

                                                                             gTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGAt                                                                        >  1:57319/1‑150 (MQ=255)

                                                                              tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                       >  1:622430/1‑150 (MQ=255)

                                                                              tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                       >  1:1098300/1‑150 (MQ=255)

                                                                              tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                       >  1:457141/1‑150 (MQ=255)

                                                                              tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                       <  1:913416/150‑1 (MQ=255)

                                                                                cAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAg                                                                     >  2:294338/1‑150 (MQ=255)

                                                                                  aaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGtt                                                                   >  2:931401/1‑150 (MQ=255)

                                                                                  aaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGtt                                                                   >  2:1080882/1‑150 (MQ=255)

                                                                                   aTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTc                                                                  >  1:1041739/1‑150 (MQ=255)

                                                                                    tAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCg                                                                 <  2:1423049/150‑1 (MQ=255)

                                                                                     aaGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGt                                                                >  2:695169/1‑150 (MQ=255)

                                                                                     aaGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGt                                                                <  2:1204170/150‑1 (MQ=255)

                                                                                      aGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTg                                                               <  2:104645/150‑1 (MQ=255)

                                                                                        cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTTTTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGAt                                                             <  1:603108/150‑1 (MQ=255)

                                                                                        cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGAt                                                             <  1:776469/150‑1 (MQ=255)

                                                                                        cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGAt                                                             >  1:150227/1‑150 (MQ=255)

                                                                                          tataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATcc                                                           <  2:886000/150‑1 (MQ=255)

                                                                                          tataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATcc                                                           >  1:250099/1‑150 (MQ=255)

                                                                                           ataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCt                                                          >  1:1121513/1‑150 (MQ=255)

                                                                                            taAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTg                                                         >  1:889960/1‑150 (MQ=255)

                                                                                                 aaaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGaa                                                    >  2:1461701/1‑150 (MQ=255)

                                                                                                    agagagAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTg                                                      >  1:538318/1‑145 (MQ=255)

                                                                                                    agagagAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTg                                                      <  2:538318/145‑1 (MQ=255)

                                                                                                     gagagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAac                                                >  1:366383/1‑150 (MQ=255)

                                                                                                       gagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAacac                                              <  1:650348/150‑1 (MQ=255)

                                                                                                         gaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACgg                                            <  1:695169/150‑1 (MQ=255)

                                                                                                           aTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACAGTc                                          <  1:806274/150‑1 (MQ=255)

                                                                                                              gCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGc                                       <  2:1090320/150‑1 (MQ=255)

                                                                                                                   ggggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGAt                                  >  1:998482/1‑150 (MQ=255)

                                                                                                                    gggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATg                                 >  1:1312988/1‑150 (MQ=255)

                                                                                                                    aacggaTCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATg                                 >  1:879795/6‑150 (MQ=255)

                                                                                                                     ggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATgg                                >  1:849339/1‑150 (MQ=255)

                                                                                                                     ggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATgg                                >  1:426514/1‑150 (MQ=255)

                                                                                                                        aaTCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGttt                             >  1:644735/1‑150 (MQ=255)

                                                                                                                           cTTTACAAGACCCGTTCCTGAACGCACTGCG‑TCTGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                          >  2:256172/1‑150 (MQ=255)

                                                                                                                             ttACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAg                        <  1:1288804/150‑1 (MQ=255)

                                                                                                                               aCAAGATCCGTTGCTGAACGCACTGGG‑TCGGGAAGGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCATGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCa                      <  2:1121513/150‑1 (MQ=255)

                                                                                                                               aCAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCa                      >  2:1518132/1‑150 (MQ=255)

                                                                                                                               aCAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCa                      <  1:1461701/150‑1 (MQ=255)

                                                                                                                                 ggaacacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcg                    >  1:833495/7‑150 (MQ=255)

                                                                                                                                 ggaacacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcg                    >  1:1242568/7‑150 (MQ=255)

                                                                                                                                   aacacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgcg                  >  2:392598/5‑150 (MQ=255)

                                                                                                                                    acacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGa                 >  2:1216072/4‑150 (MQ=255)

                                                                                                                                    acacGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGa                 >  1:800277/4‑150 (MQ=255)

                                                                                                                                      acGTTCCCG‑ACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAtt               >  2:58710/2‑150 (MQ=255)

                                                                                                                                       cGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAttt              >  1:241350/1‑150 (MQ=255)

                                                                                                                                          tCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTa           >  2:850014/1‑150 (MQ=255)

                                                                                                                                           ccTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc          >  2:1344937/1‑150 (MQ=255)

                                                                                                                                           ccTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc          <  1:44645/150‑1 (MQ=255)

                                                                                                                                            cTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACt         <  1:1344937/150‑1 (MQ=255)

                                                                                                                                              gAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtgt       <  1:1124517/150‑1 (MQ=255)

                                                                                                                                                   cACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTcc  >  1:1141723/1‑150 (MQ=255)

TGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCG‑TCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCC  >  NC_000913/4400181‑4400477

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