BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,118,630	G→C	5.6%	D139E (GAC→GAG) ‡	yceJ ←	putative cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,118,630	0	G	C	5.6%	357.0 / 10.1	107	D139E (GAC→GAG) ‡	yceJ	putative cytochrome b561
Reads supporting (aligned to +/- strand): ref base G (42/59); new base C (3/3); total (45/62)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

Only 100 of 109 total aligned reads displayed.

TATACTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGC  >  NC_000913/1118481‑1118779

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

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

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

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

    cTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTg                                                                                                                                                   >  1:304298/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)

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

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

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

               gATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACAGGCAGACTTTGCCGGTgccctg                                                                                                                                        <  1:1339413/150‑5 (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:70412/150‑1 (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                                                                                   >  2:499648/1‑150 (MQ=255)

                                                                    atgCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGttt                                                                                   <  1:1613983/150‑1 (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)

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

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

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

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

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

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

                                                                                                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:304298/150‑1 (MQ=255)

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

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

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

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

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

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

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

TATACTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTCTTGC  >  NC_000913/1118481‑1118779

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