BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,118,632	C→T	5.5%	D139N (GAC→AAC) ‡	yceJ ←	putative cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,118,632	0	C	T	5.5%	344.3 / 6.2	108	D139N (GAC→AAC) ‡	yceJ	putative cytochrome b561
Reads supporting (aligned to +/- strand): ref base C (42/60); new base T (3/3); total (45/63)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 109 total aligned reads displayed.

CTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGG  >  NC_000913/1118485‑1118781

cTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTg                                                                                                                                                     >  1:304298/1‑150 (MQ=255)

cTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTg                                                                                                                                                     <  2:1066267/150‑1 (MQ=255)

  ccATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGcc                                                                                                                                                   <  1:1600241/150‑1 (MQ=255)

  ccATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGcc                                                                                                                                                   <  2:117094/150‑1 (MQ=255)

   cATAGTCAGATGAGGACTTTCCCAGCATGCGCTTCAGAGTGTCGGCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCt                                                                                                                                                  <  1:1209508/150‑1 (MQ=255)

   cATAGTAAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCt                                                                                                                                                  <  2:1255041/150‑1 (MQ=255)

       ctcccacgacgacTTTCCCACCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCGTAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACAGGCAGACTTTGCCGGTgc                                                                                                                                              <  2:1665673/143‑1 (MQ=255)

        tCAGATGACGACTTTCCCAGCTTGCGCTTCAGAGTGGCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTgcg                                                                                                                                             <  1:327213/150‑1 (MQ=255)

         cAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACAGGCAGACTTTGCCGGTgccc                                                                                                                                            <  1:532774/150‑3 (MQ=255)

         cAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACCCTCCCCGCGAGCCAAAAATGCCGGGCACCGGCAAAGTCTGCCTGTgcgc                                                                                                                                            >  1:516018/1‑150 (MQ=255)

            aTGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAg                                                                                                                                         >  2:802781/1‑150 (MQ=255)

             tGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCCGGGCACCGGCAAAGTCTGCCTGTGCGCCAGc                                                                                                                                        >  1:1375482/1‑150 (MQ=255)

                   cTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTGCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGc                                                                                                                                  <  1:595601/150‑1 (MQ=255)

                   atttCCGAGCATGCGCTTCAGAGTTTCGTCATTAGAGATGAAATGATGCTTCAGGGCCATAAATCCGTTCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAATATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGc                                                                                                                                  <  1:1785887/149‑1 (MQ=255)

                      tCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAg                                                                                                                               <  1:1372449/150‑1 (MQ=255)

                      tCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAg                                                                                                                               <  1:70412/150‑1 (MQ=255)

                      tCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGAAGAAATTATGCTACAGGGCCATAAATACGTGCATAACGGACAGTACTACGACCCTCCACGCGAGCCAAAAATGCAGGTCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAg                                                                                                                               >  1:876044/1‑150 (MQ=255)

                      tCCCAGCATGCGCTTCAGAGTGTCGTCTTAATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATACCGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAg                                                                                                                               <  1:1167620/150‑1 (MQ=255)

                       cccAGCATGCGCTTCAGAGTGGCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAgg                                                                                                                              <  1:1025256/150‑1 (MQ=255)

                        ccAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGAACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAggg                                                                                                                             <  1:868167/150‑1 (MQ=255)

                          aGCATGCTCTTCAGAGTGTCGTCTTTAACGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTc                                                                                                                           <  1:187492/150‑1 (MQ=255)

                          aGCATGAGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTc                                                                                                                           <  1:1737998/150‑1 (MQ=255)

                              tGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCgg                                                                                                                       <  1:406260/150‑1 (MQ=255)

                               gcgcTTCAGAGTGTCGTCTTTATCGTTGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCggg                                                                                                                      <  1:1107031/150‑1 (MQ=255)

                                  cTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGAc                                                                                                                   >  2:737128/1‑150 (MQ=255)

                                     cAGAGTGTCGTCTTTATCGATGAAATGATGCTGCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTc                                                                                                                <  1:537241/150‑1 (MQ=255)

                                      agagTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATACAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCa                                                                                                               <  2:176113/150‑1 (MQ=255)

                                      agagTGTCGTCTTTATCGATGAAATGATGATTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCa                                                                                                               >  1:1289983/1‑150 (MQ=255)

                                          tgtCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAAcc                                                                                                           <  2:1825223/150‑1 (MQ=255)

                                                    aTCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                                                                                                 <  2:899417/150‑1 (MQ=255)

                                                    aTCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGAACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                                                                                                 <  1:504091/150‑1 (MQ=255)

                                                      cGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTATGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGc                                                                                               >  2:1078172/1‑150 (MQ=255)

                                                       gATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCt                                                                                              >  1:627934/1‑150 (MQ=255)

                                                        aTGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTg                                                                                             <  2:227059/150‑1 (MQ=255)

                                                        aTGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGCCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTg                                                                                             <  2:896513/150‑1 (MQ=255)

                                                        aTGAAATGATGCTTCAGGGCCATAAATCCGAGCATACCGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTg                                                                                             <  1:994925/150‑1 (MQ=255)

                                                         tGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGCGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGACAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGa                                                                                            >  2:1653294/1‑150 (MQ=255)

                                                         tGAAATGATGCTTCAGGGCCATAAATCAGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCCAAAACGCTGa                                                                                            >  1:892046/1‑150 (MQ=255)

                                                              tgatgCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGt                                                                                       >  2:790865/1‑150 (MQ=255)

                                                                atgCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGttt                                                                                     <  2:408283/150‑1 (MQ=255)

                                                                atgCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGttt                                                                                     >  2:499648/1‑150 (MQ=255)

                                                                atgCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGttt                                                                                     <  1:1613983/150‑1 (MQ=255)

                                                                   cTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTAcc                                                                                  >  2:845728/1‑150 (MQ=255)

                                                                   cTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTATAcc                                                                                  >  1:421293/1‑150 (MQ=255)

                                                                         ggCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGc                                                                            <  2:1294785/150‑1 (MQ=255)

                                                                          gCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGCCGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCa                                                                           >  1:1417209/1‑150 (MQ=255)

                                                                             aTAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCCGCCAAAAACGCTGCTCGGTTTACCATCGGCAGtt                                                                        >  1:436274/1‑150 (MQ=255)

                                                                              tAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTg                                                                       <  2:626184/150‑1 (MQ=255)

                                                                              tAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACACTGATCGGTTTACCATCGGCAGTTg                                                                       <  2:1354108/150‑1 (MQ=255)

                                                                               aaaTCCGTGCATACCGGACAGTAGTACGCCGCTCCACGAGAGCCAAAAATTCAGGGCACCTGCAAAGTCTGCCGGTGCGCCAGCGGCGGCAAGGGTCGCGGTGACGTCAAACCAGCAAAAAACGCTGATCGGTTTACCATCGGCAGTTga                                                                      <  2:287301/150‑1 (MQ=255)

                                                                                 aTCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTgaga                                                                    <  1:396689/150‑1 (MQ=255)

                                                                                   ccGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATc                                                                  >  1:30143/1‑150 (MQ=255)

                                                                                    cGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCa                                                                 >  2:1343045/1‑150 (MQ=255)

                                                                                    cGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCa                                                                 <  2:1065733/150‑1 (MQ=255)

                                                                                     gTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCAGCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAg                                                                <  2:1226913/150‑1 (MQ=255)

                                                                                        cATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATa                                                             >  2:665181/1‑150 (MQ=255)

                                                                                        cATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATa                                                             <  1:1517038/150‑1 (MQ=255)

                                                                                          tAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGc                                                           >  1:686630/1‑150 (MQ=255)

                                                                                            aCGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGCCGTCACACCAGCCAAAACCGCTGATCGGTTTACCATCGGCAGTTGAGCTCAGATAGCCg                                                         >  1:1640517/1‑150 (MQ=255)

                                                                                            aCGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTTATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCg                                                         <  2:1651106/150‑1 (MQ=255)

                                                                                            aCGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCg                                                         <  2:1117977/150‑1 (MQ=255)

                                                                                              ggACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCt                                                       <  1:1270143/150‑1 (MQ=255)

                                                                                                 cAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGAACCGGCATAGTCTGCCTGTGCGCCAGCGTCGGGTAGGGTCGCGGGGACGTCAAACCAGCCAAATACGCTGATCGGTTTACCGTCGGCAGTTGAGATCAGAAAGCCTCTGAt                                                    >  1:254720/1‑150 (MQ=255)

                                                                                                  aGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATg                                                   <  1:802781/150‑1 (MQ=255)

                                                                                                    tactacGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGcc                                                 >  1:1519372/1‑150 (MQ=255)

                                                                                                    tactacGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGcc                                                 <  2:308071/150‑1 (MQ=255)

                                                                                                    tactacGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGcc                                                 <  2:304298/150‑1 (MQ=255)

                                                                                                     actacGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCg                                                >  1:1200377/1‑150 (MQ=255)

                                                                                                      ctacGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGa                                               >  1:1403253/1‑150 (MQ=255)

                                                                                                          gacgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGc                                           <  2:516018/150‑1 (MQ=255)

                                                                                                          gacgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGc                                           <  1:498941/150‑1 (MQ=255)

                                                                                                            cgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGa                                         <  1:335998/150‑1 (MQ=255)

                                                                                                                 cACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGc                                    >  1:675430/1‑150 (MQ=255)

                                                                                                                  aCGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGGGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCa                                   >  2:378599/1‑150 (MQ=255)

                                                                                                                   cgcgAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCAGCGGCAGTTGCGATCAGATAGCCGCAGATGCCGATGGCGAAAAGaaa                                  >  1:1618626/1‑150 (MQ=255)

                                                                                                                       aGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAgg                              <  2:876044/150‑1 (MQ=255)

                                                                                                                          cAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCAACGGCAGTTGAGATCAGATATCCGCTGAAGCCGATGGCGAAAAGCAATACGTAc                           >  2:1106733/1‑150 (MQ=255)

                                                                                                                             aaaTGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAAt                        >  1:464470/1‑150 (MQ=255)

                                                                                                                                tGCAGGGCACCGGCAAAGTCTGCCTGTGCGCAAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCg                     <  1:64686/150‑1 (MQ=255)

                                                                                                                                 gCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCga                    <  1:737128/150‑1 (MQ=255)

                                                                                                                                 gCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCga                    >  1:1618158/1‑150 (MQ=255)

                                                                                                                                 gCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCga                    <  2:130216/150‑1 (MQ=255)

                                                                                                                                  ctgggCACCGGCCAAGTCTGCCTGGGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCgag                   <  2:326131/148‑1 (MQ=255)

                                                                                                                                  cAGGGCACCGGCAAAGTCTGCTTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCgag                   <  2:1507693/150‑1 (MQ=255)

                                                                                                                                  cAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCgag                   >  1:1543640/1‑150 (MQ=255)

                                                                                                                                   ggcgCACAGGCAGACTTTGCCGGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCATACCAGCCAAACACGATGATCGGTTTACCATCGGCAGGTTAGCTCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATTCgagg                  >  1:244298/4‑149 (MQ=255)

                                                                                                                                   agtgCACAGGAAAAGTCTGCCTGTGAGTCAGAGTGGGCAATTGGCGCGTGGAAGTCAAACCATCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATTCCGATGGCGAAAAGCAATAGGTACAATGCgaga                  <  1:1292344/147‑1 (MQ=255)

                                                                                                                                     ggCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATg                >  2:634778/1‑150 (MQ=255)

                                                                                                                                     cgCACAGGCAGACTTTGCCGGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAGCGCTGATCGGTTTACCATCGGCAGGTGATATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATg                >  1:726667/2‑150 (MQ=255)

                                                                                                                                      gCACAGGCAGACTTTGCCGGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGt               >  1:783820/1‑150 (MQ=255)

                                                                                                                                        aCCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTcc             >  1:1006149/1‑150 (MQ=255)

                                                                                                                                          cGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAACCCAGCCAAAAACGCTGCTCGGTTTACCCTCGGCAGTTGAGATCAGATAGCCGCGGATGCCGATGGCGAAAAGCAAGCGGTCCCATGCGAGATGTcccg           >  1:319164/1‑148 (MQ=255)

                                                                                                                                          cGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAg           >  1:1028162/1‑150 (MQ=255)

                                                                                                                                           ggCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGc          <  2:1052522/150‑1 (MQ=255)

                                                                                                                                           ggCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGc          <  2:1375482/150‑1 (MQ=255)

                                                                                                                                             cAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACAAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCAGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGctc        <  2:1289983/150‑1 (MQ=255)

                                                                                                                                                 gTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGCGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGc    >  1:1023640/1‑150 (MQ=255)

                                                                                                                                                  tCTGCCTATGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCg   <  1:1078172/150‑1 (MQ=255)

                                                                                                                                                  tCTGCCGGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCg   <  1:1675149/150‑1 (MQ=255)

                                                                                                                                                   cTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCgg  <  1:69143/150‑1 (MQ=255)

CTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGCGG  >  NC_000913/1118485‑1118781

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