BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,118,623	C→A	6.0%	G142C (GGT→TGT)	yceJ ←	putative cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,118,623	0	C	A	6.0%	319.9 / 5.9	107	G142C (GGT→TGT)	yceJ	putative cytochrome b561
Reads supporting (aligned to +/- strand): ref base C (43/57); new base A (3/4); total (46/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

CATTTTTTATACTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAG  >  NC_000913/1118474‑1118772

cATTTTTTATACTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCAcc                                                                                                                                                       <  1:1431042/150‑1 (MQ=255)

 aTTTTTTATACTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTGTTTAACGATGAAATGATGCTTCAGGGCCATAAAAACGAGCATAACGGACAGTACTACCACCCTCCACGCGAGCCAAAAATGCAGGGCACCg                                                                                                                                                      >  1:899417/1‑150 (MQ=255)

   tttttAAACTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTCTCGTCTTTATCGATGATATGATGCTTCAGGTCCATAAATCCGTGCATAACGGACAGTACTACGACCCTCAACGCTAGCCAAAAATCCAGGGCACCgcc                                                                                                                                                    >  2:311566/1‑148 (MQ=255)

    ttttATACTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCa                                                                                                                                                   <  1:667919/150‑1 (MQ=255)

    ttttATACTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCGTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATTCAGGGCACCGGCa                                                                                                                                                   >  2:331189/1‑150 (MQ=255)

       tatcCTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAg                                                                                                                                                <  2:675455/150‑1 (MQ=255)

        ataCTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGt                                                                                                                                               >  1:626184/1‑150 (MQ=255)

        ataCTCCATAGGCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTAGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGt                                                                                                                                               <  2:1417788/150‑1 (MQ=255)

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

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

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

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

             ccATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGcc                                                                                                                                          <  1:1600241/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)

                       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)

                                 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)

                                                      gtcgtcTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCCAACCa                                                                                                 >  1:8640/1‑150 (MQ=255)

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

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

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

                                                                 cGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGGCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGc                                                                                      <  2:1174220/150‑1 (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                                                                            <  1:1613983/150‑1 (MQ=255)

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

                                                                           atgCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGttt                                                                            <  2:408283/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)

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

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

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

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

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

                                                                                                   cATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATa                                                    >  2:665181/1‑150 (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)

                                                                                                          gACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTg                                             <  2:16624/150‑1 (MQ=255)

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

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

                                                                                                              gTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTGTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGc                                         <  1:331189/150‑1 (MQ=255)

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

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

                                                                                                               tactacGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGcc                                        >  1:1519372/1‑150 (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)

                                                                                                                                    ccAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTa                   >  2:50844/1‑150 (MQ=255)

                                                                                                                                     cAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTAc                  >  1:1631811/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                            gCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCga           <  1:737128/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)

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

CATTTTTTATACTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAG  >  NC_000913/1118474‑1118772

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